function plot_n
clear all;clc
format long
nx=1000;   % xÈ¡µãµÄ¸öÊý
x_data=linspace(0,250,nx);
m=length(x_data);

for k=1:m
    x=x_data(k);
n=9999;
poly(1)=8/pi/pi*(2*0+1)^(-2)*exp(-9.13e-13*x/26*pi^2*(2*0+1)^2   );
for i=1:1:n
poly(i+1)=8/pi/pi*(2*i+1)^(-2)*exp(-9.13e-13*x/26*pi^2*(2*i+1)^2   );
end

for j=1:n+1
    p(j)=sum( poly(1:1:j) );
end

y(k)=16.5-16.5*p(end);
end

figure(1)
plot([1:n+1],poly),axis([-10 n+1 -0.1 0.9] );

figure(2)
plot([1:n+1],p),axis([-10 n+1 0.8 1.05] );

figure(3)
plot(x_data,y,'r-'),xlabel('x'),ylabel('y');

3Â¥: Originally posted by ÔÂÖ»À¶ at 2016-01-18 14:25:01
x yµÄ¹ØϵÇúÏßÈ總ͼ1Ëùʾ¡£

¸½Í¼1.png

4Â¥: Originally posted by ³ÂçùÍõ´ó¿É at 2016-01-18 14:27:39
ÎÒ°Ñ´úÂëÒ²ÔËÐÐÁËһϣ¬ÕâÊÇfigure3µÄͼ°É£¬£¬£¬1ºÍ2ÊÇɶÀÕ...