n=0:40;
s=(-1).^n./sqrt(n+1);
ss=cumsum(s);
si=ss;
for i=1:30
	si=(si(1:end-1)+si(2:end))/2;			
end;
si(end)
