s1fig2.m

Industrial Mathematics

%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)