def производная(X, t, A, B, C, D): х, у = Х dotx = x * (A - B * y) doty = y * (-D + C * x) вернуть np.array( ) интеграция по определению (t,A,B,C,D,X0): res = odeint(производная, X0, t, args = (A,B,C,D)) вернуть разрешение Х0 = X = массив([ , , , ,...