Question

求matlab大神！！！！matlab计算二重积分！！！

计算D的值。下面是我写的程序，帮忙看下哪里有问题，
syms x yu=(pi/2)*sin(x)*cos(y);v=(cos(u))^2;v_1=(sin(x))^2*(cos(y))^2;s=(pi/2)*sin(x)*sin(y);k=(cos(s))^2;k_1=(sin(x))^2*(sin(y))^2;q=(sin(x))^2;q_1=2*sin(y)*cos(y);w=cos(u);w_1=cos(s);f=(v/(1-v_1)+k/(1-k_1)-((q)*(q_1)*(w)*(w_1))/((1-v_1)*(1-k_1)))*sin(x);D1=8*pi/int(int(f,x,0,pi),y,0,2*pi)D=vpa(D1,9)
我的结果是1.64，老师说不对，正确结果也大概是这个范围，大于1.38.

Answer 1

怎么说呢，MATLAB中双重积分不是这么解的，用dblquad这个函数才能得到像样的数值解
调用格式dblquad(FUN,Xm,XM,Ym,YM)%意义依次是函数，X下限，上限，Y下限，上限

追问

能否给出完整的格式呢谢谢

追答

就用f=@(x,y)...把这个被积函数打出来，这个太麻烦我就不打了，然后就按照我给出的格式调用就行了

追问

我觉得就是f这里面出问题了 但是找不到

追答

对，忘了说了，这个F函数的格式是F=@（变量列表）表达式。先不说二重积分能不能像你那么算，你的f就不对了，而且后面的表达式已经看不出来变量了。看看原式能不能简化一下，不然你自己打也容易出错。

追问

私信你了