MATLAB: Displaying output from a summation.

summation

I'm trying to find partial sums of s=1/k^2 where k= multiples of 400. I'm using a "for" statement and want to have the script output the values at those multiples. At this point, I can only get it to output the final value. I want it to output the sums for k=400, 800, 1200, … Also, I would like the output for the error for the actual sum.
Code is below.
format long
s=0;
for n=1:13
    k=400;
end
for i=1:n*k 
    s=s+1/i^2;
    y=abs(s-(pi^2/6));
end
    disp(['         Error', '               s_n'])
    disp([y,s])
Thanks.

Best Answer

fprintf('%20s%20s\n', 'Error', 's_n');
s = 0;
n = 13;    % The FOR loop is meaningless?!
k = 400;
r = pi^2/6;
for i = 1:n*k 
    s = s+1/i^2;
    if mod(i, 400) == 0
       y = abs(s-r);
       fprintf('%20g%20g\n', y, s);
    end
end
Related Question