 %  DOQCS : http://doqcs.ncbs.res.in/ 
 %  Accession Name = CaMKII_noPKA_model3 
 %  Accession Number = 62 
 %  Transcriber = Upinder S. Bhalla, NCBS 
 %  Developer = Arnold Hayer and Upinder S. Bhalla, NCBS 
 %  Species = General Mammalian 
 %  Tissue = Neuronal 
 %  Cell Compartment = Synaptic Spine, Postsynaptic Density 
 %  Notes =  This is the model of CaMKII bistability, model 3. It exhibits   bistability in CaMKII activation due to autophosphorylation at the   PSD and local saturation of PP1. This version   of model 3 does not include the full PKA regulatory pathway,   and instead has a predefined initial amount of active PKA.  
 
   % Note: All the concentrations are in micro Mole (uM)

function moose_convert
%Matlab simulation for chemical kinetics

% Saved on Wed Aug  8 16:39:59 2007
% Automatically generated from kkit model in MOOSE by conv_matlab

global buff;
global t;
global y;

% Define Tables here
tspan=[0, 2000];
% Molecule names
molnames = {
'1 /kinetics/CaM-Ca4'
'2 /kinetics/PP1-active'
'3 /kinetics/PP1-active/Deph-thr286/Deph-thr286_cplx'
'4 /kinetics/PP1-active/Deph-thr286b/Deph-thr286b_cplx'
'5 /kinetics/PP1-active/Deph-thr305/Deph-thr305_cplx'
'6 /kinetics/PP1-active/Deph-thr286c/Deph-thr286c_cplx'
'7 /kinetics/PP1-active/Deph-thr305a/Deph-thr305a_cplx'
'8 /kinetics/Ca'
'9 /kinetics/PKA-active'
'10 /kinetics/PKA-active/PKA-phosph-I1/PKA-phosph-I1_cplx'
'11 /kinetics/PKA-active/PKA-phosph-I1_PSD/PKA-phosph-I1_PSD_cplx'
'12 /kinetics/CaMKII/CaMKII'
'13 /kinetics/CaMKII/CaMKII-CaM'
'14 /kinetics/CaMKII/CaMKII-thr286*-CaM'
'15 /kinetics/CaMKII/CaMKII***'
'16 /kinetics/CaMKII/CaMKII-thr286'
'17 /kinetics/CaMKII/tot_CaM_CaMKII/CaM_act_305/CaM_act_305_cplx'
'18 /kinetics/CaMKII/tot_CaM_CaMKII/CaM_act_286/CaM_act_286_cplx'
'19 /kinetics/CaMKII/tot_autonomous_CaMKII/auton_305/auton_305_cplx'
'20 /kinetics/CaMKII/tot_autonomous_CaMKII/auton_286/auton_286_cplx'
'21 /kinetics/CaMKII/CaMK-thr305'
'22 /kinetics/CaM/CaM'
'23 /kinetics/CaM/CaM-PSD'
'24 /kinetics/CaM-Ca3'
'25 /kinetics/CaM-TR2-Ca2'
'26 /kinetics/PP1/I1'
'27 /kinetics/PP1/I1*'
'28 /kinetics/PP1/PP1-I1*'
'29 /kinetics/PP1/PP1-I1'
'30 /kinetics/PP2A'
'31 /kinetics/PP2A/PP2A-dephosph-I1/PP2A-dephosph-I1_cplx'
'32 /kinetics/PP2A/PP2A-dephosph-PP1-I*/PP2A-dephosph-PP1-I*_cplx'
'33 /kinetics/PP2A/PP2A-dephosph-I1_PSD/PP2A-dephosph-I1_PSD_cplx'
'34 /kinetics/PP2A/PP2A-dephosph-PP1-I*_PSD/PP2A-dephosph-PP1-I*_PSD_cplx'
'35 /kinetics/CaNAB-Ca4'
'36 /kinetics/CaNAB-Ca4/dephosph_inhib1_noCaM/dephosph_inhib1_noCaM_cplx'
'37 /kinetics/CaNAB-Ca4/dephosph_inhib1_noCaM_PSD/dephosph_inhib1_noCaM_PSD_cplx'
'38 /kinetics/PP2B/CaNAB'
'39 /kinetics/PP2B/CaNAB-Ca2'
'40 /kinetics/CaMKII-thr286-CaM-PSD'
'41 /kinetics/CaMKII-CaM-PSD'
'42 /kinetics/CaMKII-thr286-PSD'
'43 /kinetics/CaMKII-PSD'
'44 /kinetics/CaMKII***-PSD'
'45 /kinetics/tot-auto-PSD/auton-305-PSD/auton-305-PSD_cplx'
'46 /kinetics/tot-auto-PSD/auton-286-PSD/auton-286-PSD_cplx'
'47 /kinetics/tot-auto-PSD/auton-autoph/auton-autoph_cplx'
'48 /kinetics/CaM-TR2-Ca2-PSD'
'49 /kinetics/CaM-Ca3-PSD'
'50 /kinetics/CaM-Ca4-PSD'
'51 /kinetics/Ca-PSD'
'52 /kinetics/PP1_PSD/I1'
'53 /kinetics/PP1_PSD/I1*'
'54 /kinetics/PP1_PSD/PP1-I1*'
'55 /kinetics/PP1_PSD/PP1-I1'
'56 /kinetics/PP1-active_PSD'
'57 /kinetics/PP1-active_PSD/Deph-thr286/Deph-thr286_cplx'
'58 /kinetics/PP1-active_PSD/Deph-thr286b/Deph-thr286b_cplx'
'59 /kinetics/PP1-active_PSD/Deph-thr286c/Deph-thr286c_cplx'
'60 /kinetics/PP1-active_PSD/Deph-thr305a/Deph-thr305a_cplx'
'61 /kinetics/PP1-active_PSD/Deph-thr305/Deph-thr305_cplx'
'62 /kinetics/NMDAR'
'63 /kinetics/CaM_Ca_n-CaNAB'
'64 /kinetics/CaM_Ca_n-CaNAB/dephosph_inhib1/dephosph_inhib1_cplx'
'65 /kinetics/CaM_Ca_n-CaNAB/dephosph-PP1-I*/dephosph-PP1-I*_cplx'
'66 /kinetics/CaM_Ca_n-CaNAB/dephosph_inhib1_PSD/dephosph_inhib1_PSD_cplx'
'67 /kinetics/CaM_Ca_n-CaNAB/dephosph-PP1-I*_PSD/dephosph-PP1-I*_PSD_cplx'
'68 /kinetics/basal_CaMKII_PSD'
'69 /kinetics/tot-CaM-CaMKII-PSD/CaM-act-305-PSD/CaM-act-305-PSD_cplx'
'70 /kinetics/tot-CaM-CaMKII-PSD/CaM-act-286-PSD/CaM-act-286-PSD_cplx'
'71 /kinetics/tot-CaM-CaMKII-PSD/CaM_act_autoph/CaM_act_autoph_cplx'
'72 /kinetics/CaMKII-thr305-PSD'
};

bufnames = {
'1 /kinetics/basal_CaMKII_cyt'
'2 /kinetics/Ca_control_cyt'
'3 /kinetics/Ca_control_PSD'
'4 /kinetics/basal_CaMKII_PSD_control'
};

sumtotnames = {
'1 /kinetics/CaMKII/tot_CaM_CaMKII'
'2 /kinetics/CaMKII/tot_autonomous_CaMKII'
'3 /kinetics/tot-auto-PSD'
'4 /kinetics/286P-PSD'
'5 /kinetics/actCaMKII-PSD'
'6 /kinetics/tot_CaMKII_PSD'
'7 /kinetics/tot_CaMKII_cyt'
'8 /kinetics/act_CaMKII_cyt'
'9 /kinetics/tot-CaM-CaMKII-PSD'
};

% Initial concs of 72 molecules
y0(1) = 0;
y0(2) = 1.8;
y0(3) = 0;
y0(4) = 0;
y0(5) = 0;
y0(6) = 0;
y0(7) = 0;
y0(8) = 0.08;
y0(9) = 0.0185185;
y0(10) = 0;
y0(11) = 0;
y0(12) = 20;
y0(13) = 0;
y0(14) = 0;
y0(15) = 0;
y0(16) = 0;
y0(17) = 6.3263e-11;
y0(18) = 2.44537e-08;
y0(19) = 4.90056e-08;
y0(20) = 1.89444e-05;
y0(21) = 0;
y0(22) = 26.3333;
y0(23) = 26.3333;
y0(24) = 0;
y0(25) = 0;
y0(26) = 1.8;
y0(27) = 0;
y0(28) = 0;
y0(29) = 0;
y0(30) = 0.111111;
y0(31) = 1.99333e-06;
y0(32) = 0.00691;
y0(33) = 2.21556e-07;
y0(34) = 0.00127907;
y0(35) = 0;
y0(36) = 0;
y0(37) = 0;
y0(38) = 1;
y0(39) = 0;
y0(40) = 0;
y0(41) = 0;
y0(42) = 0;
y0(43) = 0;
y0(44) = 0;
y0(45) = 0;
y0(46) = 0;
y0(47) = 0;
y0(48) = 0;
y0(49) = 0.00254583;
y0(50) = 0;
y0(51) = 0.08;
y0(52) = 4;
y0(53) = 0;
y0(54) = 0;
y0(55) = 0;
y0(56) = 4;
y0(57) = 0;
y0(58) = 0;
y0(59) = 0;
y0(60) = 0;
y0(61) = 0;
y0(62) = 120;
y0(63) = 0;
y0(64) = 0;
y0(65) = 0;
y0(66) = 0;
y0(67) = 0;
y0(68) = 2;
y0(69) = 0;
y0(70) = 0;
y0(71) = 0;
y0(72) = 0;
buff(1) = 2;
buff(2) = 0.08;
buff(3) = 0.08;
buff(4) = 2;
% Sumtotal eqn for /kinetics/CaMKII/tot_CaM_CaMKII: sumtot(1) =  (y(13) + y(14) )
% Sumtotal eqn for /kinetics/CaMKII/tot_autonomous_CaMKII: sumtot(2) =  (y(16) + y(15) + 1 * buff(1) )
% Sumtotal eqn for /kinetics/tot-auto-PSD: sumtot(3) =  (y(42) + y(44) + y(68) )
% Sumtotal eqn for /kinetics/286P-PSD: sumtot(4) =  (y(40) + y(42) )
% Sumtotal eqn for /kinetics/actCaMKII-PSD: sumtot(5) =  (y(40) + y(41) + 1 *  (y(42) + y(44) + y(68) ) )
% Sumtotal eqn for /kinetics/tot_CaMKII_PSD: sumtot(6) =  (1 *  (y(40) + y(41) + 1 *  (y(42) + y(44) + y(68) ) ) + y(43) + y(72) )
% Sumtotal eqn for /kinetics/tot_CaMKII_cyt: sumtot(7) =  (y(13) + y(14) + y(16) + y(15) + y(21) + y(12) + 1 * buff(1) )
% Sumtotal eqn for /kinetics/act_CaMKII_cyt: sumtot(8) =  (1 *  (y(13) + y(14) ) + 1 *  (y(16) + y(15) + 1 * buff(1) ) )
% Sumtotal eqn for /kinetics/tot-CaM-CaMKII-PSD: sumtot(9) =  (y(40) + y(41) )

% Call the ODE solver and display results
[t, y] = ode23s(@f, tspan, y0);

hold on;
do_plots;

function do_plots()
global t;
global y;
global buff;
plot(t, (y(:,40) + y(:,41) + 1 *  (y(:,42) + y(:,44) + y(:,68) ) ));
plot(t, (1 *  (y(:,40) + y(:,41) + 1 *  (y(:,42) + y(:,44) + y(:,68) ) ) + y(:,43) + y(:,72) ));
plot(t, (1 *  (y(:,13) + y(:,14) ) + 1 *  (y(:,16) + y(:,15) + 1 * buff(1) ) ));
plot(t, (y(:,13) + y(:,14) + y(:,16) + y(:,15) + y(:,21) + y(:,12) + 1 * buff(:,1) ));
plot(t,y(:,63));
plot(t,y(:,50));
plot(t,y(:,56));
plot(t,y(:,9));
plot(t,y(:,8));
plot(t,y(:,2));
plot(t,y(:,1));

% Evaluation function
function dydt = f(t, y)
global buff;
dydt = [
 - 49.9997 * y(1) * y(12) + 5 * y(13) - 1000.19 * y(16) * y(1) + 0.1 * y(14) + 0.464999 * y(24) * y(8) - 10 * y(1) - 599.994 * y(1) * y(35) + 1 * y(63) + 60 * y(50) - 60 * y(1)
 - 499.981 * y(27) * y(2) + 0.1 * y(28) + 1 * y(29) - 2.45144 * y(14) * y(2) + 10 * y(3) + 2.5 * y(3) - 2.45144 * y(16) * y(2) + 10 * y(4) + 2.5 * y(4) - 2.45144 * y(15) * y(2) + 10 * y(5) + 2.5 * y(5) - 2.45144 * y(15) * y(2) + 10 * y(6) + 2.5 * y(6) - 2.45144 * y(21) * y(2) + 10 * y(7) + 2.5 * y(7)
 + 2.45144 * y(14) * y(2) - 10 * y(3) - 2.5 * y(3)
 + 2.45144 * y(16) * y(2) - 10 * y(4) - 2.5 * y(4)
 + 2.45144 * y(15) * y(2) - 10 * y(5) - 2.5 * y(5)
 + 2.45144 * y(15) * y(2) - 10 * y(6) - 2.5 * y(6)
 + 2.45144 * y(21) * y(2) - 10 * y(7) - 2.5 * y(7)
 - 143.998 * y(22) * y(8)*y(8) + 144 * y(25) - 3.60002 * y(25) * y(8) + 10 * y(24) - 0.464999 * y(24) * y(8) + 10 * y(1) - 7.20019 * y(8)*y(8) * y(39) + 2 * y(35) - 20016 * y(38) * y(8)*y(8) + 2 * y(39) + 100 * buff(2) - 100 * y(8)
 - 5.99994 * y(26) * y(9) + 36 * y(10) + 9 * y(10) - 0.66666 * y(52) * y(9) + 36 * y(11) + 9 * y(11)
 + 5.99994 * y(26) * y(9) - 36 * y(10) - 9 * y(10)
 + 0.66666 * y(52) * y(9) - 36 * y(11) - 9 * y(11)
 - 49.9997 * y(1) * y(12) + 5 * y(13) + 0.0333333 * y(43) - 6e-05 * y(12) * y(62) + 2.5 * y(4) + 2.5 * y(7)
 + 49.9997 * y(1) * y(12) - 5 * y(13) - 0.0110133 * y(13) *  (y(13) + y(14) ) + 2 * y(18) - 0.00714312 * y(13) *  (y(16) + y(15) + 1 * buff(1) ) + 2 * y(20) - 0.00012 * y(13) * y(62) + 2.5 * y(3)
 + 1000.19 * y(16) * y(1) - 0.1 * y(14) + 0.5 * y(18) + 0.5 * y(20) - 0.00012 * y(14) * y(62) - 2.45144 * y(14) * y(2) + 10 * y(3)
 + 6 * y(17) + 6 * y(19) - 2.45144 * y(15) * y(2) + 10 * y(5) - 2.45144 * y(15) * y(2) + 10 * y(6)
 - 1000.19 * y(16) * y(1) + 0.1 * y(14) - 0.13216 * y(16) *  (y(13) + y(14) ) + 24 * y(17) - 0.0857142 * y(16) *  (y(16) + y(15) + 1 * buff(1) ) + 24 * y(19) - 2.45144 * y(16) * y(2) + 10 * y(4) + 2.5 * y(5)
 + 0.13216 * y(16) *  (y(13) + y(14) ) - 24 * y(17) - 6 * y(17)
 + 0.0110133 * y(13) *  (y(13) + y(14) ) - 2 * y(18) - 0.5 * y(18)
 + 0.0857142 * y(16) *  (y(16) + y(15) + 1 * buff(1) ) - 24 * y(19) - 6 * y(19)
 + 0.00714312 * y(13) *  (y(16) + y(15) + 1 * buff(1) ) - 2 * y(20) - 0.5 * y(20)
 + 0.0333333 * y(72) - 6e-05 * y(21) * y(62) - 2.45144 * y(21) * y(2) + 10 * y(7) + 2.5 * y(6)
 - 71.999 * y(22) * y(8)*y(8) + 72 * y(25)
 - 72 * y(23) * y(51)*y(51) + 72 * y(48)
 + 3.60002 * y(25) * y(8) - 10 * y(24) - 0.464999 * y(24) * y(8) + 10 * y(1)
 + 71.999 * y(22) * y(8)*y(8) - 72 * y(25) - 3.60002 * y(25) * y(8) + 10 * y(24)
 + 0.034 * y(36) + 0.34 * y(64) + 2 * y(31) + 1 * y(29) - 5.99994 * y(26) * y(9) + 36 * y(10)
 - 499.981 * y(27) * y(2) + 0.1 * y(28) - 0.0341998 * y(27) * y(35) + 0.136 * y(36) - 0.341998 * y(27) * y(63) + 1.36 * y(64) - 0.64584 * y(27) * y(30) + 8.3334 * y(31) + 9 * y(10)
 + 499.981 * y(27) * y(2) - 0.1 * y(28) - 0.64584 * y(28) * y(30) + 8.3334 * y(32) - 0.341998 * y(28) * y(63) + 1.36 * y(65)
 - 1 * y(29) + 2 * y(32) + 0.34 * y(65)
 - 0.64584 * y(27) * y(30) + 8.3334 * y(31) + 2 * y(31) - 0.64584 * y(28) * y(30) + 8.3334 * y(32) + 2 * y(32) - 0.07176 * y(53) * y(30) + 8.3334 * y(33) + 2 * y(33) - 0.07176 * y(54) * y(30) + 8.3334 * y(34) + 2 * y(34)
 + 0.64584 * y(27) * y(30) - 8.3334 * y(31) - 2 * y(31)
 + 0.64584 * y(28) * y(30) - 8.3334 * y(32) - 2 * y(32)
 + 0.07176 * y(53) * y(30) - 8.3334 * y(33) - 2 * y(33)
 + 0.07176 * y(54) * y(30) - 8.3334 * y(34) - 2 * y(34)
 - 0.0341998 * y(27) * y(35) + 0.136 * y(36) + 0.034 * y(36) + 3.60009 * y(8)*y(8) * y(39) - 1 * y(35) - 599.994 * y(1) * y(35) + 1 * y(63) - 0.00380004 * y(53) * y(35) + 0.136 * y(37) + 0.034 * y(37)
 + 0.0341998 * y(27) * y(35) - 0.136 * y(36) - 0.034 * y(36)
 + 0.00380004 * y(53) * y(35) - 0.136 * y(37) - 0.034 * y(37)
 - 10008 * y(38) * y(8)*y(8) + 1 * y(39)
 + 10008 * y(38) * y(8)*y(8) - 1 * y(39) - 3.60009 * y(8)*y(8) * y(39) + 1 * y(35)
 + 0.00108 * y(14) * y(62) + 1000.02 * y(42) * y(50) - 0.1 * y(40) + 0.5 * y(70) + 0.5 * y(46) - 0.499998 * y(40) * y(56) + 0.8 * y(57)
 + 0.00108 * y(13) * y(62) - 0.0078126 * y(41) *  (y(40) + y(41) ) + 2 * y(70) - 0.00499998 * y(41) *  (y(42) + y(44) + y(68) ) + 2 * y(46) + 0.2 * y(57) - 5 * y(41) + 49.9998 * y(43) * y(50)
 - 0.09375 * y(42) *  (y(40) + y(41) ) + 24 * y(69) - 0.06 * y(42) *  (y(42) + y(44) + y(68) ) + 24 * y(45) - 1000.02 * y(42) * y(50) + 0.1 * y(40) - 0.499998 * y(42) * y(56) + 0.8 * y(58) + 0.2 * y(61) + 2 * y(71) + 2 * y(47)
 - 0.3 * y(43) + 0.00054 * y(12) * y(62) - 49.9998 * y(43) * y(50) + 0.2 * y(60) + 0.2 * y(58) + 5 * y(41) - 0.0312498 * y(43) *  (y(40) + y(41) ) + 8 * y(71) - 0.0199998 * y(43) *  (y(42) + y(44) + y(68) ) + 8 * y(47)
 + 6 * y(69) + 6 * y(45) - 0.499998 * y(44) * y(56) + 0.8 * y(61) - 0.499998 * y(44) * y(56) + 0.8 * y(59)
 + 0.06 * y(42) *  (y(42) + y(44) + y(68) ) - 24 * y(45) - 6 * y(45)
 + 0.00499998 * y(41) *  (y(42) + y(44) + y(68) ) - 2 * y(46) - 0.5 * y(46)
 + 0.0199998 * y(43) *  (y(42) + y(44) + y(68) ) - 8 * y(47) - 2 * y(47)
 + 72 * y(23) * y(51)*y(51) - 72 * y(48) - 3.6 * y(48) * y(51) + 10 * y(49)
 + 3.6 * y(48) * y(51) - 10 * y(49) - 0.465012 * y(49) * y(51) + 10 * y(50)
 + 0.465012 * y(49) * y(51) - 10 * y(50) - 1000.02 * y(42) * y(50) + 0.1 * y(40) - 49.9998 * y(43) * y(50) - 540 * y(50) + 540 * y(1) + 5 * y(41)
 - 144 * y(23) * y(51)*y(51) + 144 * y(48) - 3.6 * y(48) * y(51) + 10 * y(49) - 0.465012 * y(49) * y(51) + 10 * y(50) + 100 * buff(3) - 100 * y(51)
 + 1 * y(55) - 5.99994 * y(52) * y(9) + 324 * y(11) + 3.06 * y(66) + 18 * y(33) + 0.306 * y(37)
 - 499.98 * y(53) * y(56) + 0.1 * y(54) + 81 * y(11) - 0.342004 * y(53) * y(63) + 12.24 * y(66) - 0.64584 * y(53) * y(30) + 75.0006 * y(33) - 0.0342004 * y(53) * y(35) + 1.224 * y(37)
 + 499.98 * y(53) * y(56) - 0.1 * y(54) - 0.342004 * y(54) * y(63) + 12.24 * y(67) - 0.64584 * y(54) * y(30) + 75.0006 * y(34)
 + 3.06 * y(67) + 18 * y(34) - 1 * y(55)
 - 0.499998 * y(40) * y(56) + 0.8 * y(57) + 0.2 * y(57) - 0.499998 * y(42) * y(56) + 0.8 * y(58) + 0.2 * y(58) - 0.499998 * y(44) * y(56) + 0.8 * y(61) + 0.2 * y(61) - 0.499998 * y(72) * y(56) + 0.8 * y(60) + 0.2 * y(60) - 0.499998 * y(44) * y(56) + 0.8 * y(59) + 0.2 * y(59) + 1 * y(55) - 499.98 * y(53) * y(56) + 0.1 * y(54)
 + 0.499998 * y(40) * y(56) - 0.8 * y(57) - 0.2 * y(57)
 + 0.499998 * y(42) * y(56) - 0.8 * y(58) - 0.2 * y(58)
 + 0.499998 * y(44) * y(56) - 0.8 * y(59) - 0.2 * y(59)
 + 0.499998 * y(72) * y(56) - 0.8 * y(60) - 0.2 * y(60)
 + 0.499998 * y(44) * y(56) - 0.8 * y(61) - 0.2 * y(61)
 - 0.00108 * y(14) * y(62) - 0.00108 * y(13) * y(62) + 0.3 * y(43) - 0.00054 * y(12) * y(62) + 0.3 * y(72) - 0.00054 * y(21) * y(62)
 - 0.341998 * y(27) * y(63) + 1.36 * y(64) + 0.34 * y(64) - 0.341998 * y(28) * y(63) + 1.36 * y(65) + 0.34 * y(65) - 0.0380004 * y(53) * y(63) + 1.36 * y(66) + 0.34 * y(66) - 0.0380004 * y(54) * y(63) + 1.36 * y(67) + 0.34 * y(67) + 599.994 * y(1) * y(35) - 1 * y(63)
 + 0.341998 * y(27) * y(63) - 1.36 * y(64) - 0.34 * y(64)
 + 0.341998 * y(28) * y(63) - 1.36 * y(65) - 0.34 * y(65)
 + 0.0380004 * y(53) * y(63) - 1.36 * y(66) - 0.34 * y(66)
 + 0.0380004 * y(54) * y(63) - 1.36 * y(67) - 0.34 * y(67)
 + 1 * buff(4) - 1 * y(68)
 + 0.09375 * y(42) *  (y(40) + y(41) ) - 24 * y(69) - 6 * y(69)
 + 0.0078126 * y(41) *  (y(40) + y(41) ) - 2 * y(70) - 0.5 * y(70)
 + 0.0312498 * y(43) *  (y(40) + y(41) ) - 8 * y(71) - 2 * y(71)
 - 0.3 * y(72) + 0.00054 * y(21) * y(62) - 0.499998 * y(72) * y(56) + 0.8 * y(60) + 0.2 * y(59)
];

% end of model
