s1fig2.m

来自「Industrial Mathematics」· M 代码 · 共 12 行

%Figure 1.2
x=1:.01:12;
y=0*x;
f=1:12;
f=0*f
for i=1:12
 f(i+1) = f(i)+ 4/(i*i*pi*pi);
  end;
 y = f(round(x));
plot(x,y)