Integrate a system of ordinary differential equations.
Solve a system of ordinary differential equations using lsoda from the FORTRAN library odepack.
Solves the initial value problem for stiff or nonstiff systems of first order odes:
dy/dt = func(y,t0,...)
where y can be a vector.
Parameters:  func : callable(y, t0, ...)
y0 : array
t : array
args : tuple
Dfun : callable(y, t0, ...)
col_deriv : boolean
full_output : boolean
printmessg : boolean



Returns:  y : array, shape (len(y0), len(t))
infodict : dict, only returned if full_output == True
