import math
K1=0.09
K2=0.005
K3=2
W=50
P=17.8
D=0.58
M=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('характеристика каскада ферментеров')
import numpy as np
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)
print(S0)