-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy path__d2PC12_MA__.c
4419 lines (4325 loc) · 218 KB
/
__d2PC12_MA__.c
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
#include <math.h>
#include <stdio.h>
#include <float.h>
#define exponentiale M_E
#define pi M_PI
double max(double a, double b){
return a > b ? a : b;}
double min(double a, double b){
return a < b ? a : b;}
double root(double n,double x);
double root_0(double n,double x);
double root_1(double n,double x);
double cot(double x);
double cot_0(double x);
double arccot(double x);
double arccot_0(double x);
double coth(double x);
double coth_0(double x);
double csc(double x);
double csc_0(double x);
double arccsc(double x);
double arccsc_0(double x);
double csch(double x);
double csch_0(double x);
double sec(double x);
double sec_0(double x);
double arcsec(double x);
double arcsec_0(double x);
double sech(double x);
double sech_0(double x);
void var_types( int *types);
void res_function(double t, double *dynamicVars, double *yprime, double *errors, double *constants);
void jac_function(double t, double *dynamicVars, double *yprime, double *pd, double cj, double *constants);
void ic_function(double *dynamicVars, double *constants);
double root(double n,double x){
return pow(x, 1.0/n);
}
double root_0(double n,double x){
return -(log(x)*pow(x, 1.0/n)*1.0/pow(n, 2.0));
}
double root_1(double n,double x){
return pow(x, 1.0/n - 1.0)/n;
}
double cot(double x){
return 1.0/tan(x);
}
double cot_0(double x){
return -(1.0/(pow(cos(x), 2.0)*pow(tan(x), 2.0)));
}
double arccot(double x){
return atan(1.0/x);
}
double arccot_0(double x){
return -(1.0/pow(x, 2.0)/(pow(1.0/x, 2.0) + 1.0));
}
double coth(double x){
return 1.0/tanh(x);
}
double coth_0(double x){
return -(1.0/(pow(cosh(x), 2.0)*pow(tanh(x), 2.0)));
}
double csc(double x){
return 1.0/sin(x);
}
double csc_0(double x){
return -(cos(x)/pow(sin(x), 2.0));
}
double arccsc(double x){
return asin(1.0/x);
}
double arccsc_0(double x){
return -(1.0/pow(x, 2.0)/sqrt(1.0 - pow(1.0/x, 2.0)));
}
double csch(double x){
return 1.0/sinh(x);
}
double csch_0(double x){
return -(cosh(x)/pow(sinh(x), 2.0));
}
double sec(double x){
return 1.0/cos(x);
}
double sec_0(double x){
return sin(x)/pow(cos(x), 2.0);
}
double arcsec(double x){
return acos(1.0/x);
}
double arcsec_0(double x){
return -(1.0/pow(x, 2.0)/sqrt(1.0 - pow(1.0/x, 2.0)));
}
double sech(double x){
return 1.0/cosh(x);
}
double sech_0(double x){
return -(sinh(x)/pow(cosh(x), 2.0));
}
void var_types(int *types) {
types[0] = 1;
types[1] = 1;
types[2] = 1;
types[3] = 1;
types[4] = 1;
types[5] = 1;
types[6] = 1;
types[7] = 1;
types[8] = 1;
types[9] = 1;
types[10] = 1;
types[11] = 1;
types[12] = 1;
types[13] = 1;
types[14] = 1;
types[15] = 1;
types[16] = 1;
types[17] = 1;
types[18] = 1;
types[19] = 1;
types[20] = 1;
types[21] = 1;
types[22] = 1;
types[23] = 1;
types[24] = 1;
types[25] = 1;
types[26] = 1;
types[27] = 1;
types[28] = 1;
types[29] = 1;
types[30] = 1;
types[31] = 1;
types[32] = 1;
types[33] = 1;
types[34] = 1;
types[35] = 1;
types[36] = 1;
types[37] = 1;
types[38] = 1;
types[39] = 1;
types[40] = 1;
types[41] = 1;
types[42] = 1;
types[43] = 1;
types[44] = 1;
types[45] = 1;
types[46] = 1;
types[47] = 1;
types[48] = 1;
types[49] = 1;
types[50] = 1;
types[51] = 1;
types[52] = 1;
types[53] = 1;
types[54] = 1;
types[55] = 1;
types[56] = 1;
types[57] = 1;
types[58] = 1;
types[59] = 1;
types[60] = 1;
types[61] = 1;
types[62] = 1;
types[63] = 1;
types[64] = 1;
types[65] = 1;
types[66] = 1;
types[67] = 1;
types[68] = 1;
types[69] = 1;
types[70] = 1;
types[71] = 1;
types[72] = 1;
types[73] = 1;
types[74] = 1;
types[75] = 1;
types[76] = 1;
types[77] = 1;
types[78] = 1;
types[79] = 1;
types[80] = 1;
types[81] = 1;
types[82] = 1;
types[83] = 1;
types[84] = 1;
types[85] = 1;
types[86] = 1;
types[87] = 1;
types[88] = 1;
types[89] = 1;
types[90] = 1;
types[91] = 1;
types[92] = 1;
types[93] = 1;
types[94] = 1;
types[95] = 1;
types[96] = 1;
types[97] = 1;
types[98] = 1;
types[99] = 1;
types[100] = 1;
types[101] = 1;
types[102] = 1;
types[103] = 1;
types[104] = 1;
types[105] = 1;
types[106] = 1;
types[107] = 1;
types[108] = 1;
types[109] = 1;
types[110] = 1;
types[111] = 1;
types[112] = 1;
types[113] = 1;
types[114] = 1;
types[115] = 1;
types[116] = 1;
types[117] = 1;
types[118] = 1;
types[119] = 1;
types[120] = 1;
types[121] = 1;
types[122] = 1;
types[123] = 1;
types[124] = 1;
types[125] = 1;
types[126] = 1;
types[127] = 1;
types[128] = 1;
types[129] = 1;
types[130] = 1;
types[131] = 1;
types[132] = 1;
types[133] = 1;
types[134] = 1;
types[135] = 1;
types[136] = 1;
types[137] = 1;
types[138] = 1;
types[139] = 1;
types[140] = 1;
types[141] = 1;
types[142] = 1;
types[143] = 1;
types[144] = 1;
types[145] = 1;
types[146] = 1;
types[147] = 1;
types[148] = 1;
types[149] = 1;
types[150] = 1;
types[151] = 1;
types[152] = 1;
types[153] = 1;
types[154] = 1;
types[155] = 1;
types[156] = 1;
types[157] = 1;
types[158] = 1;
types[159] = 1;
types[160] = 1;
types[161] = 1;
types[162] = 1;
types[163] = 1;
types[164] = 1;
types[165] = 1;
types[166] = 1;
types[167] = 1;
types[168] = 1;
types[169] = 1;
types[170] = 1;
types[171] = 1;
types[172] = 1;
types[173] = 1;
types[174] = 1;
types[175] = 1;
types[176] = 1;
types[177] = 1;
types[178] = 1;
types[179] = 1;
types[180] = 1;
types[181] = 1;
types[182] = 1;
types[183] = 1;
types[184] = 1;
types[185] = 1;
types[186] = 1;
types[187] = 1;
types[188] = 1;
types[189] = 1;
types[190] = 1;
types[191] = 1;
types[192] = 1;
types[193] = 1;
types[194] = 1;
types[195] = 1;
types[196] = 1;
types[197] = 1;
types[198] = 1;
types[199] = 1;
types[200] = 1;
types[201] = 1;
types[202] = 1;
types[203] = 1;
types[204] = 1;
types[205] = 1;
types[206] = 1;
types[207] = 1;
types[208] = 1;
types[209] = 1;
types[210] = 1;
types[211] = 1;
types[212] = 1;
types[213] = 1;
types[214] = 1;
types[215] = 1;
}
void res_function(double t, double *dynamicVars, double *yprime, double *errors, double *constants) {
double cell = constants[0];
double RasGapTotal = constants[1];
double RapGapTotal = constants[2];
double PP2ATotal = constants[3];
double Raf1PPtaseTotal = constants[4];
double EGF_IC = constants[5];
double NGF_IC = constants[6];
double Total_EGFReceptor = constants[7];
double Total_NGFReceptor = constants[8];
double Total_Sos = constants[9];
double Total_P90Rsk = constants[10];
double Total_Ras = constants[11];
double Total_Raf1 = constants[12];
double Total_BRaf = constants[13];
double Total_Mek = constants[14];
double Total_Erk = constants[15];
double Total_PI3K = constants[16];
double Total_Akt = constants[17];
double Total_C3G = constants[18];
double Total_Rap1 = constants[19];
double Zero = constants[20];
double k00f = constants[21];
double k00b = constants[22];
double k01f = constants[23];
double k01b = constants[24];
double k02f = constants[25];
double k02b = constants[26];
double k03f = constants[27];
double k04f = constants[28];
double k04b = constants[29];
double k05f = constants[30];
double k06f = constants[31];
double k06b = constants[32];
double k07f = constants[33];
double k08f = constants[34];
double k08b = constants[35];
double k09f = constants[36];
double k10f = constants[37];
double k10b = constants[38];
double k11f = constants[39];
double k12f = constants[40];
double k12b = constants[41];
double k13f = constants[42];
double k14f = constants[43];
double k14b = constants[44];
double k15f = constants[45];
double k16f = constants[46];
double k16b = constants[47];
double k17f = constants[48];
double k18f = constants[49];
double k18b = constants[50];
double k19f = constants[51];
double k20f = constants[52];
double k20b = constants[53];
double k21f = constants[54];
double k22f = constants[55];
double k22b = constants[56];
double k23f = constants[57];
double k24f = constants[58];
double k24b = constants[59];
double k25f = constants[60];
double k26f = constants[61];
double k26b = constants[62];
double k27f = constants[63];
double k28f = constants[64];
double k28b = constants[65];
double k29f = constants[66];
double k30f = constants[67];
double k30b = constants[68];
double k31f = constants[69];
double k32f = constants[70];
double k32b = constants[71];
double k33f = constants[72];
double k34f = constants[73];
double k34b = constants[74];
double k35f = constants[75];
double k36f = constants[76];
double k36b = constants[77];
double k37f = constants[78];
double k38f = constants[79];
double k38b = constants[80];
double k39f = constants[81];
double k40f = constants[82];
double k40b = constants[83];
double k41f = constants[84];
double k42f = constants[85];
double k42b = constants[86];
double k43f = constants[87];
double k44f = constants[88];
double k44b = constants[89];
double k45f = constants[90];
double du_cell = constants[91];
double du_RasGapTotal = constants[92];
double du_RapGapTotal = constants[93];
double du_PP2ATotal = constants[94];
double du_Raf1PPtaseTotal = constants[95];
double du_EGF_IC = constants[96];
double du_NGF_IC = constants[97];
double du_Total_EGFReceptor = constants[98];
double du_Total_NGFReceptor = constants[99];
double du_Total_Sos = constants[100];
double du_Total_P90Rsk = constants[101];
double du_Total_Ras = constants[102];
double du_Total_Raf1 = constants[103];
double du_Total_BRaf = constants[104];
double du_Total_Mek = constants[105];
double du_Total_Erk = constants[106];
double du_Total_PI3K = constants[107];
double du_Total_Akt = constants[108];
double du_Total_C3G = constants[109];
double du_Total_Rap1 = constants[110];
double du_Zero = constants[111];
double du_k00f = constants[112];
double du_k00b = constants[113];
double du_k01f = constants[114];
double du_k01b = constants[115];
double du_k02f = constants[116];
double du_k02b = constants[117];
double du_k03f = constants[118];
double du_k04f = constants[119];
double du_k04b = constants[120];
double du_k05f = constants[121];
double du_k06f = constants[122];
double du_k06b = constants[123];
double du_k07f = constants[124];
double du_k08f = constants[125];
double du_k08b = constants[126];
double du_k09f = constants[127];
double du_k10f = constants[128];
double du_k10b = constants[129];
double du_k11f = constants[130];
double du_k12f = constants[131];
double du_k12b = constants[132];
double du_k13f = constants[133];
double du_k14f = constants[134];
double du_k14b = constants[135];
double du_k15f = constants[136];
double du_k16f = constants[137];
double du_k16b = constants[138];
double du_k17f = constants[139];
double du_k18f = constants[140];
double du_k18b = constants[141];
double du_k19f = constants[142];
double du_k20f = constants[143];
double du_k20b = constants[144];
double du_k21f = constants[145];
double du_k22f = constants[146];
double du_k22b = constants[147];
double du_k23f = constants[148];
double du_k24f = constants[149];
double du_k24b = constants[150];
double du_k25f = constants[151];
double du_k26f = constants[152];
double du_k26b = constants[153];
double du_k27f = constants[154];
double du_k28f = constants[155];
double du_k28b = constants[156];
double du_k29f = constants[157];
double du_k30f = constants[158];
double du_k30b = constants[159];
double du_k31f = constants[160];
double du_k32f = constants[161];
double du_k32b = constants[162];
double du_k33f = constants[163];
double du_k34f = constants[164];
double du_k34b = constants[165];
double du_k35f = constants[166];
double du_k36f = constants[167];
double du_k36b = constants[168];
double du_k37f = constants[169];
double du_k38f = constants[170];
double du_k38b = constants[171];
double du_k39f = constants[172];
double du_k40f = constants[173];
double du_k40b = constants[174];
double du_k41f = constants[175];
double du_k42f = constants[176];
double du_k42b = constants[177];
double du_k43f = constants[178];
double du_k44f = constants[179];
double du_k44b = constants[180];
double du_k45f = constants[181];
double dv_cell = constants[182];
double dv_RasGapTotal = constants[183];
double dv_RapGapTotal = constants[184];
double dv_PP2ATotal = constants[185];
double dv_Raf1PPtaseTotal = constants[186];
double dv_EGF_IC = constants[187];
double dv_NGF_IC = constants[188];
double dv_Total_EGFReceptor = constants[189];
double dv_Total_NGFReceptor = constants[190];
double dv_Total_Sos = constants[191];
double dv_Total_P90Rsk = constants[192];
double dv_Total_Ras = constants[193];
double dv_Total_Raf1 = constants[194];
double dv_Total_BRaf = constants[195];
double dv_Total_Mek = constants[196];
double dv_Total_Erk = constants[197];
double dv_Total_PI3K = constants[198];
double dv_Total_Akt = constants[199];
double dv_Total_C3G = constants[200];
double dv_Total_Rap1 = constants[201];
double dv_Zero = constants[202];
double dv_k00f = constants[203];
double dv_k00b = constants[204];
double dv_k01f = constants[205];
double dv_k01b = constants[206];
double dv_k02f = constants[207];
double dv_k02b = constants[208];
double dv_k03f = constants[209];
double dv_k04f = constants[210];
double dv_k04b = constants[211];
double dv_k05f = constants[212];
double dv_k06f = constants[213];
double dv_k06b = constants[214];
double dv_k07f = constants[215];
double dv_k08f = constants[216];
double dv_k08b = constants[217];
double dv_k09f = constants[218];
double dv_k10f = constants[219];
double dv_k10b = constants[220];
double dv_k11f = constants[221];
double dv_k12f = constants[222];
double dv_k12b = constants[223];
double dv_k13f = constants[224];
double dv_k14f = constants[225];
double dv_k14b = constants[226];
double dv_k15f = constants[227];
double dv_k16f = constants[228];
double dv_k16b = constants[229];
double dv_k17f = constants[230];
double dv_k18f = constants[231];
double dv_k18b = constants[232];
double dv_k19f = constants[233];
double dv_k20f = constants[234];
double dv_k20b = constants[235];
double dv_k21f = constants[236];
double dv_k22f = constants[237];
double dv_k22b = constants[238];
double dv_k23f = constants[239];
double dv_k24f = constants[240];
double dv_k24b = constants[241];
double dv_k25f = constants[242];
double dv_k26f = constants[243];
double dv_k26b = constants[244];
double dv_k27f = constants[245];
double dv_k28f = constants[246];
double dv_k28b = constants[247];
double dv_k29f = constants[248];
double dv_k30f = constants[249];
double dv_k30b = constants[250];
double dv_k31f = constants[251];
double dv_k32f = constants[252];
double dv_k32b = constants[253];
double dv_k33f = constants[254];
double dv_k34f = constants[255];
double dv_k34b = constants[256];
double dv_k35f = constants[257];
double dv_k36f = constants[258];
double dv_k36b = constants[259];
double dv_k37f = constants[260];
double dv_k38f = constants[261];
double dv_k38b = constants[262];
double dv_k39f = constants[263];
double dv_k40f = constants[264];
double dv_k40b = constants[265];
double dv_k41f = constants[266];
double dv_k42f = constants[267];
double dv_k42b = constants[268];
double dv_k43f = constants[269];
double dv_k44f = constants[270];
double dv_k44b = constants[271];
double dv_k45f = constants[272];
double EGF = dynamicVars[0];
double NGF = dynamicVars[1];
double boundEGFReceptor = dynamicVars[2];
double boundNGFReceptor = dynamicVars[3];
double SosActive = dynamicVars[4];
double P90RskActive = dynamicVars[5];
double RasActive = dynamicVars[6];
double Raf1Active = dynamicVars[7];
double BRafActive = dynamicVars[8];
double MekActive = dynamicVars[9];
double ErkActive = dynamicVars[10];
double PI3KActive = dynamicVars[11];
double AktActive = dynamicVars[12];
double C3GActive = dynamicVars[13];
double Rap1Active = dynamicVars[14];
double freeEGFReceptor = dynamicVars[15];
double freeNGFReceptor = dynamicVars[16];
double SosInactive = dynamicVars[17];
double P90RskInactive = dynamicVars[18];
double RasInactive = dynamicVars[19];
double Raf1Inactive = dynamicVars[20];
double BRafInactive = dynamicVars[21];
double MekInactive = dynamicVars[22];
double ErkInactive = dynamicVars[23];
double PI3KInactive = dynamicVars[24];
double AktInactive = dynamicVars[25];
double C3GInactive = dynamicVars[26];
double Rap1Inactive = dynamicVars[27];
double RasPI3K = dynamicVars[28];
double Raf1PPtaseRaf1 = dynamicVars[29];
double BRafMek = dynamicVars[30];
double boundNGFReceptorC3G = dynamicVars[31];
double P90RskSos = dynamicVars[32];
double Raf1PPtaseBRaf = dynamicVars[33];
double PP2AActive = dynamicVars[34];
double RasGapRas = dynamicVars[35];
double ErkP90Rsk = dynamicVars[36];
double RasRaf1 = dynamicVars[37];
double MekErk = dynamicVars[38];
double RapGapRap1 = dynamicVars[39];
double RapGapActive = dynamicVars[40];
double AktRaf1 = dynamicVars[41];
double Raf1PPtase = dynamicVars[42];
double C3GRap1 = dynamicVars[43];
double Raf1Mek = dynamicVars[44];
double SosRas = dynamicVars[45];
double boundEGFReceptorPI3K = dynamicVars[46];
double RasGapActive = dynamicVars[47];
double boundEGFReceptorSos = dynamicVars[48];
double boundNGFReceptorSos = dynamicVars[49];
double Rap1BRaf = dynamicVars[50];
double PP2AMek = dynamicVars[51];
double PP2AErk = dynamicVars[52];
double PI3KAkt = dynamicVars[53];
double du_EGF = dynamicVars[54];
double du_NGF = dynamicVars[55];
double du_boundEGFReceptor = dynamicVars[56];
double du_boundNGFReceptor = dynamicVars[57];
double du_SosActive = dynamicVars[58];
double du_P90RskActive = dynamicVars[59];
double du_RasActive = dynamicVars[60];
double du_Raf1Active = dynamicVars[61];
double du_BRafActive = dynamicVars[62];
double du_MekActive = dynamicVars[63];
double du_ErkActive = dynamicVars[64];
double du_PI3KActive = dynamicVars[65];
double du_AktActive = dynamicVars[66];
double du_C3GActive = dynamicVars[67];
double du_Rap1Active = dynamicVars[68];
double du_freeEGFReceptor = dynamicVars[69];
double du_freeNGFReceptor = dynamicVars[70];
double du_SosInactive = dynamicVars[71];
double du_P90RskInactive = dynamicVars[72];
double du_RasInactive = dynamicVars[73];
double du_Raf1Inactive = dynamicVars[74];
double du_BRafInactive = dynamicVars[75];
double du_MekInactive = dynamicVars[76];
double du_ErkInactive = dynamicVars[77];
double du_PI3KInactive = dynamicVars[78];
double du_AktInactive = dynamicVars[79];
double du_C3GInactive = dynamicVars[80];
double du_Rap1Inactive = dynamicVars[81];
double du_RasPI3K = dynamicVars[82];
double du_Raf1PPtaseRaf1 = dynamicVars[83];
double du_BRafMek = dynamicVars[84];
double du_boundNGFReceptorC3G = dynamicVars[85];
double du_P90RskSos = dynamicVars[86];
double du_Raf1PPtaseBRaf = dynamicVars[87];
double du_PP2AActive = dynamicVars[88];
double du_RasGapRas = dynamicVars[89];
double du_ErkP90Rsk = dynamicVars[90];
double du_RasRaf1 = dynamicVars[91];
double du_MekErk = dynamicVars[92];
double du_RapGapRap1 = dynamicVars[93];
double du_RapGapActive = dynamicVars[94];
double du_AktRaf1 = dynamicVars[95];
double du_Raf1PPtase = dynamicVars[96];
double du_C3GRap1 = dynamicVars[97];
double du_Raf1Mek = dynamicVars[98];
double du_SosRas = dynamicVars[99];
double du_boundEGFReceptorPI3K = dynamicVars[100];
double du_RasGapActive = dynamicVars[101];
double du_boundEGFReceptorSos = dynamicVars[102];
double du_boundNGFReceptorSos = dynamicVars[103];
double du_Rap1BRaf = dynamicVars[104];
double du_PP2AMek = dynamicVars[105];
double du_PP2AErk = dynamicVars[106];
double du_PI3KAkt = dynamicVars[107];
double dv_EGF = dynamicVars[108];
double dv_NGF = dynamicVars[109];
double dv_boundEGFReceptor = dynamicVars[110];
double dv_boundNGFReceptor = dynamicVars[111];
double dv_SosActive = dynamicVars[112];
double dv_P90RskActive = dynamicVars[113];
double dv_RasActive = dynamicVars[114];
double dv_Raf1Active = dynamicVars[115];
double dv_BRafActive = dynamicVars[116];
double dv_MekActive = dynamicVars[117];
double dv_ErkActive = dynamicVars[118];
double dv_PI3KActive = dynamicVars[119];
double dv_AktActive = dynamicVars[120];
double dv_C3GActive = dynamicVars[121];
double dv_Rap1Active = dynamicVars[122];
double dv_freeEGFReceptor = dynamicVars[123];
double dv_freeNGFReceptor = dynamicVars[124];
double dv_SosInactive = dynamicVars[125];
double dv_P90RskInactive = dynamicVars[126];
double dv_RasInactive = dynamicVars[127];
double dv_Raf1Inactive = dynamicVars[128];
double dv_BRafInactive = dynamicVars[129];
double dv_MekInactive = dynamicVars[130];
double dv_ErkInactive = dynamicVars[131];
double dv_PI3KInactive = dynamicVars[132];
double dv_AktInactive = dynamicVars[133];
double dv_C3GInactive = dynamicVars[134];
double dv_Rap1Inactive = dynamicVars[135];
double dv_RasPI3K = dynamicVars[136];
double dv_Raf1PPtaseRaf1 = dynamicVars[137];
double dv_BRafMek = dynamicVars[138];
double dv_boundNGFReceptorC3G = dynamicVars[139];
double dv_P90RskSos = dynamicVars[140];
double dv_Raf1PPtaseBRaf = dynamicVars[141];
double dv_PP2AActive = dynamicVars[142];
double dv_RasGapRas = dynamicVars[143];
double dv_ErkP90Rsk = dynamicVars[144];
double dv_RasRaf1 = dynamicVars[145];
double dv_MekErk = dynamicVars[146];
double dv_RapGapRap1 = dynamicVars[147];
double dv_RapGapActive = dynamicVars[148];
double dv_AktRaf1 = dynamicVars[149];
double dv_Raf1PPtase = dynamicVars[150];
double dv_C3GRap1 = dynamicVars[151];
double dv_Raf1Mek = dynamicVars[152];
double dv_SosRas = dynamicVars[153];
double dv_boundEGFReceptorPI3K = dynamicVars[154];
double dv_RasGapActive = dynamicVars[155];
double dv_boundEGFReceptorSos = dynamicVars[156];
double dv_boundNGFReceptorSos = dynamicVars[157];
double dv_Rap1BRaf = dynamicVars[158];
double dv_PP2AMek = dynamicVars[159];
double dv_PP2AErk = dynamicVars[160];
double dv_PI3KAkt = dynamicVars[161];
double duv_EGF = dynamicVars[162];
double duv_NGF = dynamicVars[163];
double duv_boundEGFReceptor = dynamicVars[164];
double duv_boundNGFReceptor = dynamicVars[165];
double duv_SosActive = dynamicVars[166];
double duv_P90RskActive = dynamicVars[167];
double duv_RasActive = dynamicVars[168];
double duv_Raf1Active = dynamicVars[169];
double duv_BRafActive = dynamicVars[170];
double duv_MekActive = dynamicVars[171];
double duv_ErkActive = dynamicVars[172];
double duv_PI3KActive = dynamicVars[173];
double duv_AktActive = dynamicVars[174];
double duv_C3GActive = dynamicVars[175];
double duv_Rap1Active = dynamicVars[176];
double duv_freeEGFReceptor = dynamicVars[177];
double duv_freeNGFReceptor = dynamicVars[178];
double duv_SosInactive = dynamicVars[179];
double duv_P90RskInactive = dynamicVars[180];
double duv_RasInactive = dynamicVars[181];
double duv_Raf1Inactive = dynamicVars[182];
double duv_BRafInactive = dynamicVars[183];
double duv_MekInactive = dynamicVars[184];
double duv_ErkInactive = dynamicVars[185];
double duv_PI3KInactive = dynamicVars[186];
double duv_AktInactive = dynamicVars[187];
double duv_C3GInactive = dynamicVars[188];
double duv_Rap1Inactive = dynamicVars[189];
double duv_RasPI3K = dynamicVars[190];
double duv_Raf1PPtaseRaf1 = dynamicVars[191];
double duv_BRafMek = dynamicVars[192];
double duv_boundNGFReceptorC3G = dynamicVars[193];
double duv_P90RskSos = dynamicVars[194];
double duv_Raf1PPtaseBRaf = dynamicVars[195];
double duv_PP2AActive = dynamicVars[196];
double duv_RasGapRas = dynamicVars[197];
double duv_ErkP90Rsk = dynamicVars[198];
double duv_RasRaf1 = dynamicVars[199];
double duv_MekErk = dynamicVars[200];
double duv_RapGapRap1 = dynamicVars[201];
double duv_RapGapActive = dynamicVars[202];
double duv_AktRaf1 = dynamicVars[203];
double duv_Raf1PPtase = dynamicVars[204];
double duv_C3GRap1 = dynamicVars[205];
double duv_Raf1Mek = dynamicVars[206];
double duv_SosRas = dynamicVars[207];
double duv_boundEGFReceptorPI3K = dynamicVars[208];
double duv_RasGapActive = dynamicVars[209];
double duv_boundEGFReceptorSos = dynamicVars[210];
double duv_boundNGFReceptorSos = dynamicVars[211];
double duv_Rap1BRaf = dynamicVars[212];
double duv_PP2AMek = dynamicVars[213];
double duv_PP2AErk = dynamicVars[214];
double duv_PI3KAkt = dynamicVars[215];
double EGF_prime = yprime[0];
double NGF_prime = yprime[1];
double boundEGFReceptor_prime = yprime[2];
double boundNGFReceptor_prime = yprime[3];
double SosActive_prime = yprime[4];
double P90RskActive_prime = yprime[5];
double RasActive_prime = yprime[6];
double Raf1Active_prime = yprime[7];
double BRafActive_prime = yprime[8];
double MekActive_prime = yprime[9];
double ErkActive_prime = yprime[10];
double PI3KActive_prime = yprime[11];
double AktActive_prime = yprime[12];
double C3GActive_prime = yprime[13];
double Rap1Active_prime = yprime[14];
double freeEGFReceptor_prime = yprime[15];
double freeNGFReceptor_prime = yprime[16];
double SosInactive_prime = yprime[17];
double P90RskInactive_prime = yprime[18];
double RasInactive_prime = yprime[19];
double Raf1Inactive_prime = yprime[20];
double BRafInactive_prime = yprime[21];
double MekInactive_prime = yprime[22];
double ErkInactive_prime = yprime[23];
double PI3KInactive_prime = yprime[24];
double AktInactive_prime = yprime[25];
double C3GInactive_prime = yprime[26];
double Rap1Inactive_prime = yprime[27];
double RasPI3K_prime = yprime[28];
double Raf1PPtaseRaf1_prime = yprime[29];
double BRafMek_prime = yprime[30];
double boundNGFReceptorC3G_prime = yprime[31];
double P90RskSos_prime = yprime[32];
double Raf1PPtaseBRaf_prime = yprime[33];
double PP2AActive_prime = yprime[34];
double RasGapRas_prime = yprime[35];
double ErkP90Rsk_prime = yprime[36];
double RasRaf1_prime = yprime[37];
double MekErk_prime = yprime[38];
double RapGapRap1_prime = yprime[39];
double RapGapActive_prime = yprime[40];
double AktRaf1_prime = yprime[41];
double Raf1PPtase_prime = yprime[42];
double C3GRap1_prime = yprime[43];
double Raf1Mek_prime = yprime[44];
double SosRas_prime = yprime[45];
double boundEGFReceptorPI3K_prime = yprime[46];
double RasGapActive_prime = yprime[47];
double boundEGFReceptorSos_prime = yprime[48];
double boundNGFReceptorSos_prime = yprime[49];
double Rap1BRaf_prime = yprime[50];
double PP2AMek_prime = yprime[51];
double PP2AErk_prime = yprime[52];
double PI3KAkt_prime = yprime[53];
double du_EGF_prime = yprime[54];
double du_NGF_prime = yprime[55];
double du_boundEGFReceptor_prime = yprime[56];
double du_boundNGFReceptor_prime = yprime[57];
double du_SosActive_prime = yprime[58];
double du_P90RskActive_prime = yprime[59];
double du_RasActive_prime = yprime[60];
double du_Raf1Active_prime = yprime[61];
double du_BRafActive_prime = yprime[62];
double du_MekActive_prime = yprime[63];
double du_ErkActive_prime = yprime[64];
double du_PI3KActive_prime = yprime[65];
double du_AktActive_prime = yprime[66];
double du_C3GActive_prime = yprime[67];
double du_Rap1Active_prime = yprime[68];
double du_freeEGFReceptor_prime = yprime[69];
double du_freeNGFReceptor_prime = yprime[70];
double du_SosInactive_prime = yprime[71];
double du_P90RskInactive_prime = yprime[72];
double du_RasInactive_prime = yprime[73];
double du_Raf1Inactive_prime = yprime[74];
double du_BRafInactive_prime = yprime[75];
double du_MekInactive_prime = yprime[76];
double du_ErkInactive_prime = yprime[77];
double du_PI3KInactive_prime = yprime[78];
double du_AktInactive_prime = yprime[79];
double du_C3GInactive_prime = yprime[80];
double du_Rap1Inactive_prime = yprime[81];
double du_RasPI3K_prime = yprime[82];
double du_Raf1PPtaseRaf1_prime = yprime[83];
double du_BRafMek_prime = yprime[84];
double du_boundNGFReceptorC3G_prime = yprime[85];
double du_P90RskSos_prime = yprime[86];
double du_Raf1PPtaseBRaf_prime = yprime[87];
double du_PP2AActive_prime = yprime[88];
double du_RasGapRas_prime = yprime[89];
double du_ErkP90Rsk_prime = yprime[90];
double du_RasRaf1_prime = yprime[91];
double du_MekErk_prime = yprime[92];
double du_RapGapRap1_prime = yprime[93];
double du_RapGapActive_prime = yprime[94];
double du_AktRaf1_prime = yprime[95];
double du_Raf1PPtase_prime = yprime[96];
double du_C3GRap1_prime = yprime[97];
double du_Raf1Mek_prime = yprime[98];
double du_SosRas_prime = yprime[99];
double du_boundEGFReceptorPI3K_prime = yprime[100];
double du_RasGapActive_prime = yprime[101];
double du_boundEGFReceptorSos_prime = yprime[102];
double du_boundNGFReceptorSos_prime = yprime[103];
double du_Rap1BRaf_prime = yprime[104];
double du_PP2AMek_prime = yprime[105];
double du_PP2AErk_prime = yprime[106];
double du_PI3KAkt_prime = yprime[107];
double dv_EGF_prime = yprime[108];
double dv_NGF_prime = yprime[109];
double dv_boundEGFReceptor_prime = yprime[110];
double dv_boundNGFReceptor_prime = yprime[111];
double dv_SosActive_prime = yprime[112];
double dv_P90RskActive_prime = yprime[113];
double dv_RasActive_prime = yprime[114];
double dv_Raf1Active_prime = yprime[115];
double dv_BRafActive_prime = yprime[116];
double dv_MekActive_prime = yprime[117];
double dv_ErkActive_prime = yprime[118];
double dv_PI3KActive_prime = yprime[119];
double dv_AktActive_prime = yprime[120];
double dv_C3GActive_prime = yprime[121];
double dv_Rap1Active_prime = yprime[122];
double dv_freeEGFReceptor_prime = yprime[123];
double dv_freeNGFReceptor_prime = yprime[124];
double dv_SosInactive_prime = yprime[125];
double dv_P90RskInactive_prime = yprime[126];
double dv_RasInactive_prime = yprime[127];
double dv_Raf1Inactive_prime = yprime[128];
double dv_BRafInactive_prime = yprime[129];
double dv_MekInactive_prime = yprime[130];
double dv_ErkInactive_prime = yprime[131];
double dv_PI3KInactive_prime = yprime[132];
double dv_AktInactive_prime = yprime[133];
double dv_C3GInactive_prime = yprime[134];
double dv_Rap1Inactive_prime = yprime[135];
double dv_RasPI3K_prime = yprime[136];
double dv_Raf1PPtaseRaf1_prime = yprime[137];
double dv_BRafMek_prime = yprime[138];
double dv_boundNGFReceptorC3G_prime = yprime[139];
double dv_P90RskSos_prime = yprime[140];
double dv_Raf1PPtaseBRaf_prime = yprime[141];
double dv_PP2AActive_prime = yprime[142];
double dv_RasGapRas_prime = yprime[143];
double dv_ErkP90Rsk_prime = yprime[144];
double dv_RasRaf1_prime = yprime[145];
double dv_MekErk_prime = yprime[146];
double dv_RapGapRap1_prime = yprime[147];
double dv_RapGapActive_prime = yprime[148];