x y的關系曲線如附圖1所示。

讓計算機做無窮加和，是不可能實現的。
實際上n的取值只要足夠大，原加和式就收斂了。
一下代碼按n=9999編寫的：
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(,poly),axis( );
figure(2)
plot(,p),axis( );
figure(3)
plot(x_data,y,'r-'),xlabel('x'),ylabel('y');