from math import sqrt def archimede(n): nb_cotes=6 u=0.5 p=3 for i in range(n): u=sqrt(0.5*(1-sqrt(1-u**2))) nb_cotes=6*2**(i+1) p=nb_cotes*u return nb_cotes,p