решение двумерных дифференциальных уравнений в питоне с помощью scipy
я новичок в питоне. У меня есть простая дифференциальная система, которая состоит из двух переменных и двух дифференциальных уравнений и начальных условийx0=1, y0=2
:
dx/dt=6*y
dy/dt=(2t-3x)/4y
Теперь я пытаюсь решить эти два дифференциальных уравнения, и я выбираюodeint
, Вот мой код:
import matplotlib.pyplot as pl
import numpy as np
from scipy.integrate import odeint
def func(z,b):
x, y=z
return [6*y, (b-3*x)/(4*y)]
z0=[1,2]
t = np.linspace(0,10,11)
b=2*t
xx=odeint(func, z0, b)
pl.figure(1)
pl.plot(t, xx[:,0])
pl.legend()
pl.show()
но результат неверный и появляется сообщение об ошибке:
Excess work done on this call (perhaps wrong Dfun type).
Run with full_output = 1 to get quantitative information.
Я не знаю, что не так с моим кодом и как я могу его решить. Любая помощь будет полезна для меня.