u=0;p=0;
while u<1.644452
  p=p+1;
  u=u+1/p^2;
end
disp(p-1);
