import tellurium as _te

r = _te.loada("""
   J1: -> x1; vm/(km + x3^n);
    x1 -> x2; k1*x1
    x2 -> x3; k2*x2
    x3 ->;    k3*x3;

    k1 = 0.5; k2 = 0.5; k3 = 0.5;
    n = 10;
    x1 = 2; x2 = 2; x3 = 2
    vm = 18; km = 0.5; 
""")

m = r.simulate (0, 30, 1000)
r.plot()
