各位高手,我现在编了一程序怎么作出的图是一条直线啊,请指教
function yout=plot112
a=zeros(20,1);
n=1;
for i=1:20
a(i)=1/(pi*500000000000*n.^6).*dblquad(@fun,0,15,0,2*pi,[],[],n);
i
n=n+1;
end
n=1:20;
plot(n,a);
function z= fun(x,y,n)
z=x.*(x.*cos(y)+n.^2*2000).^2.*(x.*sin(y)+n.^2*2000).*exp(-1/(n.^2*2000)-2*x.*cos(y)-2*x.*sin(y)-40000*n.^2);
谢谢大家