def calculs(a,b,N): u,v,n = a,b,0 while n < N: n = n + 1 u = (a+b)/2 v = ( (a*a + b*b)/2 )**0.5 a,b = u,v return u,v