← Back to branch summary

~ubuntu-branches/ubuntu/maverick/dolfin/maverick

« back to all changes in this revision

Viewing changes to demo/ode/tentusscher/cpp/tentusscher.h

  • Committer: Bazaar Package Importer
  • Author(s): Johannes Ring
  • Date: 2008-09-16 08:41:20 UTC
  • Revision ID: james.westby@ubuntu.com-20080916084120-i8k3u6lhx3mw3py3
Tags: upstream-0.9.2
Import upstream version 0.9.2

Show diffs side-by-side

added added

removed removed

Lines of Context:
demo/ode/tentusscher/cpp/tentusscher.h
 
1
#include <dolfin.h>
 
2
 
 
3
class tentusscher: public dolfin::ODE
 
4
{
 
5
  
 
6
    protected:
 
7
        const double GNa;
 
8
        const double R;
 
9
        const double F;
 
10
        const double T;
 
11
        const double RTonF;
 
12
        const double Nao;
 
13
        const double GK1;
 
14
        const double Ko;
 
15
        const double GKr;
 
16
        const double pKNa;
 
17
        const double GCaL;
 
18
        const double Cao;
 
19
        const double kNaCa;
 
20
        const double gamma;
 
21
        const double alpha;
 
22
        const double KmNai;
 
23
        const double KmCa;
 
24
        const double ksat;
 
25
        const double PNaK;
 
26
        const double KmK;
 
27
        const double KmNa;
 
28
        const double GpCa;
 
29
        const double KpCa;
 
30
        const double GpK;
 
31
        const double GbCa;
 
32
        const double GbNa;
 
33
        const double Cm;
 
34
        const double km2;
 
35
        const double maxSR;
 
36
        const double minSR;
 
37
        const double EC;
 
38
        const double k4;
 
39
        const double Vc;
 
40
        const double Vsr;
 
41
        const double Vleak;
 
42
        const double Vmaxup;
 
43
        const double Kup;
 
44
        const double Vxfer;
 
45
        const double Vrel;
 
46
        const double k3;
 
47
        const double km1;
 
48
        const double Vss;
 
49
        const double Bufc;
 
50
        const double Kbufc;
 
51
        const double Bufss;
 
52
        const double Kbufss;
 
53
        const double Bufsr;
 
54
        const double Kbufsr;
 
55
        const double Gto;
 
56
        const double GKs;
 
57
        double ist;
 
58
        double dt;
 
59
        
 
60
    public:
 
61
        tentusscher(double t_stop) : dolfin::ODE(19, t_stop), GNa(14.838), R(8314.472), F(96485.3415), T(310.0), RTonF(26.7137606597), Nao(140.0), GK1(5.405), Ko(5.4), GKr(0.153), pKNa(0.03), GCaL(3.98e-05), Cao(2.0), kNaCa(1000.0), gamma(0.35), alpha(2.5), KmNai(87.5), KmCa(1.38), ksat(0.1), PNaK(2.724), KmK(1.0), KmNa(40.0), GpCa(0.1238), KpCa(0.0005), GpK(0.0146), GbCa(0.000592), GbNa(0.00029), Cm(0.185), km2(0.045), maxSR(2.5), minSR(1.0), EC(1.5), k4(0.005), Vc(0.016404), Vsr(0.001094), Vleak(0.00036), Vmaxup(0.006375), Kup(0.00025), Vxfer(0.0038), Vrel(0.102), k3(0.06), km1(0.15), Vss(5.468e-05), Bufc(0.2), Kbufc(0.001), Bufss(0.4), Kbufss(0.00025), Bufsr(10.0), Kbufsr(0.3), Gto(0.294), GKs(0.392) 
 
62
        {
 
63
            ist = 0.0;
 
64
            dt = 0.1;
 
65
        }
 
66
 
 
67
    ~tentusscher() {}
 
68
 
 
69
        void f(const double* u, double t, double* y) {
 
70
            y[0] = -ist+4.6844433675364579e-04*1.0/( 1.0000000000000001e-01*exp(-2.4332029034784559e-02*u[0])+1.0)*( -2.0000000000000000e+00*(u[5]*u[5]*u[5])*exp(1.3101861787960915e-02*u[0])+6.8600000000000000e+06*u[2]*exp(-2.4332029034784559e-02*u[0]))+-8.0205405405405401e+01*u[8]*u[9]*( -2.6713760659695648e+01*log(1.4000000000000000e+02*1.0/u[5])+u[0])*(u[7]*u[7]*u[7])+-2.9216216216216218e+00*1.0/( ( exp( 1.0000000000000001e-01*u[0]+-2.6713760659695653e+00*log(5.4000000000000004e+00*1.0/u[6])-1.0000000000000000e+00)+3.0*exp( 2.0000000000000001e-04*u[0]+-5.3427521319391299e-03*log(5.4000000000000004e+00*1.0/u[6])+2.0000000000000000e-02))/( exp( -5.0000000000000000e-01*u[0]+1.3356880329847824e+01*log(5.4000000000000004e+00*1.0/u[6]))+1.0)+1.0000000000000001e-01*1.0/( exp( 5.9999999999999998e-02*u[0]+-1.6028256395817388e+00*log(5.4000000000000004e+00*1.0/u[6])-1.2000000000000000e+01)+1.0))/( exp( 5.9999999999999998e-02*u[0]+-1.6028256395817388e+00*log(5.4000000000000004e+00*1.0/u[6])-1.2000000000000000e+01)+1.0)*( u[0]+-2.6713760659695648e+01*log(5.4000000000000004e+00*1.0/u[6]))+-8.2702702702702702e-01*u[12]*u[13]*( u[0]+-2.6713760659695648e+01*log(5.4000000000000004e+00*1.0/u[6]))+-6.6918918918918913e-01*u[2]/( u[2]+5.0000000000000001e-04)+-1.5891891891891892e+00*u[11]*u[10]*( u[0]+-2.6713760659695648e+01*log(5.4000000000000004e+00*1.0/u[6]))+3.1081190322229455e+00*( u[0]-15.0)/( exp( 7.4867781645490947e-02*u[0]-1.1230167246823641e+00)-1.0)*u[17]*( -2.5000000000000000e-01*exp( 7.4867781645490947e-02*u[0]-1.1230167246823641e+00)*u[4]+2.0000000000000000e+00)*u[18]*u[15]*u[16]+4.1875624817901289e-02*log(1.4000000000000000e+02*1.0/u[5])+-1.2423648648648649e+01*u[5]/( u[5]+4.0000000000000000e+01)/( 3.5299999999999998e-02*exp(-3.7433890822745473e-02*u[0])+1.2450000000000000e-01*exp(-3.7433890822745472e-03*u[0])+1.0)+-2.1189189189189190e+00*(u[14]*u[14])*( u[0]+-2.6713760659695648e+01*log(9.5999999999999996e+00*1.0/( 2.9999999999999999e-02*u[5]+u[6])))+-4.7675675675675671e-03*u[0]+4.2742017055513039e-02*log(2.0000000000000000e+00*1.0/u[2])+-7.8918918918918918e-02*1.0/( exp( -1.6722408026755853e-01*u[0]+4.1806020066889626e+00)+1.0)*( u[0]+-2.6713760659695648e+01*log(5.4000000000000004e+00*1.0/u[6]));
 
71
            y[1] =  -5.0000000000000001e-03*u[1]+-4.4999999999999998e-02*( -1.5000000000000000e+00*1.0/( 2.2500000000000000e+00*1.0/(u[3]*u[3])+1.0)+2.5000000000000000e+00)*u[4]*u[1]+5.0000000000000001e-03;
 
72
            y[2] = -(1.0/8.0)*1.0/dt*( -8.0000000000000004e-01*1.0/( u[2]+1.0000000000000000e-03)*u[2]+dt*( 1.5296035113386978e-02*u[2]+2.8940849167139207e-05*1.0/( u[2]+5.0000000000000001e-04)*u[2]+4.0518218501449887e-08*( -2.0000000000000000e+00*(u[5]*u[5]*u[5])*exp(1.3101861787960915e-02*u[0])+6.8600000000000000e+06*u[2]*exp(-2.4332029034784559e-02*u[0]))/( 1.0000000000000001e-01*exp(-2.4332029034784559e-02*u[0])+1.0)+-9.6035113386978800e-05*u[3]+1.3839242897371898e-07*u[0]+1.7006217995610825e-03*1.0/( 6.2499999999999997e-08*1.0/(u[2]*u[2])+1.0)+-1.8484911123589293e-06*log(2.0000000000000000e+00*1.0/u[2])+-1.5200000000000000e-02*u[4])+4.0*u[2]+-4.0*pow( 8.0000000000000004e-04*1.0/( u[2]+1.0000000000000000e-03)*u[2]+-1.0000000000000000e-03*dt*( 1.5296035113386978e-02*u[2]+2.8940849167139207e-05*1.0/( u[2]+5.0000000000000001e-04)*u[2]+4.0518218501449887e-08*( -2.0000000000000000e+00*(u[5]*u[5]*u[5])*exp(1.3101861787960915e-02*u[0])+6.8600000000000000e+06*u[2]*exp(-2.4332029034784559e-02*u[0]))/( 1.0000000000000001e-01*exp(-2.4332029034784559e-02*u[0])+1.0)+-9.6035113386978800e-05*u[3]+1.3839242897371898e-07*u[0]+1.7006217995610825e-03*1.0/( 6.2499999999999997e-08*1.0/(u[2]*u[2])+1.0)+-1.8484911123589293e-06*log(2.0000000000000000e+00*1.0/u[2])+-1.5200000000000000e-02*u[4])+pow( -8.0000000000000004e-01*1.0/( u[2]+1.0000000000000000e-03)*u[2]+dt*( 1.5296035113386978e-02*u[2]+2.8940849167139207e-05*1.0/( u[2]+5.0000000000000001e-04)*u[2]+4.0518218501449887e-08*( -2.0000000000000000e+00*(u[5]*u[5]*u[5])*exp(1.3101861787960915e-02*u[0])+6.8600000000000000e+06*u[2]*exp(-2.4332029034784559e-02*u[0]))/( 1.0000000000000001e-01*exp(-2.4332029034784559e-02*u[0])+1.0)+-9.6035113386978800e-05*u[3]+1.3839242897371898e-07*u[0]+1.7006217995610825e-03*1.0/( 6.2499999999999997e-08*1.0/(u[2]*u[2])+1.0)+-1.8484911123589293e-06*log(2.0000000000000000e+00*1.0/u[2])+-1.5200000000000000e-02*u[4])+-4.0*u[2]+8.0400000000000005e-01,2.0)/16.0+4.0000000000000001e-03*u[2],(1.0/2.0))+8.0400000000000005e-01);
 
73
            y[3] = (1.0/2.0)*( 1.0000000000000000e+01*1.0/( u[3]+2.9999999999999999e-01)*u[3]+dt*( 6.3749999999999996e-03*1.0/( 6.2499999999999997e-08*1.0/(u[2]*u[2])+1.0)+3.6000000000000002e-04*u[2]+-3.6000000000000002e-04*u[3]+1.5299999999999998e-02*(u[4]*u[4])*u[1]*( u[4]-u[3])/( -1.5000000000000000e+00*1.0/( 2.2500000000000000e+00*1.0/(u[3]*u[3])+1.0)+2.5000000000000000e+00)/( 1.4999999999999999e-01*(u[4]*u[4])/( -1.5000000000000000e+00*1.0/( 2.2500000000000000e+00*1.0/(u[3]*u[3])+1.0)+2.5000000000000000e+00)+5.9999999999999998e-02))+pow( 1.2000000000000000e+01*1.0/( u[3]+2.9999999999999999e-01)*u[3]+1.2000000000000000e+00*dt*( 6.3749999999999996e-03*1.0/( 6.2499999999999997e-08*1.0/(u[2]*u[2])+1.0)+3.6000000000000002e-04*u[2]+-3.6000000000000002e-04*u[3]+1.5299999999999998e-02*(u[4]*u[4])*u[1]*( u[4]-u[3])/( -1.5000000000000000e+00*1.0/( 2.2500000000000000e+00*1.0/(u[3]*u[3])+1.0)+2.5000000000000000e+00)/( 1.4999999999999999e-01*(u[4]*u[4])/( -1.5000000000000000e+00*1.0/( 2.2500000000000000e+00*1.0/(u[3]*u[3])+1.0)+2.5000000000000000e+00)+5.9999999999999998e-02))+1.2000000000000000e+00*u[3]+pow( 1.0000000000000000e+01*1.0/( u[3]+2.9999999999999999e-01)*u[3]+dt*( 6.3749999999999996e-03*1.0/( 6.2499999999999997e-08*1.0/(u[2]*u[2])+1.0)+3.6000000000000002e-04*u[2]+-3.6000000000000002e-04*u[3]+1.5299999999999998e-02*(u[4]*u[4])*u[1]*( u[4]-u[3])/( -1.5000000000000000e+00*1.0/( 2.2500000000000000e+00*1.0/(u[3]*u[3])+1.0)+2.5000000000000000e+00)/( 1.4999999999999999e-01*(u[4]*u[4])/( -1.5000000000000000e+00*1.0/( 2.2500000000000000e+00*1.0/(u[3]*u[3])+1.0)+2.5000000000000000e+00)+5.9999999999999998e-02))+u[3]-1.0300000000000001e+01,2.0),(1.0/2.0))-u[3]-1.0300000000000001e+01)/dt;
 
74
            y[4] = (1.0/2.0)*( 4.0000000000000002e-01*u[4]/( u[4]+2.5000000000000001e-04)-dt*( 3.0611192392099484e-01*1.0/( 1.4999999999999999e-01*(u[4]*u[4])/( -1.5000000000000000e+00*1.0/( 2.2500000000000000e+00*1.0/(u[3]*u[3])+1.0)+2.5000000000000000e+00)+5.9999999999999998e-02)*(u[4]*u[4])*u[1]*( u[4]-u[3])/( -1.5000000000000000e+00*1.0/( 2.2500000000000000e+00*1.0/(u[3]*u[3])+1.0)+2.5000000000000000e+00)+1.1399999999999999e+00*u[4]+-1.0081409587705739e-02*u[17]*( u[0]-15.0)*u[18]/( exp( 7.4867781645490947e-02*u[0]-1.1230167246823641e+00)-1.0)*u[15]*( -2.5000000000000000e-01*u[4]*exp( 7.4867781645490947e-02*u[0]-1.1230167246823641e+00)+2.0000000000000000e+00)*u[16]+-1.1399999999999999e+00*u[2])-u[4]+pow( 4.0000000000000002e-04*u[4]/( u[4]+2.5000000000000001e-04)+-1.0000000000000000e-03*dt*( 3.0611192392099484e-01*1.0/( 1.4999999999999999e-01*(u[4]*u[4])/( -1.5000000000000000e+00*1.0/( 2.2500000000000000e+00*1.0/(u[3]*u[3])+1.0)+2.5000000000000000e+00)+5.9999999999999998e-02)*(u[4]*u[4])*u[1]*( u[4]-u[3])/( -1.5000000000000000e+00*1.0/( 2.2500000000000000e+00*1.0/(u[3]*u[3])+1.0)+2.5000000000000000e+00)+1.1399999999999999e+00*u[4]+-1.0081409587705739e-02*u[17]*( u[0]-15.0)*u[18]/( exp( 7.4867781645490947e-02*u[0]-1.1230167246823641e+00)-1.0)*u[15]*( -2.5000000000000000e-01*u[4]*exp( 7.4867781645490947e-02*u[0]-1.1230167246823641e+00)+2.0000000000000000e+00)*u[16]+-1.1399999999999999e+00*u[2])+1.0000000000000000e-03*u[4]+pow( 4.0000000000000002e-01*u[4]/( u[4]+2.5000000000000001e-04)-dt*( 3.0611192392099484e-01*1.0/( 1.4999999999999999e-01*(u[4]*u[4])/( -1.5000000000000000e+00*1.0/( 2.2500000000000000e+00*1.0/(u[3]*u[3])+1.0)+2.5000000000000000e+00)+5.9999999999999998e-02)*(u[4]*u[4])*u[1]*( u[4]-u[3])/( -1.5000000000000000e+00*1.0/( 2.2500000000000000e+00*1.0/(u[3]*u[3])+1.0)+2.5000000000000000e+00)+1.1399999999999999e+00*u[4]+-1.0081409587705739e-02*u[17]*( u[0]-15.0)*u[18]/( exp( 7.4867781645490947e-02*u[0]-1.1230167246823641e+00)-1.0)*u[15]*( -2.5000000000000000e-01*u[4]*exp( 7.4867781645490947e-02*u[0]-1.1230167246823641e+00)+2.0000000000000000e+00)*u[16]+-1.1399999999999999e+00*u[2])+u[4]-4.0024999999999999e-01,2.0),(1.0/2.0))-4.0024999999999999e-01)/dt;
 
75
            y[5] =  -1.7343470110743600e-03*( -2.6713760659695648e+01*log(1.4000000000000000e+02*1.0/u[5])+u[0])*u[8]*u[9]*(u[7]*u[7]*u[7])+3.0388663876087418e-08*( -2.0000000000000000e+00*exp(1.3101861787960915e-02*u[0])*(u[5]*u[5]*u[5])+6.8600000000000000e+06*u[2]*exp(-2.4332029034784559e-02*u[0]))/( 1.0000000000000001e-01*exp(-2.4332029034784559e-02*u[0])+1.0)+9.0551084895961074e-07*log(1.4000000000000000e+02*1.0/u[5])+-3.3896794258765630e-08*u[0]+-8.0594011556369446e-04*u[5]/( u[5]+4.0000000000000000e+01)/( 1.2450000000000000e-01*exp(-3.7433890822745472e-03*u[0])+3.5299999999999998e-02*exp(-3.7433890822745473e-02*u[0])+1.0);
 
76
            y[6] =  1.7065282626826835e-06*1.0/( exp( -1.6722408026755853e-01*u[0]+4.1806020066889626e+00)+1.0)*( 2.6713760659695648e+01*log(5.4000000000000004e+00*1.0/u[6])-u[0])+-1.1688549744401942e-04*ist+1.7883481108934970e-05*u[12]*( 2.6713760659695648e+01*log(5.4000000000000004e+00*1.0/u[6])-u[0])*u[13]+6.3176611368492499e-05*( 2.6713760659695648e+01*log(5.4000000000000004e+00*1.0/u[6])-u[0])/( 1.0/( exp( 1.3356880329847824e+01*log(5.4000000000000004e+00*1.0/u[6])+-5.0000000000000000e-01*u[0])+1.0)*( exp( -2.6713760659695653e+00*log(5.4000000000000004e+00*1.0/u[6])+1.0000000000000001e-01*u[0]-1.0000000000000000e+00)+3.0*exp( -5.3427521319391299e-03*log(5.4000000000000004e+00*1.0/u[6])+2.0000000000000001e-04*u[0]+2.0000000000000000e-02))+1.0000000000000001e-01*1.0/( exp( -1.6028256395817388e+00*log(5.4000000000000004e+00*1.0/u[6])+5.9999999999999998e-02*u[0]-1.2000000000000000e+01)+1.0))/( exp( -1.6028256395817388e+00*log(5.4000000000000004e+00*1.0/u[6])+5.9999999999999998e-02*u[0]-1.2000000000000000e+01)+1.0)+3.4364336248541711e-05*u[11]*( 2.6713760659695648e+01*log(5.4000000000000004e+00*1.0/u[6])-u[0])*u[10]+5.3729341037579627e-04*1.0/( u[5]+4.0000000000000000e+01)*u[5]/( 3.5299999999999998e-02*exp(-3.7433890822745473e-02*u[0])+1.2450000000000000e-01*exp(-3.7433890822745472e-03*u[0])+1.0)+-4.5819114998055616e-05*(u[14]*u[14])*( -2.6713760659695648e+01*log(9.5999999999999996e+00*1.0/( 2.9999999999999999e-02*u[5]+u[6]))+u[0]);
 
77
            y[7] = ( exp(-u[0]/5.0-12.0)+1.0)*( pow( exp( -1.1074197120708749e-01*u[0]-6.2967884828349945e+00)+1.0,-2.0)-u[7])/( 1.0000000000000001e-01*1.0/( exp( u[0]/5.0+7.0)+1.0)+1.0000000000000001e-01*1.0/( exp( u[0]/200.0-(1.0/4.0))+1.0));
 
78
            y[8] = -( u[8]-pow( exp( 1.3458950201884254e-01*u[0]+9.6298788694481825e+00)+1.0,-2.0))*( 5.7000000000000002e-02*1.0/( exp( u[0]+40.0)+1.0)*exp( -1.4705882352941177e-01*u[0]-1.1764705882352942e+01)+-7.7000000000000002e-01*1.0/( 1.3000000000000000e-01*exp( -9.0090090090090100e-02*u[0]-9.6036036036036043e-01)+1.3000000000000000e-01)*( 1.0/( exp( u[0]+40.0)+1.0)-1.0)+1.0/( exp( u[0]+40.0)+1.0)*( 2.7000000000000002e+00*exp(7.9000000000000001e-02*u[0])+3.1000000000000000e+05*exp(3.4849999999999998e-01*u[0])));
 
79
            y[9] = ( pow( exp( 1.3458950201884254e-01*u[0]+9.6298788694481825e+00)+1.0,-2.0)-u[9])*( ( u[0]+3.7780000000000001e+01)*( -6.9480000000000002e-06*exp(-4.3909999999999998e-02*u[0])+-2.5428000000000000e+04*exp(2.4440000000000001e-01*u[0]))/( exp( u[0]+40.0)+1.0)/( exp( 3.1100000000000000e-01*u[0]+2.4640530000000002e+01)+1.0)+2.4240000000000001e-02*1.0/( exp( -1.3780000000000001e-01*u[0]-5.5312920000000005e+00)+1.0)/( exp( u[0]+40.0)+1.0)*exp(-1.0520000000000000e-02*u[0])+-5.9999999999999998e-01*( 1.0/( exp( u[0]+40.0)+1.0)-1.0)*exp(5.7000000000000002e-02*u[0])/( exp( -1.0000000000000001e-01*u[0]-3.2000000000000002e+00)+1.0));
 
80
            y[10] = 1.0/( 9.5000000000000000e+00*exp(-(1.0/1800.0)*pow( u[0]+40.0,2.0))+8.0000000000000004e-01)*( 1.0/( exp(-u[0]/6.0+(10.0/3.0))+1.0)-u[10]);
 
81
            y[11] = 1.0/( 5.0*1.0/( exp( u[0]/5.0-4.0)+1.0)+85.0*exp(-(1.0/320.0)*pow( u[0]+45.0,2.0))+3.0)*( 1.0/( exp( u[0]/5.0+4.0)+1.0)-u[11]);
 
82
            y[12] = -(1.0/2700.0)*( u[12]-1.0/( exp(-u[0]/7.0-(26.0/7.0))+1.0))*( exp(-u[0]/10.0-(9.0/2.0))+1.0)*( exp( 8.6956521739130432e-02*u[0]+2.6086956521739131e+00)+1.0);
 
83
            y[13] = -2.9761904761904762e-01*( u[13]-1.0/( exp( u[0]/24.0+(11.0/3.0))+1.0))*( exp(-u[0]/20.0-3.0)+1.0)*( exp( u[0]/20.0-3.0)+1.0);
 
84
            y[14] = (1.0/40.0)*( 1.0/( exp(-u[0]/14.0-(5.0/14.0))+1.0)-u[14])/( 35.0*1.0/( exp( u[0]/15.0-(7.0/3.0))+1.0)*pow( exp(-u[0]/6.0+(5.0/6.0))+1.0,-(1.0/2.0))+2.0);
 
85
            y[15] = -1.0/( 1.3999999999999999e+00*( 1.3999999999999999e+00*1.0/( exp(-u[0]/13.0-(35.0/13.0))+1.0)+2.5000000000000000e-01)/( exp( u[0]/5.0+1.0)+1.0)+1.0/( exp(-u[0]/20.0+(5.0/2.0))+1.0))*( u[15]-1.0/( exp( -1.3333333333333333e-01*u[0]-1.0666666666666667e+00)+1.0));
 
86
            y[16] = -1.0/( 1.1025000000000000e+03*exp(-(1.0/225.0)*pow( u[0]+27.0,2.0))+200.0*1.0/( exp(-u[0]/10.0+(13.0/10.0))+1.0)+180.0*1.0/( exp( u[0]/10.0+3.0)+1.0)+20.0)*( u[16]-1.0/( exp( u[0]/7.0+(20.0/7.0))+1.0));
 
87
            y[17] = ( 6.7000000000000004e-01*1.0/( exp( u[0]/7.0+5.0)+1.0)-u[17]+3.3000000000000002e-01)/( 562.0*exp(-(1.0/240.0)*pow( u[0]+27.0,2.0))+31.0*1.0/( exp(-u[0]/10.0+(5.0/2.0))+1.0)+80.0*1.0/( exp( u[0]/10.0+3.0)+1.0));
 
88
            y[18] = -(1.0/2.0)*( -5.9999999999999998e-01*1.0/( 3.9999999999999994e+02*(u[4]*u[4])+1.0)+u[18]-4.0000000000000002e-01)/( 40.0*1.0/( 3.9999999999999994e+02*(u[4]*u[4])+1.0)+1.0);
 
89
        }
 
90
        
 
91
        virtual void J(const double* x, double* y, const double* u, double t) {
 
92
            y[0] = -x[10]*( u[0]-R/F*T*log(1.0/u[6]*Ko))*Gto/Cm*u[11]+x[9]*( R*log(1.0/u[5]*Nao)/F*T-u[0])*(u[7]*u[7]*u[7])/Cm*GNa*u[8]+u[9]*( R*log(1.0/u[5]*Nao)/F*T-u[0])*(u[7]*u[7]*u[7])/Cm*GNa*x[8]+-4.0000000000000000e+00*u[15]/R*x[18]*( u[0]-1.5000000000000000e+01)*u[16]*(F*F)*GCaL*u[17]/T/Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)+-4.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*u[16]*(F*F)*GCaL*u[17]/T/Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*x[15]*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-1.0/Cm*x[5]*( 3.0000000000000000e+00*1.0/( KmCa+Cao)*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)*exp(1.0/R*u[0]*F*gamma/T)/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))*Cao*(u[5]*u[5])+PNaK/( 3.5299999999999998e-02*exp(-1.0/R*u[0]*F/T)+1.2450000000000000e-01*exp(-1.0000000000000001e-01*1.0/R*u[0]*F/T)+1.0000000000000000e+00)/( KmNa+u[5])*Ko/( KmK+Ko)+R*u[9]/F*(u[7]*u[7]*u[7])*T*GNa*u[8]/u[5]+R/F*GbNa*T/u[5]-PNaK/( 3.5299999999999998e-02*exp(-1.0/R*u[0]*F/T)+1.2450000000000000e-01*exp(-1.0000000000000001e-01*1.0/R*u[0]*F/T)+1.0000000000000000e+00)/pow( KmNa+u[5],2.0)*Ko/( KmK+Ko)*u[5]+R*pKNa/F/( u[6]+pKNa*u[5])*T*(u[14]*u[14])*GKs)-u[15]/R*( u[0]-1.5000000000000000e+01)*u[16]*(F*F)*x[4]*GCaL*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[17]/T/Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*u[18]+-4.3033148291193518e-01*( u[0]-R/F*T*log(1.0/u[6]*Ko))*u[13]*sqrt(Ko)/Cm*GKr*x[12]+-5.0000000000000000e-01*x[6]/Cm*( 2.0000000000000000e+00*R/u[6]/F*Gto*u[10]*T*u[11]+8.6066296582387036e-01*R*u[12]/u[6]/F*u[13]*sqrt(Ko)*T*GKr+2.0000000000000000e+00*R/F/( u[6]+pKNa*u[5])*T*(u[14]*u[14])*GKs+2.0000000000000000e+00*R/u[6]/F*T*GpK/( exp( -1.6722408026755853e-01*u[0]+4.1806020066889626e+00)+1.0000000000000000e+00)+-5.1639777949432225e-03*R/u[6]*GK1/F*( u[0]-R/F*T*log(1.0/u[6]*Ko))*exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)*sqrt(Ko)/pow( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00,2.0)*T/( 1.0/( exp( -5.0000000000000000e-01*u[0]+5.0000000000000000e-01*R/F*T*log(1.0/u[6]*Ko))+1.0000000000000000e+00)*( 3.0000000000000000e+00*exp( 2.0000000000000001e-04*u[0]+-2.0000000000000001e-04*R/F*T*log(1.0/u[6]*Ko)+2.0000000000000000e-02)+exp( 1.0000000000000001e-01*u[0]+-1.0000000000000001e-01*R/F*T*log(1.0/u[6]*Ko)-1.0000000000000000e+00))+1.0000000000000001e-01*1.0/( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00))+8.6066296582387042e-02*R/u[6]*GK1/F*sqrt(Ko)/( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00)*T/( 1.0/( exp( -5.0000000000000000e-01*u[0]+5.0000000000000000e-01*R/F*T*log(1.0/u[6]*Ko))+1.0000000000000000e+00)*( 3.0000000000000000e+00*exp( 2.0000000000000001e-04*u[0]+-2.0000000000000001e-04*R/F*T*log(1.0/u[6]*Ko)+2.0000000000000000e-02)+exp( 1.0000000000000001e-01*u[0]+-1.0000000000000001e-01*R/F*T*log(1.0/u[6]*Ko)-1.0000000000000000e+00))+1.0000000000000001e-01*1.0/( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00))+-8.6066296582387042e-02*GK1*( u[0]-R/F*T*log(1.0/u[6]*Ko))*( 5.0000000000000000e-01*R*exp( -5.0000000000000000e-01*u[0]+5.0000000000000000e-01*R/F*T*log(1.0/u[6]*Ko))/u[6]/pow( exp( -5.0000000000000000e-01*u[0]+5.0000000000000000e-01*R/F*T*log(1.0/u[6]*Ko))+1.0000000000000000e+00,2.0)*( 3.0000000000000000e+00*exp( 2.0000000000000001e-04*u[0]+-2.0000000000000001e-04*R/F*T*log(1.0/u[6]*Ko)+2.0000000000000000e-02)+exp( 1.0000000000000001e-01*u[0]+-1.0000000000000001e-01*R/F*T*log(1.0/u[6]*Ko)-1.0000000000000000e+00))/F*T+-6.0000000000000001e-03*R/u[6]/F*exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)/pow( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00,2.0)*T+1.0/( exp( -5.0000000000000000e-01*u[0]+5.0000000000000000e-01*R/F*T*log(1.0/u[6]*Ko))+1.0000000000000000e+00)*( 6.0000000000000006e-04*R*exp( 2.0000000000000001e-04*u[0]+-2.0000000000000001e-04*R/F*T*log(1.0/u[6]*Ko)+2.0000000000000000e-02)/u[6]/F*T+1.0000000000000001e-01*R/u[6]/F*T*exp( 1.0000000000000001e-01*u[0]+-1.0000000000000001e-01*R/F*T*log(1.0/u[6]*Ko)-1.0000000000000000e+00)))*sqrt(Ko)/( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00)/pow( 1.0/( exp( -5.0000000000000000e-01*u[0]+5.0000000000000000e-01*R/F*T*log(1.0/u[6]*Ko))+1.0000000000000000e+00)*( 3.0000000000000000e+00*exp( 2.0000000000000001e-04*u[0]+-2.0000000000000001e-04*R/F*T*log(1.0/u[6]*Ko)+2.0000000000000000e-02)+exp( 1.0000000000000001e-01*u[0]+-1.0000000000000001e-01*R/F*T*log(1.0/u[6]*Ko)-1.0000000000000000e+00))+1.0000000000000001e-01*1.0/( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00),2.0))+-2.0000000000000000e+00*x[14]/Cm*u[14]*GKs*( u[0]-R*log(1.0/( u[6]+pKNa*u[5])*( pKNa*Nao+Ko))/F*T)+3.0000000000000000e+00*u[9]*x[7]*( R*log(1.0/u[5]*Nao)/F*T-u[0])*(u[7]*u[7])/Cm*GNa*u[8]-( u[0]-R/F*T*log(1.0/u[6]*Ko))*Gto*u[10]*x[11]/Cm+-5.0000000000000000e-01*x[0]*( 2.0000000000000000e+00*GbCa+-8.6066296582387042e-02*( 1.0/( exp( -5.0000000000000000e-01*u[0]+5.0000000000000000e-01*R/F*T*log(1.0/u[6]*Ko))+1.0000000000000000e+00)*( 6.0000000000000006e-04*exp( 2.0000000000000001e-04*u[0]+-2.0000000000000001e-04*R/F*T*log(1.0/u[6]*Ko)+2.0000000000000000e-02)+1.0000000000000001e-01*exp( 1.0000000000000001e-01*u[0]+-1.0000000000000001e-01*R/F*T*log(1.0/u[6]*Ko)-1.0000000000000000e+00))+-6.0000000000000001e-03*exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)/pow( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00,2.0)+5.0000000000000000e-01*exp( -5.0000000000000000e-01*u[0]+5.0000000000000000e-01*R/F*T*log(1.0/u[6]*Ko))/pow( exp( -5.0000000000000000e-01*u[0]+5.0000000000000000e-01*R/F*T*log(1.0/u[6]*Ko))+1.0000000000000000e+00,2.0)*( 3.0000000000000000e+00*exp( 2.0000000000000001e-04*u[0]+-2.0000000000000001e-04*R/F*T*log(1.0/u[6]*Ko)+2.0000000000000000e-02)+exp( 1.0000000000000001e-01*u[0]+-1.0000000000000001e-01*R/F*T*log(1.0/u[6]*Ko)-1.0000000000000000e+00)))*GK1*( u[0]-R/F*T*log(1.0/u[6]*Ko))*sqrt(Ko)/( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00)/pow( 1.0/( exp( -5.0000000000000000e-01*u[0]+5.0000000000000000e-01*R/F*T*log(1.0/u[6]*Ko))+1.0000000000000000e+00)*( 3.0000000000000000e+00*exp( 2.0000000000000001e-04*u[0]+-2.0000000000000001e-04*R/F*T*log(1.0/u[6]*Ko)+2.0000000000000000e-02)+exp( 1.0000000000000001e-01*u[0]+-1.0000000000000001e-01*R/F*T*log(1.0/u[6]*Ko)-1.0000000000000000e+00))+1.0000000000000001e-01*1.0/( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00),2.0)+8.6066296582387042e-02*GK1*sqrt(Ko)/( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00)/( 1.0/( exp( -5.0000000000000000e-01*u[0]+5.0000000000000000e-01*R/F*T*log(1.0/u[6]*Ko))+1.0000000000000000e+00)*( 3.0000000000000000e+00*exp( 2.0000000000000001e-04*u[0]+-2.0000000000000001e-04*R/F*T*log(1.0/u[6]*Ko)+2.0000000000000000e-02)+exp( 1.0000000000000001e-01*u[0]+-1.0000000000000001e-01*R/F*T*log(1.0/u[6]*Ko)-1.0000000000000000e+00))+1.0000000000000001e-01*1.0/( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00))+2.0000000000000000e+00*u[9]*(u[7]*u[7]*u[7])*GNa*u[8]+8.6066296582387036e-01*u[12]*u[13]*sqrt(Ko)*GKr+4.0000000000000000e+00*u[15]/(R*R)*( u[0]-1.5000000000000000e+01)*u[16]*(F*F*F)*GCaL*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]*u[17]/(T*T)/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*u[18]+-2.0000000000000000e+00*1.0/R/( KmCa+Cao)*kNaCa/pow( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00,2.0)*exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*F*( gamma-1.0000000000000000e+00)*( exp(1.0/R*u[0]*F*gamma/T)*Cao*(u[5]*u[5]*u[5])-exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*alpha*u[2]*(Nao*Nao*Nao))/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))/T*ksat+2.0000000000000000e+00*GbNa+-1.6000000000000000e+01*u[15]/(R*R)*( u[0]-1.5000000000000000e+01)*u[16]*(F*F*F)*GCaL*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[17]/(T*T)/pow( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00,2.0)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]+2.0000000000000000e+00*Gto*u[10]*u[11]+2.0000000000000000e+00*1.0/( KmCa+Cao)*( 1.0/R*F*gamma*exp(1.0/R*u[0]*F*gamma/T)/T*Cao*(u[5]*u[5]*u[5])-1.0/R*exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*F*( gamma-1.0000000000000000e+00)/T*alpha*u[2]*(Nao*Nao*Nao))*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))+8.0000000000000000e+00*u[15]/R*u[16]*(F*F)*GCaL*u[17]/T/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]+2.0000000000000000e+00*GpK/( exp( -1.6722408026755853e-01*u[0]+4.1806020066889626e+00)+1.0000000000000000e+00)+-2.0000000000000000e+00*PNaK/pow( 3.5299999999999998e-02*exp(-1.0/R*u[0]*F/T)+1.2450000000000000e-01*exp(-1.0000000000000001e-01*1.0/R*u[0]*F/T)+1.0000000000000000e+00,2.0)/( KmNa+u[5])*( -1.2450000000000001e-02*1.0/R*F/T*exp(-1.0000000000000001e-01*1.0/R*u[0]*F/T)+-3.5299999999999998e-02*1.0/R*exp(-1.0/R*u[0]*F/T)*F/T)*Ko/( KmK+Ko)*u[5]+3.3444816053511706e-01*( u[0]-R/F*T*log(1.0/u[6]*Ko))*GpK*exp( -1.6722408026755853e-01*u[0]+4.1806020066889626e+00)/pow( exp( -1.6722408026755853e-01*u[0]+4.1806020066889626e+00)+1.0000000000000000e+00,2.0)+2.0000000000000000e+00*(u[14]*u[14])*GKs+-5.1639777949432225e-03*GK1*( u[0]-R/F*T*log(1.0/u[6]*Ko))*exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)*sqrt(Ko)/pow( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00,2.0)/( 1.0/( exp( -5.0000000000000000e-01*u[0]+5.0000000000000000e-01*R/F*T*log(1.0/u[6]*Ko))+1.0000000000000000e+00)*( 3.0000000000000000e+00*exp( 2.0000000000000001e-04*u[0]+-2.0000000000000001e-04*R/F*T*log(1.0/u[6]*Ko)+2.0000000000000000e-02)+exp( 1.0000000000000001e-01*u[0]+-1.0000000000000001e-01*R/F*T*log(1.0/u[6]*Ko)-1.0000000000000000e+00))+1.0000000000000001e-01*1.0/( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00)))/Cm+-4.3033148291193518e-01*u[12]*x[13]*( u[0]-R/F*T*log(1.0/u[6]*Ko))*sqrt(Ko)/Cm*GKr+-4.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)*u[16]*(F*F)*GCaL/T/Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*x[17]*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]+-4.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)*x[16]*(F*F)*GCaL*u[17]/T/Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]+-5.0000000000000000e-01*1.0/Cm*x[2]*( 2.0000000000000000e+00*GpCa/( KpCa+u[2])+-2.0000000000000000e+00*1.0/( KmCa+Cao)*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)*exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))*alpha*(Nao*Nao*Nao)+-2.0000000000000000e+00*GpCa/pow( KpCa+u[2],2.0)*u[2]+R*GbCa/F*T/u[2]);
 
93
            y[1] = -x[1]*( k4+u[4]*( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))*km2)-u[1]*x[4]*( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))*km2+2.0000000000000000e+00*1.0/(u[3]*u[3]*u[3])*( maxSR-minSR)*x[3]*u[1]*u[4]/pow( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00,2.0)*(EC*EC)*km2;
 
94
            y[2] =  -6.2500000000000000e-02*x[0]/dt*( 4.0000000000000000e+00*1.0/F*dt/Vc*Cm*( GbCa+2.0000000000000000e+00*1.0/R/( KmCa+Cao)*kNaCa/pow( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00,2.0)*exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*F*( gamma-1.0000000000000000e+00)*( exp(1.0/R*u[0]*F*gamma/T)*Cao*(u[5]*u[5]*u[5])-exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*alpha*u[2]*(Nao*Nao*Nao))/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))/T*ksat+-2.0000000000000000e+00*1.0/( KmCa+Cao)*( 1.0/R*F*gamma*exp(1.0/R*u[0]*F*gamma/T)/T*Cao*(u[5]*u[5]*u[5])-1.0/R*exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*F*( gamma-1.0000000000000000e+00)/T*alpha*u[2]*(Nao*Nao*Nao))*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao)))+( 1.0/F*dt/Vc*Cm*( dt*( 4.0000000000000000e+00*Vxfer*( u[4]-u[2])+( 4.0000000000000000e+00*1.0/( KmCa+Cao)*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)*( exp(1.0/R*u[0]*F*gamma/T)*Cao*(u[5]*u[5]*u[5])-exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*alpha*u[2]*(Nao*Nao*Nao))/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))+-2.0000000000000000e+00*GpCa/( KpCa+u[2])*u[2]-GbCa*( 2.0000000000000000e+00*u[0]-R/F*log(1.0/u[2]*Cao)*T))/F/Vc*Cm+-4.0000000000000000e+00*( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2]))/Vc*Vsr)+-4.0000000000000000e+00*Bufc+4.0000000000000000e+00*u[2]+4.0000000000000000e+00*1.0/( u[2]+Kbufc)*Bufc*u[2]+-4.0000000000000000e+00*Kbufc)*( GbCa+2.0000000000000000e+00*1.0/R/( KmCa+Cao)*kNaCa/pow( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00,2.0)*exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*F*( gamma-1.0000000000000000e+00)*( exp(1.0/R*u[0]*F*gamma/T)*Cao*(u[5]*u[5]*u[5])-exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*alpha*u[2]*(Nao*Nao*Nao))/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))/T*ksat+-2.0000000000000000e+00*1.0/( KmCa+Cao)*( 1.0/R*F*gamma*exp(1.0/R*u[0]*F*gamma/T)/T*Cao*(u[5]*u[5]*u[5])-1.0/R*exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*F*( gamma-1.0000000000000000e+00)/T*alpha*u[2]*(Nao*Nao*Nao))*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao)))+8.0000000000000000e+00*1.0/F*dt/Vc*Cm*( GbCa+2.0000000000000000e+00*1.0/R/( KmCa+Cao)*kNaCa/pow( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00,2.0)*exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*F*( gamma-1.0000000000000000e+00)*( exp(1.0/R*u[0]*F*gamma/T)*Cao*(u[5]*u[5]*u[5])-exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*alpha*u[2]*(Nao*Nao*Nao))/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))/T*ksat+-2.0000000000000000e+00*1.0/( KmCa+Cao)*( 1.0/R*F*gamma*exp(1.0/R*u[0]*F*gamma/T)/T*Cao*(u[5]*u[5]*u[5])-1.0/R*exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*F*( gamma-1.0000000000000000e+00)/T*alpha*u[2]*(Nao*Nao*Nao))*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao)))*Kbufc)*pow( 6.2500000000000000e-02*pow( dt*( 4.0000000000000000e+00*Vxfer*( u[4]-u[2])+( 4.0000000000000000e+00*1.0/( KmCa+Cao)*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)*( exp(1.0/R*u[0]*F*gamma/T)*Cao*(u[5]*u[5]*u[5])-exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*alpha*u[2]*(Nao*Nao*Nao))/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))+-2.0000000000000000e+00*GpCa/( KpCa+u[2])*u[2]-GbCa*( 2.0000000000000000e+00*u[0]-R/F*log(1.0/u[2]*Cao)*T))/F/Vc*Cm+-4.0000000000000000e+00*( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2]))/Vc*Vsr)+-4.0000000000000000e+00*Bufc+4.0000000000000000e+00*u[2]+4.0000000000000000e+00*1.0/( u[2]+Kbufc)*Bufc*u[2]+-4.0000000000000000e+00*Kbufc,2.0)+( dt*( 4.0000000000000000e+00*Vxfer*( u[4]-u[2])+( 4.0000000000000000e+00*1.0/( KmCa+Cao)*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)*( exp(1.0/R*u[0]*F*gamma/T)*Cao*(u[5]*u[5]*u[5])-exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*alpha*u[2]*(Nao*Nao*Nao))/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))+-2.0000000000000000e+00*GpCa/( KpCa+u[2])*u[2]-GbCa*( 2.0000000000000000e+00*u[0]-R/F*log(1.0/u[2]*Cao)*T))/F/Vc*Cm+-4.0000000000000000e+00*( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2]))/Vc*Vsr)+4.0000000000000000e+00*u[2]+4.0000000000000000e+00*1.0/( u[2]+Kbufc)*Bufc*u[2])*Kbufc,-(1.0/2.0)))+3.7500000000000000e-01*( 4.0000000000000000e+00*1.0/( KmCa+Cao)*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)/F*dt*exp(1.0/R*u[0]*F*gamma/T)/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))/Vc*Cm*Cao*(u[5]*u[5])+( 8.0000000000000000e+00*1.0/( KmCa+Cao)*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)/F*dt*exp(1.0/R*u[0]*F*gamma/T)/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))/Vc*Cm*Cao*Kbufc*(u[5]*u[5])+1.0/( KmCa+Cao)*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)/F*dt*exp(1.0/R*u[0]*F*gamma/T)/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))/Vc*Cm*( dt*( 4.0000000000000000e+00*Vxfer*( u[4]-u[2])+( 4.0000000000000000e+00*1.0/( KmCa+Cao)*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)*( exp(1.0/R*u[0]*F*gamma/T)*Cao*(u[5]*u[5]*u[5])-exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*alpha*u[2]*(Nao*Nao*Nao))/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))+-2.0000000000000000e+00*GpCa/( KpCa+u[2])*u[2]-GbCa*( 2.0000000000000000e+00*u[0]-R/F*log(1.0/u[2]*Cao)*T))/F/Vc*Cm+-4.0000000000000000e+00*( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2]))/Vc*Vsr)+-4.0000000000000000e+00*Bufc+4.0000000000000000e+00*u[2]+4.0000000000000000e+00*1.0/( u[2]+Kbufc)*Bufc*u[2]+-4.0000000000000000e+00*Kbufc)*Cao*(u[5]*u[5]))*pow( 6.2500000000000000e-02*pow( dt*( 4.0000000000000000e+00*Vxfer*( u[4]-u[2])+( 4.0000000000000000e+00*1.0/( KmCa+Cao)*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)*( exp(1.0/R*u[0]*F*gamma/T)*Cao*(u[5]*u[5]*u[5])-exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*alpha*u[2]*(Nao*Nao*Nao))/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))+-2.0000000000000000e+00*GpCa/( KpCa+u[2])*u[2]-GbCa*( 2.0000000000000000e+00*u[0]-R/F*log(1.0/u[2]*Cao)*T))/F/Vc*Cm+-4.0000000000000000e+00*( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2]))/Vc*Vsr)+-4.0000000000000000e+00*Bufc+4.0000000000000000e+00*u[2]+4.0000000000000000e+00*1.0/( u[2]+Kbufc)*Bufc*u[2]+-4.0000000000000000e+00*Kbufc,2.0)+( dt*( 4.0000000000000000e+00*Vxfer*( u[4]-u[2])+( 4.0000000000000000e+00*1.0/( KmCa+Cao)*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)*( exp(1.0/R*u[0]*F*gamma/T)*Cao*(u[5]*u[5]*u[5])-exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*alpha*u[2]*(Nao*Nao*Nao))/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))+-2.0000000000000000e+00*GpCa/( KpCa+u[2])*u[2]-GbCa*( 2.0000000000000000e+00*u[0]-R/F*log(1.0/u[2]*Cao)*T))/F/Vc*Cm+-4.0000000000000000e+00*( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2]))/Vc*Vsr)+4.0000000000000000e+00*u[2]+4.0000000000000000e+00*1.0/( u[2]+Kbufc)*Bufc*u[2])*Kbufc,-(1.0/2.0)))/dt*x[5]+-3.1250000000000000e-02*( 1.6000000000000000e+01*1.0/pow( u[2]+Kbufc,2.0)*Bufc*u[2]+-1.6000000000000000e+01*1.0/( u[2]+Kbufc)*Bufc+4.0000000000000000e+00*dt*( 4.0000000000000000e+00*Vxfer+4.0000000000000000e+00*( Vleak+2.0000000000000000e+00*1.0/pow( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00,2.0)*Vmaxup*(Kup*Kup)/(u[2]*u[2]*u[2]))/Vc*Vsr+1.0/F/Vc*Cm*( 2.0000000000000000e+00*GpCa/( KpCa+u[2])+4.0000000000000000e+00*1.0/( KmCa+Cao)*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)*exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))*alpha*(Nao*Nao*Nao)+-2.0000000000000000e+00*GpCa/pow( KpCa+u[2],2.0)*u[2]+R*GbCa/F*T/u[2]))+pow( 6.2500000000000000e-02*pow( dt*( 4.0000000000000000e+00*Vxfer*( u[4]-u[2])+( 4.0000000000000000e+00*1.0/( KmCa+Cao)*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)*( exp(1.0/R*u[0]*F*gamma/T)*Cao*(u[5]*u[5]*u[5])-exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*alpha*u[2]*(Nao*Nao*Nao))/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))+-2.0000000000000000e+00*GpCa/( KpCa+u[2])*u[2]-GbCa*( 2.0000000000000000e+00*u[0]-R/F*log(1.0/u[2]*Cao)*T))/F/Vc*Cm+-4.0000000000000000e+00*( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2]))/Vc*Vsr)+-4.0000000000000000e+00*Bufc+4.0000000000000000e+00*u[2]+4.0000000000000000e+00*1.0/( u[2]+Kbufc)*Bufc*u[2]+-4.0000000000000000e+00*Kbufc,2.0)+( dt*( 4.0000000000000000e+00*Vxfer*( u[4]-u[2])+( 4.0000000000000000e+00*1.0/( KmCa+Cao)*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)*( exp(1.0/R*u[0]*F*gamma/T)*Cao*(u[5]*u[5]*u[5])-exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*alpha*u[2]*(Nao*Nao*Nao))/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))+-2.0000000000000000e+00*GpCa/( KpCa+u[2])*u[2]-GbCa*( 2.0000000000000000e+00*u[0]-R/F*log(1.0/u[2]*Cao)*T))/F/Vc*Cm+-4.0000000000000000e+00*( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2]))/Vc*Vsr)+4.0000000000000000e+00*u[2]+4.0000000000000000e+00*1.0/( u[2]+Kbufc)*Bufc*u[2])*Kbufc,-(1.0/2.0))*( 8.0000000000000000e+00*( 4.0000000000000000e+00*1.0/pow( u[2]+Kbufc,2.0)*Bufc*u[2]+-4.0000000000000000e+00*1.0/( u[2]+Kbufc)*Bufc+dt*( 4.0000000000000000e+00*Vxfer+4.0000000000000000e+00*( Vleak+2.0000000000000000e+00*1.0/pow( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00,2.0)*Vmaxup*(Kup*Kup)/(u[2]*u[2]*u[2]))/Vc*Vsr+1.0/F/Vc*Cm*( 2.0000000000000000e+00*GpCa/( KpCa+u[2])+4.0000000000000000e+00*1.0/( KmCa+Cao)*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)*exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))*alpha*(Nao*Nao*Nao)+-2.0000000000000000e+00*GpCa/pow( KpCa+u[2],2.0)*u[2]+R*GbCa/F*T/u[2]))-4.0000000000000000e+00)*Kbufc+( 4.0000000000000000e+00*1.0/pow( u[2]+Kbufc,2.0)*Bufc*u[2]+-4.0000000000000000e+00*1.0/( u[2]+Kbufc)*Bufc+dt*( 4.0000000000000000e+00*Vxfer+4.0000000000000000e+00*( Vleak+2.0000000000000000e+00*1.0/pow( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00,2.0)*Vmaxup*(Kup*Kup)/(u[2]*u[2]*u[2]))/Vc*Vsr+1.0/F/Vc*Cm*( 2.0000000000000000e+00*GpCa/( KpCa+u[2])+4.0000000000000000e+00*1.0/( KmCa+Cao)*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)*exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))*alpha*(Nao*Nao*Nao)+-2.0000000000000000e+00*GpCa/pow( KpCa+u[2],2.0)*u[2]+R*GbCa/F*T/u[2]))-4.0000000000000000e+00)*( dt*( 4.0000000000000000e+00*Vxfer*( u[4]-u[2])+( 4.0000000000000000e+00*1.0/( KmCa+Cao)*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)*( exp(1.0/R*u[0]*F*gamma/T)*Cao*(u[5]*u[5]*u[5])-exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*alpha*u[2]*(Nao*Nao*Nao))/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))+-2.0000000000000000e+00*GpCa/( KpCa+u[2])*u[2]-GbCa*( 2.0000000000000000e+00*u[0]-R/F*log(1.0/u[2]*Cao)*T))/F/Vc*Cm+-4.0000000000000000e+00*( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2]))/Vc*Vsr)+-4.0000000000000000e+00*Bufc+4.0000000000000000e+00*u[2]+4.0000000000000000e+00*1.0/( u[2]+Kbufc)*Bufc*u[2]+-4.0000000000000000e+00*Kbufc))+1.6000000000000000e+01)/dt*x[2]+1.2500000000000000e-01*( ( Vleak*dt/Vc*( dt*( 4.0000000000000000e+00*Vxfer*( u[4]-u[2])+( 4.0000000000000000e+00*1.0/( KmCa+Cao)*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)*( exp(1.0/R*u[0]*F*gamma/T)*Cao*(u[5]*u[5]*u[5])-exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*alpha*u[2]*(Nao*Nao*Nao))/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))+-2.0000000000000000e+00*GpCa/( KpCa+u[2])*u[2]-GbCa*( 2.0000000000000000e+00*u[0]-R/F*log(1.0/u[2]*Cao)*T))/F/Vc*Cm+-4.0000000000000000e+00*( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2]))/Vc*Vsr)+-4.0000000000000000e+00*Bufc+4.0000000000000000e+00*u[2]+4.0000000000000000e+00*1.0/( u[2]+Kbufc)*Bufc*u[2]+-4.0000000000000000e+00*Kbufc)*Vsr+8.0000000000000000e+00*Vleak*dt/Vc*Vsr*Kbufc)*pow( 6.2500000000000000e-02*pow( dt*( 4.0000000000000000e+00*Vxfer*( u[4]-u[2])+( 4.0000000000000000e+00*1.0/( KmCa+Cao)*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)*( exp(1.0/R*u[0]*F*gamma/T)*Cao*(u[5]*u[5]*u[5])-exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*alpha*u[2]*(Nao*Nao*Nao))/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))+-2.0000000000000000e+00*GpCa/( KpCa+u[2])*u[2]-GbCa*( 2.0000000000000000e+00*u[0]-R/F*log(1.0/u[2]*Cao)*T))/F/Vc*Cm+-4.0000000000000000e+00*( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2]))/Vc*Vsr)+-4.0000000000000000e+00*Bufc+4.0000000000000000e+00*u[2]+4.0000000000000000e+00*1.0/( u[2]+Kbufc)*Bufc*u[2]+-4.0000000000000000e+00*Kbufc,2.0)+( dt*( 4.0000000000000000e+00*Vxfer*( u[4]-u[2])+( 4.0000000000000000e+00*1.0/( KmCa+Cao)*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)*( exp(1.0/R*u[0]*F*gamma/T)*Cao*(u[5]*u[5]*u[5])-exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*alpha*u[2]*(Nao*Nao*Nao))/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))+-2.0000000000000000e+00*GpCa/( KpCa+u[2])*u[2]-GbCa*( 2.0000000000000000e+00*u[0]-R/F*log(1.0/u[2]*Cao)*T))/F/Vc*Cm+-4.0000000000000000e+00*( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2]))/Vc*Vsr)+4.0000000000000000e+00*u[2]+4.0000000000000000e+00*1.0/( u[2]+Kbufc)*Bufc*u[2])*Kbufc,-(1.0/2.0))+4.0000000000000000e+00*Vleak*dt/Vc*Vsr)*x[3]/dt+1.2500000000000000e-01*( ( 8.0000000000000000e+00*Vxfer*dt*Kbufc+Vxfer*dt*( dt*( 4.0000000000000000e+00*Vxfer*( u[4]-u[2])+( 4.0000000000000000e+00*1.0/( KmCa+Cao)*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)*( exp(1.0/R*u[0]*F*gamma/T)*Cao*(u[5]*u[5]*u[5])-exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*alpha*u[2]*(Nao*Nao*Nao))/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))+-2.0000000000000000e+00*GpCa/( KpCa+u[2])*u[2]-GbCa*( 2.0000000000000000e+00*u[0]-R/F*log(1.0/u[2]*Cao)*T))/F/Vc*Cm+-4.0000000000000000e+00*( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2]))/Vc*Vsr)+-4.0000000000000000e+00*Bufc+4.0000000000000000e+00*u[2]+4.0000000000000000e+00*1.0/( u[2]+Kbufc)*Bufc*u[2]+-4.0000000000000000e+00*Kbufc))*pow( 6.2500000000000000e-02*pow( dt*( 4.0000000000000000e+00*Vxfer*( u[4]-u[2])+( 4.0000000000000000e+00*1.0/( KmCa+Cao)*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)*( exp(1.0/R*u[0]*F*gamma/T)*Cao*(u[5]*u[5]*u[5])-exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*alpha*u[2]*(Nao*Nao*Nao))/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))+-2.0000000000000000e+00*GpCa/( KpCa+u[2])*u[2]-GbCa*( 2.0000000000000000e+00*u[0]-R/F*log(1.0/u[2]*Cao)*T))/F/Vc*Cm+-4.0000000000000000e+00*( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2]))/Vc*Vsr)+-4.0000000000000000e+00*Bufc+4.0000000000000000e+00*u[2]+4.0000000000000000e+00*1.0/( u[2]+Kbufc)*Bufc*u[2]+-4.0000000000000000e+00*Kbufc,2.0)+( dt*( 4.0000000000000000e+00*Vxfer*( u[4]-u[2])+( 4.0000000000000000e+00*1.0/( KmCa+Cao)*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)*( exp(1.0/R*u[0]*F*gamma/T)*Cao*(u[5]*u[5]*u[5])-exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*alpha*u[2]*(Nao*Nao*Nao))/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))+-2.0000000000000000e+00*GpCa/( KpCa+u[2])*u[2]-GbCa*( 2.0000000000000000e+00*u[0]-R/F*log(1.0/u[2]*Cao)*T))/F/Vc*Cm+-4.0000000000000000e+00*( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2]))/Vc*Vsr)+4.0000000000000000e+00*u[2]+4.0000000000000000e+00*1.0/( u[2]+Kbufc)*Bufc*u[2])*Kbufc,-(1.0/2.0))+4.0000000000000000e+00*Vxfer*dt)/dt*x[4];
 
95
            y[3] =  -5.0000000000000000e-01*( km1*( u[3]-u[4])*Vrel*dt*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+( km1*( u[3]-u[4])*Vrel*dt*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))*( ( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2])-km1*( u[3]-u[4])*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*dt+u[3]-Kbufsr+u[3]*Bufsr/( u[3]+Kbufsr)-Bufsr)+2.0000000000000000e+00*km1*( u[3]-u[4])*Kbufsr*Vrel*dt*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*pow( 4.0000000000000000e+00*Kbufsr*( ( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2])-km1*( u[3]-u[4])*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*dt+u[3]+u[3]*Bufsr/( u[3]+Kbufsr))+pow( ( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2])-km1*( u[3]-u[4])*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*dt+u[3]-Kbufsr+u[3]*Bufsr/( u[3]+Kbufsr)-Bufsr,2.0),-(1.0/2.0)))/dt*x[1]+5.0000000000000000e-01*x[3]*( Bufsr/( u[3]+Kbufsr)-dt*( 2.0000000000000000e+00*km1/(u[3]*u[3]*u[3])*( maxSR-minSR)*( u[3]-u[4])*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/pow( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00,2.0)/pow( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00),2.0)*(EC*EC)+Vleak+-2.0000000000000000e+00*(km1*km1)/(u[3]*u[3]*u[3])*( maxSR-minSR)*( u[3]-u[4])*Vrel*u[1]*((u[4]*u[4])*(u[4]*u[4]))/pow( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)),2.0)/pow( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00,2.0)/pow( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00),3.0)*(EC*EC)+km1*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))-u[3]*Bufsr/pow( u[3]+Kbufsr,2.0)+( 2.0000000000000000e+00*Kbufsr*( Bufsr/( u[3]+Kbufsr)-dt*( 2.0000000000000000e+00*km1/(u[3]*u[3]*u[3])*( maxSR-minSR)*( u[3]-u[4])*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/pow( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00,2.0)/pow( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00),2.0)*(EC*EC)+Vleak+-2.0000000000000000e+00*(km1*km1)/(u[3]*u[3]*u[3])*( maxSR-minSR)*( u[3]-u[4])*Vrel*u[1]*((u[4]*u[4])*(u[4]*u[4]))/pow( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)),2.0)/pow( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00,2.0)/pow( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00),3.0)*(EC*EC)+km1*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))-u[3]*Bufsr/pow( u[3]+Kbufsr,2.0)+1.0000000000000000e+00)+( Bufsr/( u[3]+Kbufsr)-dt*( 2.0000000000000000e+00*km1/(u[3]*u[3]*u[3])*( maxSR-minSR)*( u[3]-u[4])*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/pow( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00,2.0)/pow( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00),2.0)*(EC*EC)+Vleak+-2.0000000000000000e+00*(km1*km1)/(u[3]*u[3]*u[3])*( maxSR-minSR)*( u[3]-u[4])*Vrel*u[1]*((u[4]*u[4])*(u[4]*u[4]))/pow( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)),2.0)/pow( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00,2.0)/pow( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00),3.0)*(EC*EC)+km1*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))-u[3]*Bufsr/pow( u[3]+Kbufsr,2.0)+1.0000000000000000e+00)*( ( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2])-km1*( u[3]-u[4])*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*dt+u[3]-Kbufsr+u[3]*Bufsr/( u[3]+Kbufsr)-Bufsr))*pow( 4.0000000000000000e+00*Kbufsr*( ( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2])-km1*( u[3]-u[4])*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*dt+u[3]+u[3]*Bufsr/( u[3]+Kbufsr))+pow( ( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2])-km1*( u[3]-u[4])*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*dt+u[3]-Kbufsr+u[3]*Bufsr/( u[3]+Kbufsr)-Bufsr,2.0),-(1.0/2.0))-1.0000000000000000e+00)/dt+5.0000000000000000e-01*( ( Vleak+2.0000000000000000e+00*1.0/pow( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00,2.0)*Vmaxup*(Kup*Kup)/(u[2]*u[2]*u[2]))*dt+( ( Vleak+2.0000000000000000e+00*1.0/pow( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00,2.0)*Vmaxup*(Kup*Kup)/(u[2]*u[2]*u[2]))*dt*( ( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2])-km1*( u[3]-u[4])*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*dt+u[3]-Kbufsr+u[3]*Bufsr/( u[3]+Kbufsr)-Bufsr)+2.0000000000000000e+00*Kbufsr*( Vleak+2.0000000000000000e+00*1.0/pow( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00,2.0)*Vmaxup*(Kup*Kup)/(u[2]*u[2]*u[2]))*dt)*pow( 4.0000000000000000e+00*Kbufsr*( ( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2])-km1*( u[3]-u[4])*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*dt+u[3]+u[3]*Bufsr/( u[3]+Kbufsr))+pow( ( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2])-km1*( u[3]-u[4])*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*dt+u[3]-Kbufsr+u[3]*Bufsr/( u[3]+Kbufsr)-Bufsr,2.0),-(1.0/2.0)))/dt*x[2]+-5.0000000000000000e-01*1.0/dt*x[4]*( dt*( 2.0000000000000000e+00*km1*( u[3]-u[4])*Vrel*u[1]*u[4]/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))-km1*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+-2.0000000000000000e+00*(km1*km1)*( u[3]-u[4])*Vrel*u[1]*(u[4]*u[4]*u[4])/pow( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)),2.0)/pow( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00),2.0))+( 2.0000000000000000e+00*Kbufsr*dt*( 2.0000000000000000e+00*km1*( u[3]-u[4])*Vrel*u[1]*u[4]/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))-km1*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+-2.0000000000000000e+00*(km1*km1)*( u[3]-u[4])*Vrel*u[1]*(u[4]*u[4]*u[4])/pow( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)),2.0)/pow( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00),2.0))+dt*( 2.0000000000000000e+00*km1*( u[3]-u[4])*Vrel*u[1]*u[4]/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))-km1*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+-2.0000000000000000e+00*(km1*km1)*( u[3]-u[4])*Vrel*u[1]*(u[4]*u[4]*u[4])/pow( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)),2.0)/pow( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00),2.0))*( ( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2])-km1*( u[3]-u[4])*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*dt+u[3]-Kbufsr+u[3]*Bufsr/( u[3]+Kbufsr)-Bufsr))*pow( 4.0000000000000000e+00*Kbufsr*( ( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2])-km1*( u[3]-u[4])*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*dt+u[3]+u[3]*Bufsr/( u[3]+Kbufsr))+pow( ( 1.0/( (Kup*Kup)/(u[2]*u[2])+1.0000000000000000e+00)*Vmaxup-Vleak*( u[3]-u[2])-km1*( u[3]-u[4])*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*dt+u[3]-Kbufsr+u[3]*Bufsr/( u[3]+Kbufsr)-Bufsr,2.0),-(1.0/2.0)));
 
96
            y[4] =  5.0000000000000000e-01*1.0/dt*( Vxfer/Vss*dt*Vc+pow( 4.0000000000000000e+00*Kbufss*( Bufss*u[4]/( Kbufss+u[4])-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)+u[4])+pow( Bufss*u[4]/( Kbufss+u[4])-Bufss-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)-Kbufss+u[4],2.0),-(1.0/2.0))*( Vxfer*( Bufss*u[4]/( Kbufss+u[4])-Bufss-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)-Kbufss+u[4])/Vss*dt*Vc+2.0000000000000000e+00*Vxfer/Vss*dt*Vc*Kbufss))*x[2]-x[18]*( u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*dt*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)+( u[15]/R*( u[0]-1.5000000000000000e+01)*( Bufss*u[4]/( Kbufss+u[4])-Bufss-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)-Kbufss+u[4])/Vss*u[16]*F*dt*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)+2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*dt*GCaL*Kbufss*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao))*pow( 4.0000000000000000e+00*Kbufss*( Bufss*u[4]/( Kbufss+u[4])-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)+u[4])+pow( Bufss*u[4]/( Kbufss+u[4])-Bufss-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)-Kbufss+u[4],2.0),-(1.0/2.0)))/dt+5.0000000000000000e-01*( dt*( 4.0000000000000000e+00*u[15]/(R*R)*( u[0]-1.5000000000000000e+01)/Vss*u[16]*(F*F)*GCaL*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[17]/(T*T)*Cm/pow( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00,2.0)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]+-2.0000000000000000e+00*u[15]/R/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-u[15]/(R*R)*( u[0]-1.5000000000000000e+01)/Vss*u[16]*(F*F)*GCaL*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]*u[17]/(T*T)*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*u[18])+pow( 4.0000000000000000e+00*Kbufss*( Bufss*u[4]/( Kbufss+u[4])-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)+u[4])+pow( Bufss*u[4]/( Kbufss+u[4])-Bufss-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)-Kbufss+u[4],2.0),-(1.0/2.0))*( 2.0000000000000000e+00*dt*Kbufss*( 4.0000000000000000e+00*u[15]/(R*R)*( u[0]-1.5000000000000000e+01)/Vss*u[16]*(F*F)*GCaL*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[17]/(T*T)*Cm/pow( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00,2.0)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]+-2.0000000000000000e+00*u[15]/R/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-u[15]/(R*R)*( u[0]-1.5000000000000000e+01)/Vss*u[16]*(F*F)*GCaL*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]*u[17]/(T*T)*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*u[18])+( Bufss*u[4]/( Kbufss+u[4])-Bufss-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)-Kbufss+u[4])*dt*( 4.0000000000000000e+00*u[15]/(R*R)*( u[0]-1.5000000000000000e+01)/Vss*u[16]*(F*F)*GCaL*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[17]/(T*T)*Cm/pow( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00,2.0)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]+-2.0000000000000000e+00*u[15]/R/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-u[15]/(R*R)*( u[0]-1.5000000000000000e+01)/Vss*u[16]*(F*F)*GCaL*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]*u[17]/(T*T)*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*u[18])))*x[0]/dt-( u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*F*dt*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]+pow( 4.0000000000000000e+00*Kbufss*( Bufss*u[4]/( Kbufss+u[4])-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)+u[4])+pow( Bufss*u[4]/( Kbufss+u[4])-Bufss-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)-Kbufss+u[4],2.0),-(1.0/2.0))*( u[15]/R*( u[0]-1.5000000000000000e+01)*( Bufss*u[4]/( Kbufss+u[4])-Bufss-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)-Kbufss+u[4])/Vss*F*dt*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]+2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*F*dt*GCaL*Kbufss*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]))*x[16]/dt-1.0/dt*( 1.0/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*dt*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]+pow( 4.0000000000000000e+00*Kbufss*( Bufss*u[4]/( Kbufss+u[4])-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)+u[4])+pow( Bufss*u[4]/( Kbufss+u[4])-Bufss-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)-Kbufss+u[4],2.0),-(1.0/2.0))*( 1.0/R*( u[0]-1.5000000000000000e+01)*( Bufss*u[4]/( Kbufss+u[4])-Bufss-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)-Kbufss+u[4])/Vss*u[16]*F*dt*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]+2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*dt*GCaL*Kbufss*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]))*x[15]+5.0000000000000000e-01*1.0/dt*( km1*( u[3]-u[4])/Vss*Vrel*dt*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+pow( 4.0000000000000000e+00*Kbufss*( Bufss*u[4]/( Kbufss+u[4])-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)+u[4])+pow( Bufss*u[4]/( Kbufss+u[4])-Bufss-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)-Kbufss+u[4],2.0),-(1.0/2.0))*( km1*( u[3]-u[4])*( Bufss*u[4]/( Kbufss+u[4])-Bufss-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)-Kbufss+u[4])/Vss*Vrel*dt*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+2.0000000000000000e+00*km1*( u[3]-u[4])/Vss*Vrel*dt*Kbufss*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))))*x[1]-( u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*dt*GCaL/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]+( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*dt*GCaL*Kbufss/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]+u[15]/R*( u[0]-1.5000000000000000e+01)*( Bufss*u[4]/( Kbufss+u[4])-Bufss-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)-Kbufss+u[4])/Vss*u[16]*F*dt*GCaL/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18])*pow( 4.0000000000000000e+00*Kbufss*( Bufss*u[4]/( Kbufss+u[4])-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)+u[4])+pow( Bufss*u[4]/( Kbufss+u[4])-Bufss-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)-Kbufss+u[4],2.0),-(1.0/2.0)))/dt*x[17]+5.0000000000000000e-01*x[3]*( ( km1/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+-2.0000000000000000e+00*(km1*km1)/(u[3]*u[3]*u[3])*( maxSR-minSR)*( u[3]-u[4])/Vss*Vrel*u[1]*((u[4]*u[4])*(u[4]*u[4]))/pow( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)),2.0)*Vsr/pow( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00,2.0)/pow( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00),3.0)*(EC*EC)+2.0000000000000000e+00*km1/(u[3]*u[3]*u[3])*( maxSR-minSR)*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/pow( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00,2.0)/pow( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00),2.0)*(EC*EC))*dt+( 2.0000000000000000e+00*( km1/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+-2.0000000000000000e+00*(km1*km1)/(u[3]*u[3]*u[3])*( maxSR-minSR)*( u[3]-u[4])/Vss*Vrel*u[1]*((u[4]*u[4])*(u[4]*u[4]))/pow( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)),2.0)*Vsr/pow( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00,2.0)/pow( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00),3.0)*(EC*EC)+2.0000000000000000e+00*km1/(u[3]*u[3]*u[3])*( maxSR-minSR)*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/pow( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00,2.0)/pow( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00),2.0)*(EC*EC))*dt*Kbufss+( Bufss*u[4]/( Kbufss+u[4])-Bufss-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)-Kbufss+u[4])*( km1/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+-2.0000000000000000e+00*(km1*km1)/(u[3]*u[3]*u[3])*( maxSR-minSR)*( u[3]-u[4])/Vss*Vrel*u[1]*((u[4]*u[4])*(u[4]*u[4]))/pow( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)),2.0)*Vsr/pow( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00,2.0)/pow( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00),3.0)*(EC*EC)+2.0000000000000000e+00*km1/(u[3]*u[3]*u[3])*( maxSR-minSR)*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/pow( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00,2.0)/pow( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00),2.0)*(EC*EC))*dt)*pow( 4.0000000000000000e+00*Kbufss*( Bufss*u[4]/( Kbufss+u[4])-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)+u[4])+pow( Bufss*u[4]/( Kbufss+u[4])-Bufss-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)-Kbufss+u[4],2.0),-(1.0/2.0)))/dt+5.0000000000000000e-01*( Bufss/( Kbufss+u[4])+( ( Bufss*u[4]/( Kbufss+u[4])-Bufss-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)-Kbufss+u[4])*( Bufss/( Kbufss+u[4])-( 2.0000000000000000e+00*(km1*km1)*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4]*u[4])/pow( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)),2.0)*Vsr/pow( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00),2.0)+km1/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer/Vss*Vc+5.0000000000000000e-01*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*u[18]+-2.0000000000000000e+00*km1*( u[3]-u[4])/Vss*Vrel*u[1]*u[4]/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*dt-Bufss*u[4]/pow( Kbufss+u[4],2.0)+1.0000000000000000e+00)+2.0000000000000000e+00*( Bufss/( Kbufss+u[4])-( 2.0000000000000000e+00*(km1*km1)*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4]*u[4])/pow( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)),2.0)*Vsr/pow( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00),2.0)+km1/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer/Vss*Vc+5.0000000000000000e-01*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*u[18]+-2.0000000000000000e+00*km1*( u[3]-u[4])/Vss*Vrel*u[1]*u[4]/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*dt-Bufss*u[4]/pow( Kbufss+u[4],2.0)+1.0000000000000000e+00)*Kbufss)*pow( 4.0000000000000000e+00*Kbufss*( Bufss*u[4]/( Kbufss+u[4])-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)+u[4])+pow( Bufss*u[4]/( Kbufss+u[4])-Bufss-dt*( 2.0000000000000000e+00*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*( 2.5000000000000000e-01*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[4]-Cao)*u[18]-km1*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer*( u[4]-u[2])/Vss*Vc)-Kbufss+u[4],2.0),-(1.0/2.0))-( 2.0000000000000000e+00*(km1*km1)*( u[3]-u[4])/Vss*Vrel*u[1]*(u[4]*u[4]*u[4])/pow( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)),2.0)*Vsr/pow( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00),2.0)+km1/Vss*Vrel*u[1]*(u[4]*u[4])/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00))+Vxfer/Vss*Vc+5.0000000000000000e-01*u[15]/R*( u[0]-1.5000000000000000e+01)/Vss*u[16]*F*GCaL*exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)*u[17]/T*Cm/( exp(2.0000000000000000e+00*1.0/R*( u[0]-1.5000000000000000e+01)*F/T)-1.0000000000000000e+00)*u[18]+-2.0000000000000000e+00*km1*( u[3]-u[4])/Vss*Vrel*u[1]*u[4]/( k3+km1*(u[4]*u[4])/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*Vsr/( maxSR-( maxSR-minSR)/( 1.0/(u[3]*u[3])*(EC*EC)+1.0000000000000000e+00)))*dt-Bufss*u[4]/pow( Kbufss+u[4],2.0)-1.0000000000000000e+00)/dt*x[4];
 
97
            y[5] = -( 9.0000000000000000e+00*1.0/( KmCa+Cao)*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)*exp(1.0/R*u[0]*F*gamma/T)/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))*Cao*(u[5]*u[5])+3.0000000000000000e+00*PNaK/( 3.5299999999999998e-02*exp(-1.0/R*u[0]*F/T)+1.2450000000000000e-01*exp(-1.0000000000000001e-01*1.0/R*u[0]*F/T)+1.0000000000000000e+00)/( KmNa+u[5])*Ko/( KmK+Ko)+R*u[9]/F*(u[7]*u[7]*u[7])*T*GNa*u[8]/u[5]+R/F*GbNa*T/u[5]+-3.0000000000000000e+00*PNaK/( 3.5299999999999998e-02*exp(-1.0/R*u[0]*F/T)+1.2450000000000000e-01*exp(-1.0000000000000001e-01*1.0/R*u[0]*F/T)+1.0000000000000000e+00)/pow( KmNa+u[5],2.0)*Ko/( KmK+Ko)*u[5])/F/Vc*Cm*x[5]-x[0]*( u[9]*(u[7]*u[7]*u[7])*GNa*u[8]+-3.0000000000000000e+00*1.0/R/( KmCa+Cao)*kNaCa/pow( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00,2.0)*exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*F*( gamma-1.0000000000000000e+00)*( exp(1.0/R*u[0]*F*gamma/T)*Cao*(u[5]*u[5]*u[5])-exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*alpha*u[2]*(Nao*Nao*Nao))/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))/T*ksat+GbNa+3.0000000000000000e+00*1.0/( KmCa+Cao)*( 1.0/R*F*gamma*exp(1.0/R*u[0]*F*gamma/T)/T*Cao*(u[5]*u[5]*u[5])-1.0/R*exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*F*( gamma-1.0000000000000000e+00)/T*alpha*u[2]*(Nao*Nao*Nao))*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))+-3.0000000000000000e+00*PNaK/pow( 3.5299999999999998e-02*exp(-1.0/R*u[0]*F/T)+1.2450000000000000e-01*exp(-1.0000000000000001e-01*1.0/R*u[0]*F/T)+1.0000000000000000e+00,2.0)/( KmNa+u[5])*( -1.2450000000000001e-02*1.0/R*F/T*exp(-1.0000000000000001e-01*1.0/R*u[0]*F/T)+-3.5299999999999998e-02*1.0/R*exp(-1.0/R*u[0]*F/T)*F/T)*Ko/( KmK+Ko)*u[5])/F/Vc*Cm+u[9]/F*( R*log(1.0/u[5]*Nao)/F*T-u[0])/Vc*(u[7]*u[7]*u[7])*Cm*GNa*x[8]+3.0000000000000000e+00*u[9]/F*x[7]*( R*log(1.0/u[5]*Nao)/F*T-u[0])/Vc*(u[7]*u[7])*Cm*GNa*u[8]+x[9]/F*( R*log(1.0/u[5]*Nao)/F*T-u[0])/Vc*(u[7]*u[7]*u[7])*Cm*GNa*u[8]+3.0000000000000000e+00*1.0/( KmCa+Cao)*kNaCa/( exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)*ksat+1.0000000000000000e+00)*exp(1.0/R*u[0]*F*( gamma-1.0000000000000000e+00)/T)/F/( (KmNai*KmNai*KmNai)+(Nao*Nao*Nao))/Vc*Cm*alpha*x[2]*(Nao*Nao*Nao);
 
98
            y[6] =  -2.0000000000000000e+00*1.0/F/Vc*x[14]*Cm*u[14]*GKs*( u[0]-R*log(1.0/( u[6]+pKNa*u[5])*( pKNa*Nao+Ko))/F*T)+-4.3033148291193518e-01*u[12]*x[13]/F*( u[0]-R/F*T*log(1.0/u[6]*Ko))/Vc*sqrt(Ko)*Cm*GKr+-4.3033148291193518e-01*1.0/F*( u[0]-R/F*T*log(1.0/u[6]*Ko))*u[13]/Vc*sqrt(Ko)*Cm*GKr*x[12]-1.0/F*x[10]*( u[0]-R/F*T*log(1.0/u[6]*Ko))*Gto/Vc*Cm*u[11]-x[6]*( R/u[6]/F*Gto*u[10]*T*u[11]+4.3033148291193518e-01*R*u[12]/u[6]/F*u[13]*sqrt(Ko)*T*GKr+R/F/( u[6]+pKNa*u[5])*T*(u[14]*u[14])*GKs+R/u[6]/F*T*GpK/( exp( -1.6722408026755853e-01*u[0]+4.1806020066889626e+00)+1.0000000000000000e+00)+-2.5819888974716113e-03*R/u[6]*GK1/F*( u[0]-R/F*T*log(1.0/u[6]*Ko))*exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)*sqrt(Ko)/pow( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00,2.0)*T/( 1.0/( exp( -5.0000000000000000e-01*u[0]+5.0000000000000000e-01*R/F*T*log(1.0/u[6]*Ko))+1.0000000000000000e+00)*( 3.0000000000000000e+00*exp( 2.0000000000000001e-04*u[0]+-2.0000000000000001e-04*R/F*T*log(1.0/u[6]*Ko)+2.0000000000000000e-02)+exp( 1.0000000000000001e-01*u[0]+-1.0000000000000001e-01*R/F*T*log(1.0/u[6]*Ko)-1.0000000000000000e+00))+1.0000000000000001e-01*1.0/( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00))+4.3033148291193521e-02*R/u[6]*GK1/F*sqrt(Ko)/( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00)*T/( 1.0/( exp( -5.0000000000000000e-01*u[0]+5.0000000000000000e-01*R/F*T*log(1.0/u[6]*Ko))+1.0000000000000000e+00)*( 3.0000000000000000e+00*exp( 2.0000000000000001e-04*u[0]+-2.0000000000000001e-04*R/F*T*log(1.0/u[6]*Ko)+2.0000000000000000e-02)+exp( 1.0000000000000001e-01*u[0]+-1.0000000000000001e-01*R/F*T*log(1.0/u[6]*Ko)-1.0000000000000000e+00))+1.0000000000000001e-01*1.0/( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00))+-4.3033148291193521e-02*GK1*( u[0]-R/F*T*log(1.0/u[6]*Ko))*( 5.0000000000000000e-01*R*exp( -5.0000000000000000e-01*u[0]+5.0000000000000000e-01*R/F*T*log(1.0/u[6]*Ko))/u[6]/pow( exp( -5.0000000000000000e-01*u[0]+5.0000000000000000e-01*R/F*T*log(1.0/u[6]*Ko))+1.0000000000000000e+00,2.0)*( 3.0000000000000000e+00*exp( 2.0000000000000001e-04*u[0]+-2.0000000000000001e-04*R/F*T*log(1.0/u[6]*Ko)+2.0000000000000000e-02)+exp( 1.0000000000000001e-01*u[0]+-1.0000000000000001e-01*R/F*T*log(1.0/u[6]*Ko)-1.0000000000000000e+00))/F*T+-6.0000000000000001e-03*R/u[6]/F*exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)/pow( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00,2.0)*T+1.0/( exp( -5.0000000000000000e-01*u[0]+5.0000000000000000e-01*R/F*T*log(1.0/u[6]*Ko))+1.0000000000000000e+00)*( 6.0000000000000006e-04*R*exp( 2.0000000000000001e-04*u[0]+-2.0000000000000001e-04*R/F*T*log(1.0/u[6]*Ko)+2.0000000000000000e-02)/u[6]/F*T+1.0000000000000001e-01*R/u[6]/F*T*exp( 1.0000000000000001e-01*u[0]+-1.0000000000000001e-01*R/F*T*log(1.0/u[6]*Ko)-1.0000000000000000e+00)))*sqrt(Ko)/( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00)/pow( 1.0/( exp( -5.0000000000000000e-01*u[0]+5.0000000000000000e-01*R/F*T*log(1.0/u[6]*Ko))+1.0000000000000000e+00)*( 3.0000000000000000e+00*exp( 2.0000000000000001e-04*u[0]+-2.0000000000000001e-04*R/F*T*log(1.0/u[6]*Ko)+2.0000000000000000e-02)+exp( 1.0000000000000001e-01*u[0]+-1.0000000000000001e-01*R/F*T*log(1.0/u[6]*Ko)-1.0000000000000000e+00))+1.0000000000000001e-01*1.0/( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00),2.0))/F/Vc*Cm+( 2.0000000000000000e+00*PNaK/( 3.5299999999999998e-02*exp(-1.0/R*u[0]*F/T)+1.2450000000000000e-01*exp(-1.0000000000000001e-01*1.0/R*u[0]*F/T)+1.0000000000000000e+00)/( KmNa+u[5])*Ko/( KmK+Ko)+-2.0000000000000000e+00*PNaK/( 3.5299999999999998e-02*exp(-1.0/R*u[0]*F/T)+1.2450000000000000e-01*exp(-1.0000000000000001e-01*1.0/R*u[0]*F/T)+1.0000000000000000e+00)/pow( KmNa+u[5],2.0)*Ko/( KmK+Ko)*u[5]-R*pKNa/F/( u[6]+pKNa*u[5])*T*(u[14]*u[14])*GKs)/F/Vc*Cm*x[5]-x[0]/F/Vc*Cm*( -4.3033148291193521e-02*( 1.0/( exp( -5.0000000000000000e-01*u[0]+5.0000000000000000e-01*R/F*T*log(1.0/u[6]*Ko))+1.0000000000000000e+00)*( 6.0000000000000006e-04*exp( 2.0000000000000001e-04*u[0]+-2.0000000000000001e-04*R/F*T*log(1.0/u[6]*Ko)+2.0000000000000000e-02)+1.0000000000000001e-01*exp( 1.0000000000000001e-01*u[0]+-1.0000000000000001e-01*R/F*T*log(1.0/u[6]*Ko)-1.0000000000000000e+00))+-6.0000000000000001e-03*exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)/pow( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00,2.0)+5.0000000000000000e-01*exp( -5.0000000000000000e-01*u[0]+5.0000000000000000e-01*R/F*T*log(1.0/u[6]*Ko))/pow( exp( -5.0000000000000000e-01*u[0]+5.0000000000000000e-01*R/F*T*log(1.0/u[6]*Ko))+1.0000000000000000e+00,2.0)*( 3.0000000000000000e+00*exp( 2.0000000000000001e-04*u[0]+-2.0000000000000001e-04*R/F*T*log(1.0/u[6]*Ko)+2.0000000000000000e-02)+exp( 1.0000000000000001e-01*u[0]+-1.0000000000000001e-01*R/F*T*log(1.0/u[6]*Ko)-1.0000000000000000e+00)))*GK1*( u[0]-R/F*T*log(1.0/u[6]*Ko))*sqrt(Ko)/( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00)/pow( 1.0/( exp( -5.0000000000000000e-01*u[0]+5.0000000000000000e-01*R/F*T*log(1.0/u[6]*Ko))+1.0000000000000000e+00)*( 3.0000000000000000e+00*exp( 2.0000000000000001e-04*u[0]+-2.0000000000000001e-04*R/F*T*log(1.0/u[6]*Ko)+2.0000000000000000e-02)+exp( 1.0000000000000001e-01*u[0]+-1.0000000000000001e-01*R/F*T*log(1.0/u[6]*Ko)-1.0000000000000000e+00))+1.0000000000000001e-01*1.0/( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00),2.0)+4.3033148291193521e-02*GK1*sqrt(Ko)/( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00)/( 1.0/( exp( -5.0000000000000000e-01*u[0]+5.0000000000000000e-01*R/F*T*log(1.0/u[6]*Ko))+1.0000000000000000e+00)*( 3.0000000000000000e+00*exp( 2.0000000000000001e-04*u[0]+-2.0000000000000001e-04*R/F*T*log(1.0/u[6]*Ko)+2.0000000000000000e-02)+exp( 1.0000000000000001e-01*u[0]+-1.0000000000000001e-01*R/F*T*log(1.0/u[6]*Ko)-1.0000000000000000e+00))+1.0000000000000001e-01*1.0/( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00))+4.3033148291193518e-01*u[12]*u[13]*sqrt(Ko)*GKr+Gto*u[10]*u[11]+GpK/( exp( -1.6722408026755853e-01*u[0]+4.1806020066889626e+00)+1.0000000000000000e+00)+2.0000000000000000e+00*PNaK/pow( 3.5299999999999998e-02*exp(-1.0/R*u[0]*F/T)+1.2450000000000000e-01*exp(-1.0000000000000001e-01*1.0/R*u[0]*F/T)+1.0000000000000000e+00,2.0)/( KmNa+u[5])*( -1.2450000000000001e-02*1.0/R*F/T*exp(-1.0000000000000001e-01*1.0/R*u[0]*F/T)+-3.5299999999999998e-02*1.0/R*exp(-1.0/R*u[0]*F/T)*F/T)*Ko/( KmK+Ko)*u[5]+1.6722408026755853e-01*( u[0]-R/F*T*log(1.0/u[6]*Ko))*GpK*exp( -1.6722408026755853e-01*u[0]+4.1806020066889626e+00)/pow( exp( -1.6722408026755853e-01*u[0]+4.1806020066889626e+00)+1.0000000000000000e+00,2.0)+(u[14]*u[14])*GKs+-2.5819888974716113e-03*GK1*( u[0]-R/F*T*log(1.0/u[6]*Ko))*exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)*sqrt(Ko)/pow( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00,2.0)/( 1.0/( exp( -5.0000000000000000e-01*u[0]+5.0000000000000000e-01*R/F*T*log(1.0/u[6]*Ko))+1.0000000000000000e+00)*( 3.0000000000000000e+00*exp( 2.0000000000000001e-04*u[0]+-2.0000000000000001e-04*R/F*T*log(1.0/u[6]*Ko)+2.0000000000000000e-02)+exp( 1.0000000000000001e-01*u[0]+-1.0000000000000001e-01*R/F*T*log(1.0/u[6]*Ko)-1.0000000000000000e+00))+1.0000000000000001e-01*1.0/( exp( 5.9999999999999998e-02*u[0]+-5.9999999999999998e-02*R/F*T*log(1.0/u[6]*Ko)-1.2000000000000000e+01)+1.0000000000000000e+00)))-1.0/F*( u[0]-R/F*T*log(1.0/u[6]*Ko))*Gto*u[10]/Vc*x[11]*Cm;
 
99
            y[7] = -1.0/( 1.0000000000000001e-01*1.0/( exp( 5.0000000000000001e-03*u[0]-2.5000000000000000e-01)+1.0000000000000000e+00)+1.0000000000000001e-01*1.0/( exp( 2.0000000000000001e-01*u[0]+7.0000000000000000e+00)+1.0000000000000000e+00))*x[7]*( exp( -2.0000000000000001e-01*u[0]-1.2000000000000000e+01)+1.0000000000000000e+00)+-2.0000000000000001e-01*x[0]*( -1.1074197120708749e+00*1.0/( 1.0000000000000001e-01*1.0/( exp( 5.0000000000000001e-03*u[0]-2.5000000000000000e-01)+1.0000000000000000e+00)+1.0000000000000001e-01*1.0/( exp( 2.0000000000000001e-01*u[0]+7.0000000000000000e+00)+1.0000000000000000e+00))/pow( exp( -1.1074197120708749e-01*u[0]-6.2967884828349945e+00)+1.0000000000000000e+00,3.0)*exp( -1.1074197120708749e-01*u[0]-6.2967884828349945e+00)*( exp( -2.0000000000000001e-01*u[0]-1.2000000000000000e+01)+1.0000000000000000e+00)-1.0/( 1.0000000000000001e-01*1.0/( exp( 5.0000000000000001e-03*u[0]-2.5000000000000000e-01)+1.0000000000000000e+00)+1.0000000000000001e-01*1.0/( exp( 2.0000000000000001e-01*u[0]+7.0000000000000000e+00)+1.0000000000000000e+00))*exp( -2.0000000000000001e-01*u[0]-1.2000000000000000e+01)*( u[7]-pow( exp( -1.1074197120708749e-01*u[0]-6.2967884828349945e+00)+1.0000000000000000e+00,-2.0))+-5.0000000000000000e+00*1.0/pow( 1.0000000000000001e-01*1.0/( exp( 5.0000000000000001e-03*u[0]-2.5000000000000000e-01)+1.0000000000000000e+00)+1.0000000000000001e-01*1.0/( exp( 2.0000000000000001e-01*u[0]+7.0000000000000000e+00)+1.0000000000000000e+00),2.0)*( u[7]-pow( exp( -1.1074197120708749e-01*u[0]-6.2967884828349945e+00)+1.0000000000000000e+00,-2.0))*( -2.0000000000000000e-02*1.0/pow( exp( 2.0000000000000001e-01*u[0]+7.0000000000000000e+00)+1.0000000000000000e+00,2.0)*exp( 2.0000000000000001e-01*u[0]+7.0000000000000000e+00)+-5.0000000000000001e-04*1.0/pow( exp( 5.0000000000000001e-03*u[0]-2.5000000000000000e-01)+1.0000000000000000e+00,2.0)*exp( 5.0000000000000001e-03*u[0]-2.5000000000000000e-01))*( exp( -2.0000000000000001e-01*u[0]-1.2000000000000000e+01)+1.0000000000000000e+00));
 
100
            y[8] =  ( ( -8.3823529411764710e-03*exp( -1.4705882352941177e-01*u[0]-1.1764705882352942e+01)/( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00)+7.7000000000000002e-01*1.0/( 1.3000000000000000e-01*exp( -9.0090090090090100e-02*u[0]-9.6036036036036043e-01)+1.3000000000000000e-01)/pow( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00,2.0)*exp( u[0]+4.0000000000000000e+01)+-5.7000000000000002e-02*exp( -1.4705882352941177e-01*u[0]-1.1764705882352942e+01)/pow( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00,2.0)*exp( u[0]+4.0000000000000000e+01)+( 1.0803499999999999e+05*exp(3.4849999999999998e-01*u[0])+2.1330000000000002e-01*exp(7.9000000000000001e-02*u[0]))/( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00)-( 3.1000000000000000e+05*exp(3.4849999999999998e-01*u[0])+2.7000000000000002e+00*exp(7.9000000000000001e-02*u[0]))/pow( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00,2.0)*exp( u[0]+4.0000000000000000e+01)+-9.0180180180180185e-03*1.0/pow( 1.3000000000000000e-01*exp( -9.0090090090090100e-02*u[0]-9.6036036036036043e-01)+1.3000000000000000e-01,2.0)*( 1.0/( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00)-1.0000000000000000e+00)*exp( -9.0090090090090100e-02*u[0]-9.6036036036036043e-01))*( pow( exp( 1.3458950201884254e-01*u[0]+9.6298788694481825e+00)+1.0000000000000000e+00,-2.0)-u[8])+-2.6917900403768508e-01*( 5.7000000000000002e-02*exp( -1.4705882352941177e-01*u[0]-1.1764705882352942e+01)/( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00)+-7.7000000000000002e-01*1.0/( 1.3000000000000000e-01*exp( -9.0090090090090100e-02*u[0]-9.6036036036036043e-01)+1.3000000000000000e-01)*( 1.0/( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00)-1.0000000000000000e+00)+( 3.1000000000000000e+05*exp(3.4849999999999998e-01*u[0])+2.7000000000000002e+00*exp(7.9000000000000001e-02*u[0]))/( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00))*exp( 1.3458950201884254e-01*u[0]+9.6298788694481825e+00)/pow( exp( 1.3458950201884254e-01*u[0]+9.6298788694481825e+00)+1.0000000000000000e+00,3.0))*x[0]+( -5.7000000000000002e-02*exp( -1.4705882352941177e-01*u[0]-1.1764705882352942e+01)/( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00)+7.7000000000000002e-01*1.0/( 1.3000000000000000e-01*exp( -9.0090090090090100e-02*u[0]-9.6036036036036043e-01)+1.3000000000000000e-01)*( 1.0/( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00)-1.0000000000000000e+00)-( 3.1000000000000000e+05*exp(3.4849999999999998e-01*u[0])+2.7000000000000002e+00*exp(7.9000000000000001e-02*u[0]))/( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00))*x[8];
 
101
            y[9] = -x[9]*( ( u[0]+3.7780000000000001e+01)*( -6.9480000000000002e-06*exp(-4.3909999999999998e-02*u[0])+-2.5428000000000000e+04*exp(2.4440000000000001e-01*u[0]))/( exp( 3.1100000000000000e-01*u[0]+2.4640530000000002e+01)+1.0000000000000000e+00)/( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00)+2.4240000000000001e-02*1.0/( exp( -1.3780000000000001e-01*u[0]-5.5312920000000005e+00)+1.0000000000000000e+00)/( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00)*exp(-1.0520000000000000e-02*u[0])+-5.9999999999999998e-01*( 1.0/( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00)-1.0000000000000000e+00)*exp(5.7000000000000002e-02*u[0])/( exp( -1.0000000000000001e-01*u[0]-3.2000000000000002e+00)+1.0000000000000000e+00))-x[0]*( ( u[9]-pow( exp( 1.3458950201884254e-01*u[0]+9.6298788694481825e+00)+1.0000000000000000e+00,-2.0))*( ( -6.9480000000000002e-06*exp(-4.3909999999999998e-02*u[0])+-2.5428000000000000e+04*exp(2.4440000000000001e-01*u[0]))/( exp( 3.1100000000000000e-01*u[0]+2.4640530000000002e+01)+1.0000000000000000e+00)/( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00)+3.3402720000000005e-03*exp( -1.3780000000000001e-01*u[0]-5.5312920000000005e+00)/pow( exp( -1.3780000000000001e-01*u[0]-5.5312920000000005e+00)+1.0000000000000000e+00,2.0)/( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00)*exp(-1.0520000000000000e-02*u[0])+-2.5500480000000003e-04*1.0/( exp( -1.3780000000000001e-01*u[0]-5.5312920000000005e+00)+1.0000000000000000e+00)/( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00)*exp(-1.0520000000000000e-02*u[0])+( u[0]+3.7780000000000001e+01)*( 3.0508668000000001e-07*exp(-4.3909999999999998e-02*u[0])+-6.2146032000000005e+03*exp(2.4440000000000001e-01*u[0]))/( exp( 3.1100000000000000e-01*u[0]+2.4640530000000002e+01)+1.0000000000000000e+00)/( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00)+-3.1100000000000000e-01*( u[0]+3.7780000000000001e+01)*( -6.9480000000000002e-06*exp(-4.3909999999999998e-02*u[0])+-2.5428000000000000e+04*exp(2.4440000000000001e-01*u[0]))*exp( 3.1100000000000000e-01*u[0]+2.4640530000000002e+01)/pow( exp( 3.1100000000000000e-01*u[0]+2.4640530000000002e+01)+1.0000000000000000e+00,2.0)/( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00)+5.9999999999999998e-01*exp(5.7000000000000002e-02*u[0])/( exp( -1.0000000000000001e-01*u[0]-3.2000000000000002e+00)+1.0000000000000000e+00)/pow( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00,2.0)*exp( u[0]+4.0000000000000000e+01)-( u[0]+3.7780000000000001e+01)*( -6.9480000000000002e-06*exp(-4.3909999999999998e-02*u[0])+-2.5428000000000000e+04*exp(2.4440000000000001e-01*u[0]))/( exp( 3.1100000000000000e-01*u[0]+2.4640530000000002e+01)+1.0000000000000000e+00)/pow( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00,2.0)*exp( u[0]+4.0000000000000000e+01)+-5.9999999999999998e-02*( 1.0/( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00)-1.0000000000000000e+00)*exp(5.7000000000000002e-02*u[0])/pow( exp( -1.0000000000000001e-01*u[0]-3.2000000000000002e+00)+1.0000000000000000e+00,2.0)*exp( -1.0000000000000001e-01*u[0]-3.2000000000000002e+00)+-2.4240000000000001e-02*1.0/( exp( -1.3780000000000001e-01*u[0]-5.5312920000000005e+00)+1.0000000000000000e+00)/pow( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00,2.0)*exp(-1.0520000000000000e-02*u[0])*exp( u[0]+4.0000000000000000e+01)+-3.4200000000000001e-02*( 1.0/( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00)-1.0000000000000000e+00)*exp(5.7000000000000002e-02*u[0])/( exp( -1.0000000000000001e-01*u[0]-3.2000000000000002e+00)+1.0000000000000000e+00))+2.6917900403768508e-01*( ( u[0]+3.7780000000000001e+01)*( -6.9480000000000002e-06*exp(-4.3909999999999998e-02*u[0])+-2.5428000000000000e+04*exp(2.4440000000000001e-01*u[0]))/( exp( 3.1100000000000000e-01*u[0]+2.4640530000000002e+01)+1.0000000000000000e+00)/( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00)+2.4240000000000001e-02*1.0/( exp( -1.3780000000000001e-01*u[0]-5.5312920000000005e+00)+1.0000000000000000e+00)/( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00)*exp(-1.0520000000000000e-02*u[0])+-5.9999999999999998e-01*( 1.0/( exp( u[0]+4.0000000000000000e+01)+1.0000000000000000e+00)-1.0000000000000000e+00)*exp(5.7000000000000002e-02*u[0])/( exp( -1.0000000000000001e-01*u[0]-3.2000000000000002e+00)+1.0000000000000000e+00))*exp( 1.3458950201884254e-01*u[0]+9.6298788694481825e+00)/pow( exp( 1.3458950201884254e-01*u[0]+9.6298788694481825e+00)+1.0000000000000000e+00,3.0));
 
102
            y[10] = -x[10]/( 9.5000000000000000e+00*exp(-5.5555555555555556e-04*pow( u[0]+4.0000000000000000e+01,2.0))+8.0000000000000004e-01)+-1.6666666666666666e-01*x[0]*(-1.0/( 9.5000000000000000e+00*exp(-5.5555555555555556e-04*pow( u[0]+4.0000000000000000e+01,2.0))+8.0000000000000004e-01)*exp( -1.6666666666666666e-01*u[0]+3.3333333333333335e+00)/pow( exp( -1.6666666666666666e-01*u[0]+3.3333333333333335e+00)+1.0000000000000000e+00,2.0)+6.3333333333333339e-02*1.0/pow( 9.5000000000000000e+00*exp(-5.5555555555555556e-04*pow( u[0]+4.0000000000000000e+01,2.0))+8.0000000000000004e-01,2.0)*(-1.0/( exp( -1.6666666666666666e-01*u[0]+3.3333333333333335e+00)+1.0000000000000000e+00)+u[10])*exp(-5.5555555555555556e-04*pow( u[0]+4.0000000000000000e+01,2.0))*( u[0]+4.0000000000000000e+01));
 
103
            y[11] =  -6.2500000000000003e-03*x[0]*( 3.2000000000000000e+01*1.0/pow( exp( 2.0000000000000001e-01*u[0]+4.0000000000000000e+00)+1.0000000000000000e+00,2.0)/( 5.0000000000000000e+00*1.0/( exp( 2.0000000000000001e-01*u[0]-4.0000000000000000e+00)+1.0000000000000000e+00)+8.5000000000000000e+01*exp(-3.1250000000000002e-03*pow( u[0]+4.5000000000000000e+01,2.0))+3.0000000000000000e+00)*exp( 2.0000000000000001e-01*u[0]+4.0000000000000000e+00)+-5.0000000000000000e+00*( 1.0/( exp( 2.0000000000000001e-01*u[0]+4.0000000000000000e+00)+1.0000000000000000e+00)-u[11])/pow( 5.0000000000000000e+00*1.0/( exp( 2.0000000000000001e-01*u[0]-4.0000000000000000e+00)+1.0000000000000000e+00)+8.5000000000000000e+01*exp(-3.1250000000000002e-03*pow( u[0]+4.5000000000000000e+01,2.0))+3.0000000000000000e+00,2.0)*( 1.7000000000000000e+01*exp(-3.1250000000000002e-03*pow( u[0]+4.5000000000000000e+01,2.0))*( u[0]+4.5000000000000000e+01)+3.2000000000000000e+01*exp( 2.0000000000000001e-01*u[0]-4.0000000000000000e+00)/pow( exp( 2.0000000000000001e-01*u[0]-4.0000000000000000e+00)+1.0000000000000000e+00,2.0)))-1.0/( 5.0000000000000000e+00*1.0/( exp( 2.0000000000000001e-01*u[0]-4.0000000000000000e+00)+1.0000000000000000e+00)+8.5000000000000000e+01*exp(-3.1250000000000002e-03*pow( u[0]+4.5000000000000000e+01,2.0))+3.0000000000000000e+00)*x[11];
 
104
            y[12] =  -3.7037037037037035e-04*( exp( -1.0000000000000001e-01*u[0]-4.5000000000000000e+00)+1.0000000000000000e+00)*( exp( 8.6956521739130432e-02*u[0]+2.6086956521739131e+00)+1.0000000000000000e+00)*x[12]+5.2910052910052911e-06*x[0]*( -6.0869565217391308e+00*( exp( -1.0000000000000001e-01*u[0]-4.5000000000000000e+00)+1.0000000000000000e+00)*( u[12]-1.0/( exp( -1.4285714285714285e-01*u[0]-3.7142857142857144e+00)+1.0000000000000000e+00))*exp( 8.6956521739130432e-02*u[0]+2.6086956521739131e+00)+1.0000000000000000e+01*( exp( -1.0000000000000001e-01*u[0]-4.5000000000000000e+00)+1.0000000000000000e+00)*exp( -1.4285714285714285e-01*u[0]-3.7142857142857144e+00)/pow( exp( -1.4285714285714285e-01*u[0]-3.7142857142857144e+00)+1.0000000000000000e+00,2.0)*( exp( 8.6956521739130432e-02*u[0]+2.6086956521739131e+00)+1.0000000000000000e+00)+7.0000000000000000e+00*( u[12]-1.0/( exp( -1.4285714285714285e-01*u[0]-3.7142857142857144e+00)+1.0000000000000000e+00))*( exp( 8.6956521739130432e-02*u[0]+2.6086956521739131e+00)+1.0000000000000000e+00)*exp( -1.0000000000000001e-01*u[0]-4.5000000000000000e+00));
 
105
            y[13] =  x[0]*( 1.4880952380952380e-02*( u[13]-1.0/( exp( 4.1666666666666664e-02*u[0]+3.6666666666666665e+00)+1.0000000000000000e+00))*exp( -5.0000000000000003e-02*u[0]-3.0000000000000000e+00)*( exp( 5.0000000000000003e-02*u[0]-3.0000000000000000e+00)+1.0000000000000000e+00)+-1.4880952380952380e-02*( u[13]-1.0/( exp( 4.1666666666666664e-02*u[0]+3.6666666666666665e+00)+1.0000000000000000e+00))*( exp( -5.0000000000000003e-02*u[0]-3.0000000000000000e+00)+1.0000000000000000e+00)*exp( 5.0000000000000003e-02*u[0]-3.0000000000000000e+00)+-1.2400793650793650e-02*exp( 4.1666666666666664e-02*u[0]+3.6666666666666665e+00)/pow( exp( 4.1666666666666664e-02*u[0]+3.6666666666666665e+00)+1.0000000000000000e+00,2.0)*( exp( -5.0000000000000003e-02*u[0]-3.0000000000000000e+00)+1.0000000000000000e+00)*( exp( 5.0000000000000003e-02*u[0]-3.0000000000000000e+00)+1.0000000000000000e+00))+-2.9761904761904762e-01*x[13]*( exp( -5.0000000000000003e-02*u[0]-3.0000000000000000e+00)+1.0000000000000000e+00)*( exp( 5.0000000000000003e-02*u[0]-3.0000000000000000e+00)+1.0000000000000000e+00);
 
106
            y[14] =  -2.5000000000000001e-02*x[14]/( 3.5000000000000000e+01*pow( exp( -1.6666666666666666e-01*u[0]+8.3333333333333337e-01)+1.0000000000000000e+00,-(1.0/2.0))/( exp( 6.6666666666666666e-02*u[0]-2.3333333333333335e+00)+1.0000000000000000e+00)+2.0000000000000000e+00)+2.9761904761904765e-04*x[0]*( 6.0000000000000000e+00*exp( -7.1428571428571425e-02*u[0]-3.5714285714285715e-01)/( 3.5000000000000000e+01*pow( exp( -1.6666666666666666e-01*u[0]+8.3333333333333337e-01)+1.0000000000000000e+00,-(1.0/2.0))/( exp( 6.6666666666666666e-02*u[0]-2.3333333333333335e+00)+1.0000000000000000e+00)+2.0000000000000000e+00)/pow( exp( -7.1428571428571425e-02*u[0]-3.5714285714285715e-01)+1.0000000000000000e+00,2.0)+4.9000000000000000e+01*( 1.0/( exp( -7.1428571428571425e-02*u[0]-3.5714285714285715e-01)+1.0000000000000000e+00)-u[14])*( -5.0000000000000000e+00*pow( exp( -1.6666666666666666e-01*u[0]+8.3333333333333337e-01)+1.0000000000000000e+00,-(3.0/2.0))*exp( -1.6666666666666666e-01*u[0]+8.3333333333333337e-01)/( exp( 6.6666666666666666e-02*u[0]-2.3333333333333335e+00)+1.0000000000000000e+00)+4.0000000000000000e+00*pow( exp( -1.6666666666666666e-01*u[0]+8.3333333333333337e-01)+1.0000000000000000e+00,-(1.0/2.0))*exp( 6.6666666666666666e-02*u[0]-2.3333333333333335e+00)/pow( exp( 6.6666666666666666e-02*u[0]-2.3333333333333335e+00)+1.0000000000000000e+00,2.0))/pow( 3.5000000000000000e+01*pow( exp( -1.6666666666666666e-01*u[0]+8.3333333333333337e-01)+1.0000000000000000e+00,-(1.0/2.0))/( exp( 6.6666666666666666e-02*u[0]-2.3333333333333335e+00)+1.0000000000000000e+00)+2.0000000000000000e+00,2.0));
 
107
            y[15] = -1.0/( 1.3999999999999999e+00*1.0/( exp( 2.0000000000000001e-01*u[0]+1.0000000000000000e+00)+1.0000000000000000e+00)*( 1.3999999999999999e+00*1.0/( exp( -7.6923076923076927e-02*u[0]-2.6923076923076925e+00)+1.0000000000000000e+00)+2.5000000000000000e-01)+1.0/( exp( -5.0000000000000003e-02*u[0]+2.5000000000000000e+00)+1.0000000000000000e+00))*x[15]+-5.0000000000000003e-02*x[0]*( -2.6666666666666665e+00*1.0/( 1.3999999999999999e+00*1.0/( exp( 2.0000000000000001e-01*u[0]+1.0000000000000000e+00)+1.0000000000000000e+00)*( 1.3999999999999999e+00*1.0/( exp( -7.6923076923076927e-02*u[0]-2.6923076923076925e+00)+1.0000000000000000e+00)+2.5000000000000000e-01)+1.0/( exp( -5.0000000000000003e-02*u[0]+2.5000000000000000e+00)+1.0000000000000000e+00))*exp( -1.3333333333333333e-01*u[0]-1.0666666666666667e+00)/pow( exp( -1.3333333333333333e-01*u[0]-1.0666666666666667e+00)+1.0000000000000000e+00,2.0)+1.0/pow( 1.3999999999999999e+00*1.0/( exp( 2.0000000000000001e-01*u[0]+1.0000000000000000e+00)+1.0000000000000000e+00)*( 1.3999999999999999e+00*1.0/( exp( -7.6923076923076927e-02*u[0]-2.6923076923076925e+00)+1.0000000000000000e+00)+2.5000000000000000e-01)+1.0/( exp( -5.0000000000000003e-02*u[0]+2.5000000000000000e+00)+1.0000000000000000e+00),2.0)*( -3.0153846153846149e+00*1.0/( exp( 2.0000000000000001e-01*u[0]+1.0000000000000000e+00)+1.0000000000000000e+00)*exp( -7.6923076923076927e-02*u[0]-2.6923076923076925e+00)/pow( exp( -7.6923076923076927e-02*u[0]-2.6923076923076925e+00)+1.0000000000000000e+00,2.0)-1.0/pow( exp( -5.0000000000000003e-02*u[0]+2.5000000000000000e+00)+1.0000000000000000e+00,2.0)*exp( -5.0000000000000003e-02*u[0]+2.5000000000000000e+00)+5.5999999999999996e+00*1.0/pow( exp( 2.0000000000000001e-01*u[0]+1.0000000000000000e+00)+1.0000000000000000e+00,2.0)*exp( 2.0000000000000001e-01*u[0]+1.0000000000000000e+00)*( 1.3999999999999999e+00*1.0/( exp( -7.6923076923076927e-02*u[0]-2.6923076923076925e+00)+1.0000000000000000e+00)+2.5000000000000000e-01))*( u[15]-1.0/( exp( -1.3333333333333333e-01*u[0]-1.0666666666666667e+00)+1.0000000000000000e+00)));
 
108
            y[16] =  -1.4285714285714285e-01*( 7.0000000000000000e+00*( -1.8000000000000000e+01*1.0/pow( exp( 1.0000000000000001e-01*u[0]+3.0000000000000000e+00)+1.0000000000000000e+00,2.0)*exp( 1.0000000000000001e-01*u[0]+3.0000000000000000e+00)+2.0000000000000000e+01*exp( -1.0000000000000001e-01*u[0]+1.3000000000000000e+00)/pow( exp( -1.0000000000000001e-01*u[0]+1.3000000000000000e+00)+1.0000000000000000e+00,2.0)+-9.8000000000000007e+00*exp(-4.4444444444444444e-03*pow( u[0]+2.7000000000000000e+01,2.0))*( u[0]+2.7000000000000000e+01))*( 1.0/( exp( 1.4285714285714285e-01*u[0]+2.8571428571428572e+00)+1.0000000000000000e+00)-u[16])/pow( 1.8000000000000000e+02*1.0/( exp( 1.0000000000000001e-01*u[0]+3.0000000000000000e+00)+1.0000000000000000e+00)+1.1025000000000000e+03*exp(-4.4444444444444444e-03*pow( u[0]+2.7000000000000000e+01,2.0))+2.0000000000000000e+02*1.0/( exp( -1.0000000000000001e-01*u[0]+1.3000000000000000e+00)+1.0000000000000000e+00)+2.0000000000000000e+01,2.0)+exp( 1.4285714285714285e-01*u[0]+2.8571428571428572e+00)/pow( exp( 1.4285714285714285e-01*u[0]+2.8571428571428572e+00)+1.0000000000000000e+00,2.0)/( 1.8000000000000000e+02*1.0/( exp( 1.0000000000000001e-01*u[0]+3.0000000000000000e+00)+1.0000000000000000e+00)+1.1025000000000000e+03*exp(-4.4444444444444444e-03*pow( u[0]+2.7000000000000000e+01,2.0))+2.0000000000000000e+02*1.0/( exp( -1.0000000000000001e-01*u[0]+1.3000000000000000e+00)+1.0000000000000000e+00)+2.0000000000000000e+01))*x[0]-x[16]/( 1.8000000000000000e+02*1.0/( exp( 1.0000000000000001e-01*u[0]+3.0000000000000000e+00)+1.0000000000000000e+00)+1.1025000000000000e+03*exp(-4.4444444444444444e-03*pow( u[0]+2.7000000000000000e+01,2.0))+2.0000000000000000e+02*1.0/( exp( -1.0000000000000001e-01*u[0]+1.3000000000000000e+00)+1.0000000000000000e+00)+2.0000000000000000e+01);
 
109
            y[17] = -1.0/( 8.0000000000000000e+01*1.0/( exp( 1.0000000000000001e-01*u[0]+3.0000000000000000e+00)+1.0000000000000000e+00)+3.1000000000000000e+01*1.0/( exp( -1.0000000000000001e-01*u[0]+2.5000000000000000e+00)+1.0000000000000000e+00)+5.6200000000000000e+02*exp(-4.1666666666666666e-03*pow( u[0]+2.7000000000000000e+01,2.0)))*x[17]+-1.6666666666666666e-02*x[0]*( 1.0/pow( 8.0000000000000000e+01*1.0/( exp( 1.0000000000000001e-01*u[0]+3.0000000000000000e+00)+1.0000000000000000e+00)+3.1000000000000000e+01*1.0/( exp( -1.0000000000000001e-01*u[0]+2.5000000000000000e+00)+1.0000000000000000e+00)+5.6200000000000000e+02*exp(-4.1666666666666666e-03*pow( u[0]+2.7000000000000000e+01,2.0)),2.0)*( 4.8000000000000000e+02*1.0/pow( exp( 1.0000000000000001e-01*u[0]+3.0000000000000000e+00)+1.0000000000000000e+00,2.0)*exp( 1.0000000000000001e-01*u[0]+3.0000000000000000e+00)+-1.8600000000000000e+02*exp( -1.0000000000000001e-01*u[0]+2.5000000000000000e+00)/pow( exp( -1.0000000000000001e-01*u[0]+2.5000000000000000e+00)+1.0000000000000000e+00,2.0)+2.8100000000000000e+02*( u[0]+2.7000000000000000e+01)*exp(-4.1666666666666666e-03*pow( u[0]+2.7000000000000000e+01,2.0)))*( u[17]+-6.7000000000000004e-01*1.0/( exp( 1.4285714285714285e-01*u[0]+5.0000000000000000e+00)+1.0000000000000000e+00)-3.3000000000000002e-01)+5.7428571428571429e+00*1.0/( 8.0000000000000000e+01*1.0/( exp( 1.0000000000000001e-01*u[0]+3.0000000000000000e+00)+1.0000000000000000e+00)+3.1000000000000000e+01*1.0/( exp( -1.0000000000000001e-01*u[0]+2.5000000000000000e+00)+1.0000000000000000e+00)+5.6200000000000000e+02*exp(-4.1666666666666666e-03*pow( u[0]+2.7000000000000000e+01,2.0)))*exp( 1.4285714285714285e-01*u[0]+5.0000000000000000e+00)/pow( exp( 1.4285714285714285e-01*u[0]+5.0000000000000000e+00)+1.0000000000000000e+00,2.0));
 
110
            y[18] =  -5.0000000000000000e-01*x[18]/( 4.0000000000000000e+01*1.0/( 3.9999999999999994e+02*(u[4]*u[4])+1.0000000000000000e+00)+1.0000000000000000e+00)+x[4]*( -1.5999999999999998e+04*1.0/pow( 4.0000000000000000e+01*1.0/( 3.9999999999999994e+02*(u[4]*u[4])+1.0000000000000000e+00)+1.0000000000000000e+00,2.0)/pow( 3.9999999999999994e+02*(u[4]*u[4])+1.0000000000000000e+00,2.0)*u[4]*( -5.9999999999999998e-01*1.0/( 3.9999999999999994e+02*(u[4]*u[4])+1.0000000000000000e+00)+u[18]-4.0000000000000002e-01)+-2.3999999999999997e+02*1.0/( 4.0000000000000000e+01*1.0/( 3.9999999999999994e+02*(u[4]*u[4])+1.0000000000000000e+00)+1.0000000000000000e+00)/pow( 3.9999999999999994e+02*(u[4]*u[4])+1.0000000000000000e+00,2.0)*u[4]); 
 
111
        }
 
112
 
 
113
        bool update(const double* u, double t, bool end)
 
114
        {
 
115
            return true;
 
116
        }
 
117
 
 
118
        void u0(double* u) {
 
119
            u[0] = -8.6200000000000003e+01;
 
120
            u[1] = 1.0000000000000000e+00;
 
121
            u[2] = 6.9999999999999994e-05;
 
122
            u[3] = 1.3000000000000000e+00;
 
123
            u[4] = 6.9999999999999994e-05;
 
124
            u[5] = 7.6699999999999999e+00;
 
125
            u[6] = 1.3830000000000001e+02;
 
126
            u[7] = 0.0000000000000000e+00;
 
127
            u[8] = 7.5000000000000000e-01;
 
128
            u[9] = 7.5000000000000000e-01;
 
129
            u[10] = 0.0000000000000000e+00;
 
130
            u[11] = 1.0000000000000000e+00;
 
131
            u[12] = 0.0000000000000000e+00;
 
132
            u[13] = 1.0000000000000000e+00;
 
133
            u[14] = 0.0000000000000000e+00;
 
134
            u[15] = 0.0000000000000000e+00;
 
135
            u[16] = 1.0000000000000000e+00;
 
136
            u[17] = 1.0000000000000000e+00;
 
137
            u[18] = 1.0000000000000000e+00;
 
138
 
 
139
            // Initial kick
 
140
            u[0] = -25.0;
 
141
        }
 
142
};