/*
 * Modelling insulin kinetics in response to oral glucose administration
 * 
 * Model Status
 * 
 * This CellML model runs in both OpenCell and COR to recreate
 * the published results. Thank you to the model author Jorrit
 * Hornberg for providing us with the original mathematica code,
 * and also to Joe Cursons for the hours he has spent debugging
 * the CellMl model. Please note there are several differences
 * between the original CellML version of this model (which was
 * based on the published paper), and the current version of the
 * model (which was based on the Mathematica code). These differences
 * are listed here as a PDF document.
 * 
 * Model Structure
 * 
 * ABSTRACT:
 * 
 * ABSTRACT: Oncogenesis results from changes in kinetics or in
 * abundance of proteins in signal transduction networks. Recently,
 * it was shown that control of signalling cannot reside in a single
 * gene product, and might well be dispersed over many components.
 * Which of the reactions in these complex networks are most important,
 * and how can the existing molecular information be used to understand
 * why particular genes are oncogenes whereas others are not? We
 * implement a new method to help address such questions. We apply
 * control analysis to a detailed kinetic model of the epidermal
 * growth factor-induced mitogen-activated protein kinase network.
 * We determine the control of each reaction with respect to three
 * biologically relevant characteristics of the output of this
 * network: the amplitude, duration and integrated output of the
 * transient phosphorylation of extracellular signal-regulated
 * kinase (ERK). We confirm that control is distributed, but far
 * from randomly: a small proportion of reactions substantially
 * control signalling. In particular, the activity of Raf is in
 * control of all characteristics of the transient profile of ERK
 * phosphorylation, which may clarify why Raf is an oncogene. Most
 * reactions that really matter for one signalling characteristic
 * are also important for the other characteristics. Our analysis
 * also predicts the effects of mutations and changes in gene expression.
 * 
 * The original paper reference is cited below:
 * 
 * Control of MAPK signalling: from complexity to what really matters,
 * Jorrit J Hornberg, Bernd Binder, Frank J Bruggeman, Birgit Schoeberl,
 * Reinhart Heinrich and Hans V Westerhoff, 2005, Oncogene, 26,
 * 117-125. PubMed ID: 16007170
 * 
 * reaction diagram
 * 
 * [[Image file: hornberg_2005.png]]
 * 
 * Schematic diagram of the signal transduction nectwork described
 * by the model.
 */

import nsrunit;
unit conversion on;
// unit molar predefined
unit flux=1E3 meter^(-3)*second^(-1)*mole^1;
unit first_order_rate_constant=1 second^(-1);
unit second_order_rate_constant=.001 meter^3*second^(-1)*mole^(-1);

math main {
	realDomain time second;
	time.min=0;
	extern time.max;
	extern time.delta;
	real c1(time) molar;
	when(time=time.min) c1=5.0E-8;
	real v1(time) flux;
	real c2(time) molar;
	when(time=time.min) c2=5.0E4;
	real v6(time) flux;
	real v13(time) flux;
	real c3(time) molar;
	when(time=time.min) c3=0.0;
	real v2(time) flux;
	real c4(time) molar;
	when(time=time.min) c4=0.0;
	real v3(time) flux;
	real c5(time) molar;
	when(time=time.min) c5=0.0;
	real v7(time) flux;
	real v8(time) flux;
	real c6(time) molar;
	when(time=time.min) c6=0.0;
	real v10(time) flux;
	real v60(time) flux;
	real c7(time) molar;
	when(time=time.min) c7=0.0;
	real v4(time) flux;
	real v5(time) flux;
	real c8(time) molar;
	when(time=time.min) c8=0.0;
	real v12(time) flux;
	real v14(time) flux;
	real v62(time) flux;
	real c9(time) molar;
	when(time=time.min) c9=0.0;
	real v15(time) flux;
	real v107(time) flux;
	real v110(time) flux;
	real v113(time) flux;
	real v116(time) flux;
	real v119(time) flux;
	real v122(time) flux;
	real v125(time) flux;
	real c10(time) molar;
	when(time=time.min) c10=0.0;
	real v11(time) flux;
	real c11(time) molar;
	when(time=time.min) c11=0.0;
	real c12(time) molar;
	when(time=time.min) c12=8.1E4;
	real v106(time) flux;
	real v109(time) flux;
	real v112(time) flux;
	real v115(time) flux;
	real v118(time) flux;
	real v121(time) flux;
	real v124(time) flux;
	real c13(time) molar;
	when(time=time.min) c13=0.0;
	real v61(time) flux;
	real c14(time) molar;
	when(time=time.min) c14=1.2E4;
	real c15(time) molar;
	when(time=time.min) c15=0.0;
	real v16(time) flux;
	real v22(time) flux;
	real v32(time) flux;
	real v34(time) flux;
	real v37(time) flux;
	real v39(time) flux;
	real v102(time) flux;
	real c16(time) molar;
	when(time=time.min) c16=0.0;
	real c17(time) molar;
	when(time=time.min) c17=0.0;
	real v63(time) flux;
	real v69(time) flux;
	real v79(time) flux;
	real v80(time) flux;
	real v81(time) flux;
	real v82(time) flux;
	real v132(time) flux;
	real c18(time) molar;
	when(time=time.min) c18=0.0;
	real v9(time) flux;
	real v64(time) flux;
	real v133(time) flux;
	real c19(time) molar;
	when(time=time.min) c19=0.0;
	real v65(time) flux;
	real v66(time) flux;
	real v67(time) flux;
	real v68(time) flux;
	real v105(time) flux;
	real v127(time) flux;
	real v134(time) flux;
	real c20(time) molar;
	when(time=time.min) c20=0.0;
	real v108(time) flux;
	real v135(time) flux;
	real c21(time) molar;
	when(time=time.min) c21=0.0;
	real v111(time) flux;
	real v136(time) flux;
	real c22(time) molar;
	when(time=time.min) c22=1.1E4;
	real v24(time) flux;
	real v35(time) flux;
	real v38(time) flux;
	real v71(time) flux;
	real c23(time) molar;
	when(time=time.min) c23=0.0;
	real v17(time) flux;
	real c24(time) molar;
	when(time=time.min) c24=2.63E4;
	real v72(time) flux;
	real v130(time) flux;
	real v131(time) flux;
	real v25(time) flux;
	real v40(time) flux;
	real c25(time) molar;
	when(time=time.min) c25=0.0;
	real v18(time) flux;
	real v19(time) flux;
	real v20(time) flux;
	real v21(time) flux;
	real v126(time) flux;
	real c26(time) molar;
	when(time=time.min) c26=7.2E4;
	real v26(time) flux;
	real v31(time) flux;
	real v78(time) flux;
	real v73(time) flux;
	real c27(time) molar;
	when(time=time.min) c27=0.0;
	real c28(time) molar;
	when(time=time.min) c28=0.0;
	real v27(time) flux;
	real v28(time) flux;
	real c29(time) molar;
	when(time=time.min) c29=0.0;
	real c30(time) molar;
	when(time=time.min) c30=4.0E4;
	real v33(time) flux;
	real v41(time) flux;
	real v83(time) flux;
	real c31(time) molar;
	when(time=time.min) c31=1.01E5;
	real v36(time) flux;
	real c32(time) molar;
	when(time=time.min) c32=0.0;
	real v23(time) flux;
	real v103(time) flux;
	real c33(time) molar;
	when(time=time.min) c33=0.0;
	real v104(time) flux;
	real c34(time) molar;
	when(time=time.min) c34=0.0;
	real v114(time) flux;
	real c35(time) molar;
	when(time=time.min) c35=0.0;
	real v30(time) flux;
	real v117(time) flux;
	real v128(time) flux;
	real c36(time) molar;
	when(time=time.min) c36=0.0;
	real v120(time) flux;
	real c37(time) molar;
	when(time=time.min) c37=0.0;
	real v123(time) flux;
	real c38(time) molar;
	when(time=time.min) c38=0.0;
	real c39(time) molar;
	when(time=time.min) c39=0.0;
	real c40(time) molar;
	when(time=time.min) c40=0.0;
	real c41(time) molar;
	when(time=time.min) c41=4.0E4;
	real v43(time) flux;
	real v75(time) flux;
	real v85(time) flux;
	real c42(time) molar;
	when(time=time.min) c42=0.0;
	real v29(time) flux;
	real c43(time) molar;
	when(time=time.min) c43=0.0;
	real c44(time) molar;
	when(time=time.min) c44=4.0E4;
	real v42(time) flux;
	real v84(time) flux;
	real c45(time) molar;
	when(time=time.min) c45=0.0;
	real v44(time) flux;
	real v45(time) flux;
	real v46(time) flux;
	real v47(time) flux;
	real c46(time) molar;
	when(time=time.min) c46=0.0;
	real c47(time) molar;
	when(time=time.min) c47=2.1E7;
	real v51(time) flux;
	real v86(time) flux;
	real v93(time) flux;
	real c48(time) molar;
	when(time=time.min) c48=0.0;
	real c49(time) molar;
	when(time=time.min) c49=0.0;
	real v49(time) flux;
	real v50(time) flux;
	real c50(time) molar;
	when(time=time.min) c50=0.0;
	real c51(time) molar;
	when(time=time.min) c51=0.0;
	real v48(time) flux;
	real v52(time) flux;
	real v53(time) flux;
	real v54(time) flux;
	real v55(time) flux;
	real c52(time) molar;
	when(time=time.min) c52=0.0;
	real c53(time) molar;
	when(time=time.min) c53=4.0E4;
	real v90(time) flux;
	real v91(time) flux;
	real v92(time) flux;
	real c54(time) molar;
	when(time=time.min) c54=0.0;
	real c55(time) molar;
	when(time=time.min) c55=2.21E7;
	real v59(time) flux;
	real v94(time) flux;
	real v101(time) flux;
	real c56(time) molar;
	when(time=time.min) c56=0.0;
	real c57(time) molar;
	when(time=time.min) c57=0.0;
	real v57(time) flux;
	real v58(time) flux;
	real c58(time) molar;
	when(time=time.min) c58=0.0;
	real c59(time) molar;
	when(time=time.min) c59=0.0;
	real v56(time) flux;
	real v143(time) flux;
	real v144(time) flux;
	real v145(time) flux;
	real c60(time) molar;
	when(time=time.min) c60=1.0E7;
	real v98(time) flux;
	real v99(time) flux;
	real v100(time) flux;
	real c61(time) molar;
	when(time=time.min) c61=0.0;
	real c62(time) molar;
	when(time=time.min) c62=0.0;
	real c63(time) molar;
	when(time=time.min) c63=0.0;
	real v70(time) flux;
	real v137(time) flux;
	real c64(time) molar;
	when(time=time.min) c64=0.0;
	real v138(time) flux;
	real c65(time) molar;
	when(time=time.min) c65=0.0;
	real v139(time) flux;
	real c66(time) molar;
	when(time=time.min) c66=0.0;
	real v74(time) flux;
	real v77(time) flux;
	real v129(time) flux;
	real v140(time) flux;
	real c67(time) molar;
	when(time=time.min) c67=0.0;
	real v141(time) flux;
	real c68(time) molar;
	when(time=time.min) c68=0.0;
	real v142(time) flux;
	real c69(time) molar;
	when(time=time.min) c69=0.0;
	real c70(time) molar;
	when(time=time.min) c70=0.0;
	real v76(time) flux;
	real c71(time) molar;
	when(time=time.min) c71=0.0;
	real c72(time) molar;
	when(time=time.min) c72=0.0;
	real v87(time) flux;
	real v88(time) flux;
	real v89(time) flux;
	real c73(time) molar;
	when(time=time.min) c73=0.0;
	real c74(time) molar;
	when(time=time.min) c74=0.0;
	real c75(time) molar;
	when(time=time.min) c75=0.0;
	real c76(time) molar;
	when(time=time.min) c76=0.0;
	real c77(time) molar;
	when(time=time.min) c77=0.0;
	real v95(time) flux;
	real v96(time) flux;
	real v97(time) flux;
	real c78(time) molar;
	when(time=time.min) c78=0.0;
	real c79(time) molar;
	when(time=time.min) c79=0.0;
	real c80(time) molar;
	when(time=time.min) c80=0.0;
	real c81(time) molar;
	when(time=time.min) c81=0.0;
	real c82(time) molar;
	when(time=time.min) c82=0.0;
	real c83(time) molar;
	when(time=time.min) c83=0.0;
	real v146(time) flux;
	real v147(time) flux;
	real v148(time) flux;
	real c84(time) molar;
	when(time=time.min) c84=0.0;
	real c85(time) molar;
	when(time=time.min) c85=0.0;
	real c86(time) molar;
	when(time=time.min) c86=0.0;
	real c87(time) molar;
	when(time=time.min) c87=0.0;
	real c88(time) molar;
	when(time=time.min) c88=0.0;
	real c89(time) molar;
	when(time=time.min) c89=0.0;
	real c90(time) molar;
	when(time=time.min) c90=0.0;
	real c91(time) molar;
	when(time=time.min) c91=0.0;
	real c92(time) molar;
	when(time=time.min) c92=0.0;
	real c93(time) molar;
	when(time=time.min) c93=0.0;
	real c94(time) molar;
	when(time=time.min) c94=0.0;
	real c95(time) molar;
	when(time=time.min) c95=0.0;
	real c96(time) molar;
	when(time=time.min) c96=0.0;
	real c97(time) molar;
	when(time=time.min) c97=0.0;
	real c98(time) molar;
	when(time=time.min) c98=0.0;
	real c99(time) molar;
	when(time=time.min) c99=0.0;
	real c100(time) molar;
	when(time=time.min) c100=0.0;
	real c101(time) molar;
	when(time=time.min) c101=0.0;
	real c102(time) molar;
	when(time=time.min) c102=0.0;
	real c103(time) molar;
	when(time=time.min) c103=0.0;
	real k1 second_order_rate_constant;
	k1=3.0E7;
	real kd1 first_order_rate_constant;
	kd1=3.84E-3;
	real k2 second_order_rate_constant;
	k2=1.66E-5;
	real kd2 first_order_rate_constant;
	kd2=1E-1;
	real k3 first_order_rate_constant;
	k3=1E0;
	real kd3 first_order_rate_constant;
	kd3=1E-2;
	real k4 second_order_rate_constant;
	k4=1.73E-7;
	real kd4 first_order_rate_constant;
	kd4=1.66E-3;
	real k5 second_order_rate_constant;
	k5=0E0;
	real kd5 first_order_rate_constant;
	kd5=1.46E-2;
	real k6 first_order_rate_constant;
	k6=5E-4;
	real kd6 first_order_rate_constant;
	kd6=5E-3;
	real k8 second_order_rate_constant;
	k8=1.66E-6;
	real kd8 first_order_rate_constant;
	kd8=2E-1;
	real k10b second_order_rate_constant;
	k10b=5.43E-2;
	real kd10 first_order_rate_constant;
	kd10=1.1E-2;
	real k13 first_order_rate_constant;
	k13=2.17E0;
	real kd13 first_order_rate_constant;
	kd13=0E0;
	real k15 first_order_rate_constant;
	k15=1E4;
	real kd15 first_order_rate_constant;
	kd15=0E0;
	real k16 second_order_rate_constant;
	k16=1.66E-5;
	real kd63 first_order_rate_constant;
	kd63=2.75E-1;
	real k17 second_order_rate_constant;
	k17=1.66E-5;
	real kd17 first_order_rate_constant;
	kd17=6E-2;
	real k18 second_order_rate_constant;
	k18=2.5E-5;
	real kd18 first_order_rate_constant;
	kd18=1.3E0;
	real k19 second_order_rate_constant;
	k19=1.66E-7;
	real kd19 first_order_rate_constant;
	kd19=5E-1;
	real k20 second_order_rate_constant;
	k20=3.5E-6;
	real kd20 first_order_rate_constant;
	kd20=4E-1;
	real k21 second_order_rate_constant;
	k21=3.66E-7;
	real kd21 first_order_rate_constant;
	kd21=2.3E-2;
	real k22 second_order_rate_constant;
	k22=3.5E-5;
	real kd22 first_order_rate_constant;
	kd22=1E-1;
	real k23 first_order_rate_constant;
	k23=6E0;
	real kd23 first_order_rate_constant;
	kd23=6E-2;
	real kd24 first_order_rate_constant;
	kd24=5.5E-1;
	real k25 second_order_rate_constant;
	k25=1.66E-5;
	real kd25 first_order_rate_constant;
	kd25=2.14E-2;
	real k28 second_order_rate_constant;
	k28=1.66E-6;
	real kd28 first_order_rate_constant;
	kd28=5.3E-3;
	real k29 second_order_rate_constant;
	k29=1.17E-6;
	real kd29 first_order_rate_constant;
	kd29=1E0;
	real k32 second_order_rate_constant;
	k32=4E-7;
	real kd32 first_order_rate_constant;
	kd32=1E-1;
	real k33 second_order_rate_constant;
	k33=3.5E-5;
	real kd33 first_order_rate_constant;
	kd33=2E-1;
	real k34 second_order_rate_constant;
	k34=7.5E-6;
	real kd34 first_order_rate_constant;
	kd34=3E-2;
	real k35 second_order_rate_constant;
	k35=7.5E-6;
	real kd35 first_order_rate_constant;
	kd35=1.5E-3;
	real k36 first_order_rate_constant;
	k36=5E-3;
	real kd36 first_order_rate_constant;
	kd36=0E0;
	real k37 second_order_rate_constant;
	k37=1.5E-6;
	real kd37 first_order_rate_constant;
	kd37=3E-1;
	real k40 second_order_rate_constant;
	k40=5E-5;
	real kd40 first_order_rate_constant;
	kd40=6.4E-2;
	real k41 second_order_rate_constant;
	k41=5E-5;
	real kd41 first_order_rate_constant;
	kd41=4.29E-2;
	real k42 second_order_rate_constant;
	k42=1.18E-4;
	real kd42 first_order_rate_constant;
	kd42=2E-1;
	real k43 second_order_rate_constant;
	k43=0E0;
	real kd43 first_order_rate_constant;
	kd43=1E0;
	real k44 second_order_rate_constant;
	k44=1.95E-5;
	real kd52 first_order_rate_constant;
	kd52=3.3E-2;
	real k45 second_order_rate_constant;
	k45=0E0;
	real kd45 first_order_rate_constant;
	kd45=3.5E0;
	real k47 second_order_rate_constant;
	k47=0E0;
	real kd47 first_order_rate_constant;
	kd47=2.9E0;
	real k48 second_order_rate_constant;
	k48=2.38E-5;
	real kd48 first_order_rate_constant;
	kd48=8E-1;
	real k49 second_order_rate_constant;
	k49=0E0;
	real kd49 first_order_rate_constant;
	kd49=5.68E-2;
	real k50 second_order_rate_constant;
	k50=4.5E-7;
	real kd50 first_order_rate_constant;
	kd50=5E-1;
	real k52 second_order_rate_constant;
	k52=8.91E-5;
	real kd44 first_order_rate_constant;
	kd44=1.83E-2;
	real k53 second_order_rate_constant;
	k53=0E0;
	real kd53 first_order_rate_constant;
	kd53=1.6E1;
	real k55 second_order_rate_constant;
	k55=0E0;
	real kd55 first_order_rate_constant;
	kd55=5.7E0;
	real k56 second_order_rate_constant;
	k56=2.35E-5;
	real kd56 first_order_rate_constant;
	kd56=6E-1;
	real k57 second_order_rate_constant;
	k57=0E0;
	real kd57 first_order_rate_constant;
	kd57=2.46E-1;
	real k58 second_order_rate_constant;
	k58=8.33E-6;
	real kd58 first_order_rate_constant;
	kd58=5E-1;
	real k60 first_order_rate_constant;
	k60=5.5E-3;
	real kd60 first_order_rate_constant;
	kd60=0E0;
	real k61 first_order_rate_constant;
	k61=6.7E-4;
	real kd61 first_order_rate_constant;
	kd61=0E0;
	real k126 second_order_rate_constant;
	k126=1.66E-7;
	real kd126 first_order_rate_constant;
	kd126=2E0;
	real k127 second_order_rate_constant;
	k127=0E0;
	real kd127 first_order_rate_constant;
	kd127=1E-4;

	// <component name="environment">

	// <component name="c1">
	c1:time=(0 flux);

	// <component name="c2">
	c2:time=(v13-(v1+v6));

	// <component name="c3">
	c3:time=(v1-2*v2);

	// <component name="c4">
	c4:time=(v2-v3);

	// <component name="c5">
	c5:time=(v3-(v7+v8));

	// <component name="c6">
	c6:time=(v6-(v10+v60));

	// <component name="c7">
	c7:time=(v4+v5);

	// <component name="c8">
	c8:time=(v7+v12-(v14+v62));

	// <component name="c9">
	c9:time=((-1)*(v5+v15+v107+v110+v113+v116+v119+v122+v125));

	// <component name="c10">
	c10:time=(v10-2*v11);

	// <component name="c11">
	c11:time=(v11-v12);

	// <component name="c12">
	c12:time=(v15-(v4+v106+v109+v112+v115+v118+v121+v124));

	// <component name="c13">
	c13:time=v61;

	// <component name="c14">
	c14:time=((-1)*(v8+v14));

	// <component name="c15">
	c15:time=(v8-(v16+v22+v32+v34+v37+v39+v102));

	// <component name="c16">
	c16:time=((-1)*(v10+v61));

	// <component name="c17">
	c17:time=(v14+v102-(v63+v69+v79+v80+v81+v82+v132));

	// <component name="c18">
	c18:time=(v9+v63-(v5+v64+v133));

	// <component name="c19">
	c19:time=(v80+v64+v105-(v65+v66+v67+v68+v107+v127+v134));

	// <component name="c20">
	c20:time=(v65+v66+v108-(v110+v135));

	// <component name="c21">
	c21:time=(v67+v68+v111-(v113+v136));

	// <component name="c22">
	c22:time=((-1)*(v16+v24+v35+v38+v63+v71));

	// <component name="c23">
	c23:time=(v16-(v4+v9+v17));

	// <component name="c24">
	c24:time=((-1)*(v17+v25+v35+v40+v64+v72+v130+v131));

	// <component name="c25">
	c25:time=(v17+v34-(v18+v19+v20+v21+v105+v106+v126));

	// <component name="c26">
	c26:time=((-1)*(v18+v21+v26+v31+v65+v68+v73+v78));

	// <component name="c27">
	c27:time=(v18+v19-(v108+v109));

	// <component name="c28">
	c28:time=((-1)*(v19+v27+v28));

	// <component name="c29">
	c29:time=(v20+v21-(v111+v112));

	// <component name="c30">
	c30:time=(v35-(v33+v34+v41+v80+v83));

	// <component name="c31">
	c31:time=(v36-(v22+v69));

	// <component name="c32">
	c32:time=(v22-(v23+v103));

	// <component name="c33">
	c33:time=(v23+v37-(v24+v41+v104));

	// <component name="c34">
	c34:time=(v39+v24-(v25+v114+v115));

	// <component name="c35">
	c35:time=(v25+v32+v41-(v117+v118+v128+v26+v27+v30+v31));

	// <component name="c36">
	c36:time=(v27+v26-(v120+v121));

	// <component name="c37">
	c37:time=(v30+v31-(v123+v124));

	// <component name="c38">
	c38:time=(v33+v40-(v32+v79));

	// <component name="c39">
	c39:time=(v38-(v39+v40+v82));

	// <component name="c40">
	c40:time=((-1)*(v33+v36+v37+v38+v81));

	// <component name="c41">
	c41:time=((-1)*(v28+v43+v75+v85));

	// <component name="c42">
	c42:time=(v28+v29);

	// <component name="c43">
	c43:time=((-1)*(v29+v20+v30));

	// <component name="c44">
	c44:time=((-1)*(v42+v43+v84+v85));

	// <component name="c45">
	c45:time=((-1)*(v29+v42+v44+v45+v46+v47));

	// <component name="c46">
	c46:time=(v42+v43);

	// <component name="c47">
	c47:time=((-1)*(v44+v51+v86+v93));

	// <component name="c48">
	c48:time=(v44+v45);

	// <component name="c49">
	c49:time=((-1)*(v45+v46+v49+v50));

	// <component name="c50">
	c50:time=(v46+v47);

	// <component name="c51">
	c51:time=((-1)*(v47+v53+v55+v48+v52+v54));

	// <component name="c52">
	c52:time=(v48+v49);

	// <component name="c53">
	c53:time=((-1)*(v48+v49+v50+v51+v90+v91+v92+v93));

	// <component name="c54">
	c54:time=(v50+v51);

	// <component name="c55">
	c55:time=((-1)*(v52+v59+v94+v101));

	// <component name="c56">
	c56:time=(v52+v53);

	// <component name="c57">
	c57:time=((-1)*(v53+v54+v57+v58));

	// <component name="c58">
	c58:time=(v54+v55);

	// <component name="c59">
	c59:time=((-1)*(v55+v56+v126+v128+v130+v143+v144+v145));

	// <component name="c60">
	c60:time=((-1)*(v56+v57+v58+v59+v98+v99+v100+v101));

	// <component name="c61">
	c61:time=(v56+v57);

	// <component name="c62">
	c62:time=(v58+v59);

	// <component name="c63">
	c63:time=(v69+v103-(v70+v137));

	// <component name="c64">
	c64:time=(v70+v81+v104-(v71+v83+v138));

	// <component name="c65">
	c65:time=(v71+v82+v114-(v72+v116+v139));

	// <component name="c66">
	c66:time=(v72+v79+v83+v117-(v73+v74+v77+v78+v119+v129+v140));

	// <component name="c67">
	c67:time=(v73+v74+v120-(v122+v141));

	// <component name="c68">
	c68:time=(v77+v78+v123-(v125+v142));

	// <component name="c69">
	c69:time=((-1)*(v66+v74+v75));

	// <component name="c70">
	c70:time=(v75+v76);

	// <component name="c71">
	c71:time=((-1)*(v67+v76+v77));

	// <component name="c72">
	c72:time=((-1)*(v76+v84+v86+v87+v88+v89));

	// <component name="c73">
	c73:time=(v84+v85);

	// <component name="c74">
	c74:time=(v86+v87);

	// <component name="c75">
	c75:time=((-1)*(v87+v88+v91+v92));

	// <component name="c76">
	c76:time=(v88+v89);

	// <component name="c77">
	c77:time=((-1)*(v89+v90+v94+v95+v96+v97));

	// <component name="c78">
	c78:time=(v90+v91);

	// <component name="c79">
	c79:time=(v92+v93);

	// <component name="c80">
	c80:time=(v94+v95);

	// <component name="c81">
	c81:time=((-1)*(v95+v96+v99+v100));

	// <component name="c82">
	c82:time=(v96+v97);

	// <component name="c83">
	c83:time=((-1)*(v97+v98+v127+v129+v131+v146+v147+v148));

	// <component name="c84">
	c84:time=(v98+v99);

	// <component name="c85">
	c85:time=(v100+v101);

	// <component name="c86">
	c86:time=v60;

	// <component name="c87">
	c87:time=(v62+v132+v133+v134+v135+v136+v137+v138+v139+v140+v141+v142);

	// <component name="c88">
	c88:time=(v106+v107);

	// <component name="c89">
	c89:time=(v109+v110);

	// <component name="c90">
	c90:time=(v112+v113);

	// <component name="c91">
	c91:time=(v115+v116);

	// <component name="c92">
	c92:time=(v118+v119);

	// <component name="c93">
	c93:time=(v121+v122);

	// <component name="c94">
	c94:time=(v124+v125);

	// <component name="c95">
	c95:time=(v126+v143);

	// <component name="c96">
	c96:time=(v127+v146);

	// <component name="c97">
	c97:time=(v128+v144);

	// <component name="c98">
	c98:time=(v129+v147);

	// <component name="c99">
	c99:time=((-1)*(v143+v144));

	// <component name="c100">
	c100:time=((-1)*(v146+v147));

	// <component name="c101">
	c101:time=(v130+v145);

	// <component name="c102">
	c102:time=(v131+v148);

	// <component name="c103">
	c103:time=((-1)*(v145+v148));

	// <component name="v1">
	v1=(k1*c1*c2-kd1*c3);

	// <component name="v2">
	v2=(k2*c3*c3-kd2*c4);

	// <component name="v3">
	v3=(k3*c4*1-kd3*c5);

	// <component name="v4">
	v4=(k4*c23*c12-kd4*c7);

	// <component name="v5">
	v5=(k5*c18*c9-kd5*c7);

	// <component name="v6">
	v6=(k6*c2-kd6*c6);

	// <component name="v7">
	v7=(k6*c5-kd6*c8);

	// <component name="v8">
	v8=(k8*c5*c14-kd8*c15);

	// <component name="v9">
	v9=(k6*c23-kd6*c18);

	// <component name="v10">
	v10=(k10b*c6*c16-kd10*c10);

	// <component name="v11">
	v11=(k2*c10*c10-kd2*c11);

	// <component name="v12">
	v12=(k3*c11-kd3*c8);

	// <component name="v13">
	v13=(k13*(1 molar)-kd13*c2);

	// <component name="v14">
	v14=(k8*c8*c14-kd8*c17);

	// <component name="v15">
	v15=(k15*c9-kd15*c12);

	// <component name="v16">
	v16=(k16*c22*c15-kd63*c23);

	// <component name="v17">
	v17=(k17*c24*c23-kd17*c25);

	// <component name="v18">
	v18=(k18*c26*c25-kd18*c27);

	// <component name="v19">
	v19=(k19*c28*c25-kd19*c27);

	// <component name="v20">
	v20=(k20*c25*c43-kd20*c29);

	// <component name="v21">
	v21=(k21*c25*c26-kd21*c29);

	// <component name="v22">
	v22=(k22*c31*c15-kd22*c32);

	// <component name="v23">
	v23=(k23*c32-kd23*c33);

	// <component name="v24">
	v24=(k16*c22*c33-kd24*c34);

	// <component name="v25">
	v25=(k25*c24*c34-kd25*c35);

	// <component name="v26">
	v26=(k18*c26*c35-kd18*c36);

	// <component name="v27">
	v27=(k19*c35*c28-kd19*c36);

	// <component name="v28">
	v28=(k28*c28*c41-kd28*c42);

	// <component name="v29">
	v29=(k29*c43*c45-kd29*c42);

	// <component name="v30">
	v30=(k20*c35*c43-kd20*c37);

	// <component name="v31">
	v31=(k21*c35*c26-kd21*c37);

	// <component name="v32">
	v32=(k32*c38*c15-kd32*c35);

	// <component name="v33">
	v33=(k33*c40*c30-kd33*c38);

	// <component name="v34">
	v34=(k34*c15*c30-kd34*c25);

	// <component name="v35">
	v35=(k35*c24*c22-kd35*c30);

	// <component name="v36">
	v36=(k36*c40-kd36*c31);

	// <component name="v37">
	v37=(k37*c15*c40-kd37*c33);

	// <component name="v38">
	v38=(k16*c22*c40-kd24*c39);

	// <component name="v39">
	v39=(k37*c15*c39-kd37*c34);

	// <component name="v40">
	v40=(k40*c24*c39-kd40*c38);

	// <component name="v41">
	v41=(k41*c30*c33-kd41*c35);

	// <component name="v42">
	v42=(k42*c44*c45-kd42*c46);

	// <component name="v43">
	v43=(k43*c41*c44-kd43*c46);

	// <component name="v44">
	v44=(k44*c47*c45-kd52*c48);

	// <component name="v45">
	v45=(k45*c49*c45-kd45*c48);

	// <component name="v46">
	v46=(k44*c49*c45-kd52*c50);

	// <component name="v47">
	v47=(k47*c51*c45-kd47*c50);

	// <component name="v48">
	v48=(k48*c51*c53-kd48*c52);

	// <component name="v49">
	v49=(k49*c49*c53-kd49*c52);

	// <component name="v50">
	v50=(k50*c53*c49-kd50*c54);

	// <component name="v51">
	v51=(k49*c47*c53-kd49*c54);

	// <component name="v52">
	v52=(k52*c55*c51-kd44*c56);

	// <component name="v53">
	v53=(k53*c51*c57-kd53*c56);

	// <component name="v54">
	v54=(k52*c51*c57-kd44*c58);

	// <component name="v55">
	v55=(k55*c59*c51-kd55*c58);

	// <component name="v56">
	v56=(k56*c59*c60-kd56*c61);

	// <component name="v57">
	v57=(k57*c57*c60-kd57*c61);

	// <component name="v58">
	v58=(k58*c60*c57-kd58*c62);

	// <component name="v59">
	v59=(k57*c55*c60-kd57*c62);

	// <component name="v60">
	v60=(k60*c6-kd60*c86);

	// <component name="v61">
	v61=(k61*c16-kd61*c13);

	// <component name="v62">
	v62=(k60*c8-kd60*c87);

	// <component name="v63">
	v63=(k16*c17*c22-kd63*c18);

	// <component name="v64">
	v64=(k17*c24*c18-kd17*c19);

	// <component name="v65">
	v65=(k18*c26*c19-kd18*c20);

	// <component name="v66">
	v66=(k19*c69*c19-kd19*c20);

	// <component name="v67">
	v67=(k20*c71*c19-kd20*c21);

	// <component name="v68">
	v68=(k21*c19*c26-kd21*c21);

	// <component name="v69">
	v69=(k22*c31*c17-kd22*c63);

	// <component name="v70">
	v70=(k23*c63-kd23*c64);

	// <component name="v71">
	v71=(k16*c22*c64-kd24*c65);

	// <component name="v72">
	v72=(k25*c24*c65-kd25*c66);

	// <component name="v73">
	v73=(k18*c26*c66-kd18*c67);

	// <component name="v74">
	v74=(k19*c66*c69-kd19*c67);

	// <component name="v75">
	v75=(k28*c69*c41-kd28*c70);

	// <component name="v76">
	v76=(k29*c71*c72-kd29*c70);

	// <component name="v77">
	v77=(k20*c71*c66-kd20*c68);

	// <component name="v78">
	v78=(k21*c66*c26-kd21*c68);

	// <component name="v79">
	v79=(k32*c17*c38-kd32*c66);

	// <component name="v80">
	v80=(k34*c17*c30-kd34*c19);

	// <component name="v81">
	v81=(k37*c17*c40-kd37*c64);

	// <component name="v82">
	v82=(k37*c17*c39-kd37*c65);

	// <component name="v83">
	v83=(k41*c30*c64-kd41*c66);

	// <component name="v84">
	v84=(k42*c44*c72-kd42*c73);

	// <component name="v85">
	v85=(k43*c41*c44-kd43*c73);

	// <component name="v86">
	v86=(k44*c47*c72-kd52*c74);

	// <component name="v87">
	v87=(k45*c75*c72-kd45*c74);

	// <component name="v88">
	v88=(k44*c72*c75-kd52*c76);

	// <component name="v89">
	v89=(k47*c72*c77-kd47*c76);

	// <component name="v90">
	v90=(k48*c77*c53-kd48*c78);

	// <component name="v91">
	v91=(k49*c75*c53-kd49*c78);

	// <component name="v92">
	v92=(k50*c53*c75-kd50*c79);

	// <component name="v93">
	v93=(k49*c47*c53-kd49*c79);

	// <component name="v94">
	v94=(k52*c55*c77-kd44*c80);

	// <component name="v95">
	v95=(k53*c81*c77-kd53*c80);

	// <component name="v96">
	v96=(k52*c77*c81-kd44*c82);

	// <component name="v97">
	v97=(k55*c83*c77-kd55*c82);

	// <component name="v98">
	v98=(k56*c83*c60-kd56*c84);

	// <component name="v99">
	v99=(k57*c81*c60-kd57*c84);

	// <component name="v100">
	v100=(k58*c60*c81-kd58*c85);

	// <component name="v101">
	v101=(k57*c55*c60-kd57*c85);

	// <component name="v102">
	v102=(k6*c15-kd6*c17);

	// <component name="v103">
	v103=(k6*c32-kd6*c63);

	// <component name="v104">
	v104=(k6*c33-kd6*c64);

	// <component name="v105">
	v105=(k6*c25-kd6*c19);

	// <component name="v106">
	v106=(k4*c25*c12-kd4*c88);

	// <component name="v107">
	v107=(k5*c9*c19-kd5*c88);

	// <component name="v108">
	v108=(k6*c27-kd6*c20);

	// <component name="v109">
	v109=(k4*c27*c12-kd4*c89);

	// <component name="v110">
	v110=(k5*c9*c20-kd5*c89);

	// <component name="v111">
	v111=(k6*c29-kd6*c21);

	// <component name="v112">
	v112=(k4*c29*c12-kd4*c90);

	// <component name="v113">
	v113=(k5*c9*c21-kd5*c90);

	// <component name="v114">
	v114=(k6*c34-kd6*c65);

	// <component name="v115">
	v115=(k4*c34*c12-kd4*c91);

	// <component name="v116">
	v116=(k5*c9*c65-kd5*c91);

	// <component name="v117">
	v117=(k6*c35-kd6*c66);

	// <component name="v118">
	v118=(k4*c35*c12-kd4*c92);

	// <component name="v119">
	v119=(k5*c9*c66-kd5*c92);

	// <component name="v120">
	v120=(k6*c36-kd6*c67);

	// <component name="v121">
	v121=(k4*c36*c12-kd4*c93);

	// <component name="v122">
	v122=(k5*c9*c67-kd5*c93);

	// <component name="v123">
	v123=(k6*c37-kd6*c68);

	// <component name="v124">
	v124=(k4*c37*c12-kd4*c94);

	// <component name="v125">
	v125=(k5*c68*c9-kd5*c94);

	// <component name="v126">
	v126=(k126*c59*c25-kd126*c95);

	// <component name="v127">
	v127=(k126*c83*c19-kd126*c96);

	// <component name="v128">
	v128=(k126*c59*c35-kd126*c97);

	// <component name="v129">
	v129=(k126*c83*c66-kd126*c98);

	// <component name="v130">
	v130=(k126*c59*c24-kd126*c101);

	// <component name="v131">
	v131=(k126*c83*c24-kd126*c102);

	// <component name="v132">
	v132=(k60*c17-kd60*c87);

	// <component name="v133">
	v133=(k60*c18-kd60*c87);

	// <component name="v134">
	v134=(k60*c19-kd60*c87);

	// <component name="v135">
	v135=(k60*c20-kd60*c87);

	// <component name="v136">
	v136=(k60*c21-kd60*c87);

	// <component name="v137">
	v137=(k60*c63-kd60*c87);

	// <component name="v138">
	v138=(k60*c64-kd60*c87);

	// <component name="v139">
	v139=(k60*c65-kd60*c87);

	// <component name="v140">
	v140=(k60*c66-kd60*c87);

	// <component name="v141">
	v141=(k60*c67-kd60*c87);

	// <component name="v142">
	v142=(k60*c68-kd60*c87);

	// <component name="v143">
	v143=(k127*c59*c99-kd127*c95);

	// <component name="v144">
	v144=(k127*c59*c99-kd127*c97);

	// <component name="v145">
	v145=(k127*c59*c103-kd127*c101);

	// <component name="v146">
	v146=(k127*c83*c100-kd127*c96);

	// <component name="v147">
	v147=(k127*c83*c100-kd127*c98);

	// <component name="v148">
	v148=(k127*c83*c103-kd127*c102);

	// <component name="reaction_rates">
}