带有二个变参数的二重积分问题，以两变参数数值积分后绘图！

ljh2011

二元函数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)

不知错在哪，请各位指点！！

qibbxxt

1. ezsurf(@(k,m)quad2d(@(x,y)sin(k.*x).*cos(m.*y).*x.^2,1,2,1,5))

复制代码