function dy = boxoODE(t,y,theta) % ode function for the boxo example A = y(1); B = y(2); k1mean = theta(1); e1 = theta(2); k2mean = theta(3); e2 = theta(4); tmean = theta(5); temp = theta(6); e1 = e1*1e6; e2 = e2*1e6; R = 8.314; z = 1./R *(1./temp-1./tmean); k1 = k1mean * exp(-e1*z); k2 = k2mean * exp(-e2*z); dy =[ - k1 * A ; k1 * A - k2 * B];