Matlab diary for pi^2/6 lecture
% How to estimate pi^2/6 from N terms of series 1 + 1/2^2 + 1/3^2 + 1/4^2 + ...
% Using a for-loop
sum=0;
N=100;
for n=1:N
sum = sum + 1/n^2;
end
sum
sum = 1.6350
pi^2/6
ans = 1.6449
sum-pi^2/6
ans = -0.0100 % off by about 0.01 for N=100 terms --not so great
% To try different values of number of terms, write a function
% see file pisquared6.m
pisquared6(100)
ans =
1.6350
pisquared6(100) -pi^2/6
ans =
-0.0100
pisquared6(1000) -pi^2/6
ans =
-9.9950e-04
pisquared6(10000) -pi^2/6
ans =
-9.9995e-05
% Even better, produce a plot of value of series versus N
% See pigraph.m
% Now, do the same calculations with one line of Matlab
sum((1:10).^(-2)) % sum 1/n^2 over n=1 through 10
ans =
1.5498
sum((1:100).^(-2)) % sum 1/n^2 over n=1 through 100
ans =
1.6350
sum((1:100).^(-2)) - pi^2/6
ans =
-0.0100
sum((1:1000).^(-2)) - pi^2/6
ans =
-9.9950e-04
sum((1:10000).^(-2)) - pi^2/6
ans =
-9.9995e-05
sum((1:100000).^(-2)) - pi^2/6
ans =
-9.9999e-06