import numpy as np
import matplotlib.pyplot as plt
import monedo as mo


def f(t,x):
    k1, k2, k3 =0.04, 3*10**7, 10**4
    return np.array([-k1*x[0]+k3*x[1]*x[2], k1*x[0]-k3*x[1]*x[2]-k2*x[1]**2,
        k2*x[1]**2])

y0=np.array([1. ,0., 0.])

Y=mo.euler_all(f,0,.3,y0,300)
plt.plot(np.linspace(0,.3,301),Y[1,:])
plt.show()