ModelingToolkit Simple Chemical Reaction ODE (B)

using ModelingToolkit
using ModelingToolkit: t_nounits as t, D_nounits as D
using DifferentialEquations, Plots

@parameters k1=1.71 k2=280.0 k3=8.32 k4=0.69 k5=0.43 k6=1.81
@variables y1(t)=1.0 y2(t)=0.0 y3(t)=0.0 y4(t)=0.0 y5(t)=0.0 y6(t)=0.0 y7(t)=0.0 y8(t)=0.0057

eqs = [
    D(y1) ~ (-k1*y1 + k5*y2 + k3*y3 + 0.0007)
    D(y2) ~ (k1*y1 - 8.75*y2)
    D(y3) ~ (-10.03*y3 + k5*y4 + 0.035*y5)
    D(y4) ~ (k3*y2 + k1*y3 - 1.12*y4)
    D(y5) ~ (-1.745*y5 + k5*y6 + k5*y7)
    D(y6) ~ (-k2*y6*y8 + k4*y4 + k1*y5 - k5*y6 + k4*y7)
    D(y7) ~ (k2*y6*y8 - k6*y7)
    D(y8) ~ (-k2*y6*y8 + k6*y7)
]

@named hires_sys = ODESystem(eqs, t)
hires_model = mtkcompile(hires_sys)

prob = ODEProblem(hires_model, [], (0.0, 321.8122))
sol = solve(prob)
pp = plot(sol, dpi=400)