Eu sou um novato em python. Eu tenho um sistema diferencial simples, que consiste em duas variáveis e duas equações diferenciais e condições iniciaisx0=1, y0=2:

dx/dt=6*y
dy/dt=(2t-3x)/4y

agora estou tentando resolver essas duas equações diferenciais e escolhoodeint. Aqui está o meu código:

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()

mas o resultado está incorreto e há uma mensagem de erro:

Excess work done on this call (perhaps wrong Dfun type).
Run with full_output = 1 to get quantitative information.

Não sei o que há de errado com meu código e como posso resolvê-lo. Qualquer ajuda será útil para mim.