boxoode.m

蒙特卡罗方法及其例子

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];