n = 10;
pn = 0;
format long
k=0:n;
termivekt=(1./16.^k).*(4./(8*k+1)-2./(8*k+4)-1./(8*k+5)-1./(8*k+6))
sum(termivekt)
[k' cumsum(termivekt)']

format  %Palautetaan oletusformaatti
