from math import factorial
x=[1.50,1.51,1.52,1.53,1.54,1.55,1.56,1.57,1.58,1.59,1.60]
y=[0.51183,0.50624,0.50064,0.49503,0.48940,0.48376,0.47811,0.47245,0.46678,0.46100,0.45540]
a=(float(input()))
b=round(a,2)
print(b)
q=abs((a-b)/0.01)
print(q)
d1=[]
d2=[]
d3=[]
d4=[]
d5=[]
d6=[]
d7=[]
d8=[]
d9=[]
d10=[]
d11=[]

vy=[]

for i in range (0,len(y)-1):
    d1.append(abs(y[i+1]-y[i]))
print(d1)

for i in range (0,len(d1)-1):
    d2.append(abs(d1[i+1]-d1[i]))
print(d2)

for i in range (0,len(d2)-1):
    d3.append(abs(d2[i+1]-d2[i]))
print(d3)

for i in range (0,len(d3)-1):
    d4.append(abs(d3[i+1]-d3[i]))
print(d4)

for i in range (0,len(d4)-1):
    d5.append(abs(d4[i+1]-d4[i]))
print(d5)

for i in range (0,len(d5)-1):
    d6.append(abs(d5[i+1]-d5[i]))
print(d6)

for i in range (0,len(d6)-1):
    d7.append(abs(d6[i+1]-d6[i]))
print(d7)

for i in range (0,len(d7)-1):
    d8.append(abs(d7[i+1]-d7[i]))
print(d8)

for i in range (0,len(d8)-1):
    d9.append(abs(d8[i+1]-d8[i]))
print(d9)

for i in range (0,len(d9)-1):
    d10.append(abs(d9[i+1]-d9[i]))
print(d10)

Pnx=y[0] +d1[0]*q + d2[0]*(q-1)/factorial(2) + d3[0]*(q-2)*(q-1)/factorial(3) + d4[0]*(q-3)*(q-2)*(q-1)/factorial(4) + d5[0]*(q-4)*(q-3)*(q-2)*(q-1)/factorial(5) + d6[0]*(q-5)*(q-4)*(q-3)*(q-2)*(q-1)/factorial(6) + d7[0]*(q-6)*(q-5)*(q-4)*(q-3)*(q-2)*(q-1)/factorial(7) + d8[0]*(q-7)*(q-6)*(q-5)*(q-4)*(q-3)*(q-2)*(q-1)/factorial(8) + d9[0]*(q-8)*(q-7)*(q-6)*(q-5)*(q-4)*(q-3)*(q-2)*(q-1)/factorial(9) + d10[0]*(q-9)*(q-8)*(q-7)*(q-6)*(q-5)*(q-4)*(q-3)*(q-2)*(q-1)/factorial(10)
print(Pnx)