[fanxing39]

主程序：clear
clc
[x,y]=meshgrid(-10:0.5:10);
[m,n]=size(x);
for i=1:m
for j=1:n
P(i,j)=quad(@(s)fun2(s,x(i,j),y(i,j)),0,2);
end
end
surf(P)
副程序：function f=fun2(s,x,y)
f=sin(s).^2.*exp((x.*sin(s).^2+y.*cos(s).^2)./sqrt(1+sin(s)));