import numpy as np def main(): n=input("Numero di iterazioni:\n") n=int(n) l0_buono=np.sqrt(np.float128(2.0)) for i in range(n): l1_buono=np.sqrt((l0_buono/2.0)**2+(1.0-np.sqrt(1.0-(l0_buono/2.0)**2))**2) l0_buono=l1_buono pi_buono=l1_buono*4*2.0**(i) print(str(i)+ " "+str(pi_buono)) if __name__ == "__main__": main()