带有二个变参数的二重积分问题,以两变参数数值积分后绘图!
二元函数f(x,y) =sin(k.*x).*cos(m.*y).*x.^2,y从1积分到5,x从1积分到2;代码:
kk=linspace(0,5);
mm=linspace(0,5);
y=zeros(size(kk));
ff=@(k,m) ['sin(',num2str(k),'*x).*cos(',num2str(m),'*y).*x.^2'];
f=@(k,m) quadl(@(x) arrayfun(@(xx) quadl(ff(k,m),1,5),x),1,2);
for ii=1:length(kk)
for jj=1:length(mm)
y(ii,jj)=f(kk(ii),mm(jj));
end
end
plot3(kk,mm,y)
不知错在哪,请各位指点!! ezsurf(@(k,m)quad2d(@(x,y)sin(k.*x).*cos(m.*y).*x.^2,1,2,1,5))
页:
[1]