p * (1-0.5*(n*p)^2 + 1/12 * (n*p)^4 - 1/144*(n*p)^6)*(1-0.5*(p* 2 * pi * cos(f)*sin(123)/630*10^9)^2 + (2 * pi * cos(f)*sin(123)/630*10^9 * p)^4/24     - (2 * pi * cos(f)*sin(123)/630*10^9*p)^6/720)

这个函数，关于p 和 f 双重积分，咋搞？

dblquad(@(p,f)p * (1-0.5*(n*p)^2 + 1/12 * (n*p)^4 - 1/144*(n*p)^6)*(1-0.5*(p* 2 * pi * cos(f)*sin(123)/630*10^9)^2 + (2 * pi * cos(f)*sin(123)/630*10^9 * p)^4/24     - (2 * pi * cos(f)*sin(123)/630*10^9*p)^6/720) ,0.01,0.19,0,2*pi  )
    

这样会报错：

??? Error using ==> mpower
Matrix must be square.

Error in ==>
@(p,f)p*(1-0.5*(n*p)^2+1/12*(n*p)^4-1/144*(n*p)^6)*(1-0.5*(p*2*pi*cos(f)*sin(123)/630*10^9)^2+(2*pi*cos(f)*sin(123)/630*10^9*p)^4/24-(2*pi*cos(f)*sin(123)/630*10^9*p)^6/720)


Error in ==> quad at 77
y = f(x, varargin{:});

Error in ==> dblquad>innerintegral at 84
    Q(i) = quadf(intfcn, xmin, xmax, tol, trace, y(i), varargin{:});

Error in ==> quad at 77
y = f(x, varargin{:});

Error in ==> dblquad at 60
Q = quadf(@innerintegral, ymin, ymax, tol, trace, intfcn, ...

Error in ==> wang at 3
dblquad(@(p,f)p * (1-0.5*(n*p)^2 + 1/12 * (n*p)^4 - 1/144*(n*p)^6)*(1-0.5*(p* 2 * pi *
cos(f)*sin(123)/630*10^9)^2 + (2 * pi * cos(f)*sin(123)/630*10^9 * p)^4/24     - (2 *
pi * cos(f)*sin(123)/630*1

谢谢指教