import numpy as np def main(): n=input("Numero di iterazioni:\n") n=int(n) l0_buono=np.sqrt(2.0) l0_ciofeca=np.sqrt(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) l1_ciofeca=np.sqrt(2-2*np.sqrt(1-(l0_ciofeca/2)**2)) l0_ciofeca=l1_ciofeca pi_ciofeca=l1_ciofeca*4*2.0**(i) print(str(i)+ " "+str(pi_buono)+" "+str(pi_ciofeca)) if __name__ == "__main__": main()