p=int(input("p =")) pas=10**-p h=1 while h>0: c=1/2*h**2+1 print("h =",h,"c =",c) h=round(h-pas,p)