本人写了一个程序，目的主要是用一个函数的数值积分去拟合实际数据，确定函数中的参数。现在程序中出现了两个关键问题，求各位大神帮忙~程序如下:



如果最后一张图不清楚的话可以见附件：
主要问题如下：
1、主程序中显示，Function value and YDATA sizes are not equal.
这里xdata是一个400*2的矩阵，P是一个400*1的列向量。
2、关于omega的处理，最内层的函数firststage的值是和omega有关的，curvefun函数是对omega的一个广义数值积分，不知道这种情况下omega应该怎样定义和处理才能使值能够正常的传递，使输入x和xdata后curvefun可以得出结果呢？我现在没有对omega进行任何定义，因为能想到的处理方法都不对...
谢谢大家！