参考代码:
syms x y uf=x^3*(100-2*u)/(pi*(x^2+(y-u)^2)^2);X=-10:10;Y=-10:10;Z=zeros(length(X),length(Y));for i=1:length(X) for j=1:length(Y) z=int(subs(f,{x,y},{X(i),Y(j)}),u,-50,50); Z(i,j)=z; endendmesh(X,Y,Z)
注意自行选择变量x和y的范围。