import math
K1=0.09
K2=0.005
K3=2
W=50
P=17.8
D0=0.58
M0=0
N=int(input("N="))
print('Шаг интегрирования',N)
V0=0
S0=0
print('Каскада ферментеров P9')
P9=float(input("P9="))
print('Конечное значение',P9)
X0=P9/P
X=X0**(1/N)
print('P9=',P, ",", 'X0=', X0, ",", 'X=', X)
print('Характеристика каскада ферментеров')

for i in range(N):
    P = X*P
    M = M0+(K2/K1)*((1/X)-1)
    D = D0+(K1*P-K2)*(M0-M)/(-K2)
    T = (M0-M)/(-K2*D)
    V = W*T*K3
    H=(4*V/3.14)**(0.33)
    S=1.5*3.14*(H/2)
    V0=V0+V
    S0=S0+S
    D0=D
    M0=M
    P0=P
    print('P=', P, ",", 'D=', D, ",",  'M=', M, ",", 'H=', H, ",", 'S=', S, ",", 'V=', V, "." )
    print()

print("V0=", V0)
print("S0=", S0)