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');