In [ ]:

    

from pyodeint import integrate_predefined
import numpy as np
from van_der_pol import get_f_and_j


    

In [ ]:

    

help(integrate_predefined)


    

In [ ]:

    

rhs, jac = get_f_and_j(1.0)
xout = np.linspace(0, 10, 200)
yout, info = integrate_predefined(rhs, jac, [0, 1], xout, dx0=1e-12, atol=1e-8, rtol=1e-8)
print(info['success'])


    

In [ ]:

    

import matplotlib.pyplot as plt
%matplotlib inline
_ = plt.plot(xout, yout)