## Set file for gap.ode on Tue Nov 3 12:28:02 2020 391 Number of equations and auxiliaries 43 Number of parameters # Numerical stuff 1 nout 40 nullcline mesh 8 Qual RK 50000 total 0.5 DeltaT 0 T0 0 Transient 700000 Bound 1e-12 DtMin 1 DtMax 1e-07 Tolerance 0.001 Abs. Tolerance 0 Max Delay 100 Eigenvector iterates 0.001 Eigenvector tolerance 0.001 Newton tolerance 0 Poincare plane 1e-05 Boundary value tolerance 1e-05 Boundary value epsilon 20 Boundary value iterates 0 Poincare None 1 Poincare variable 1 Poincare sign 0 Stop on Section 0 Delay flag 50000 Current time 50000 Last Time 0 MyStart 1 INFLAG # Delays 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 # Bndry conds 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 # Old ICs -70.53 V1 -70.53 V2 -70.53 V3 -70.53 V4 -70.53 V5 -70.53 V6 -70.53 V7 -70.53 V8 -70.53 V9 -70.53 V10 -70.53 V11 -70.53 V12 -70.53 V13 -70.53 V14 -70.53 V15 -70.53 V16 -70.53 V17 -70.53 V18 -70.53 V19 -70.53 V20 -70.53 V21 -70.53 V22 -70.53 V23 -70.53 V24 -70.53 V25 -70.53 V26 -70.53 V27 -70.53 V28 -70.53 V29 -70.53 V30 0.25 N1 0.25 N2 0.25 N3 0.25 N4 0.25 N5 0.25 N6 0.25 N7 0.25 N8 0.25 N9 0.25 N10 0.25 N11 0.25 N12 0.25 N13 0.25 N14 0.25 N15 0.25 N16 0.25 N17 0.25 N18 0.25 N19 0.25 N20 0.25 N21 0.25 N22 0.25 N23 0.25 N24 0.25 N25 0.25 N26 0.25 N27 0.25 N28 0.25 N29 0.25 N30 0.976 HP1 0.976 HP2 0.976 HP3 0.976 HP4 0.976 HP5 0.976 HP6 0.976 HP7 0.976 HP8 0.976 HP9 0.976 HP10 0.976 HP11 0.976 HP12 0.976 HP13 0.976 HP14 0.976 HP15 0.976 HP16 0.976 HP17 0.976 HP18 0.976 HP19 0.976 HP20 0.976 HP21 0.976 HP22 0.976 HP23 0.976 HP24 0.976 HP25 0.976 HP26 0.976 HP27 0.976 HP28 0.976 HP29 0.976 HP30 5.8 KE1 5.8 KE2 5.8 KE3 5.8 KE4 5.8 KE5 5.8 KE6 5.8 KE7 5.8 KE8 5.8 KE9 5.8 KE10 5.8 KE11 5.8 KE12 5.8 KE13 5.8 KE14 5.8 KE15 5.8 KE16 5.8 KE17 5.8 KE18 5.8 KE19 5.8 KE20 5.8 KE21 5.8 KE22 5.8 KE23 5.8 KE24 5.8 KE25 5.8 KE26 5.8 KE27 5.8 KE28 5.8 KE29 5.8 KE30 3.7 NAI1 3.7 NAI2 3.7 NAI3 3.7 NAI4 3.7 NAI5 3.7 NAI6 3.7 NAI7 3.7 NAI8 3.7 NAI9 3.7 NAI10 3.7 NAI11 3.7 NAI12 3.7 NAI13 3.7 NAI14 3.7 NAI15 3.7 NAI16 3.7 NAI17 3.7 NAI18 3.7 NAI19 3.7 NAI20 3.7 NAI21 3.7 NAI22 3.7 NAI23 3.7 NAI24 3.7 NAI25 3.7 NAI26 3.7 NAI27 3.7 NAI28 3.7 NAI29 3.7 NAI30 145 NAE1 145 NAE2 145 NAE3 145 NAE4 145 NAE5 145 NAE6 145 NAE7 145 NAE8 145 NAE9 145 NAE10 145 NAE11 145 NAE12 145 NAE13 145 NAE14 145 NAE15 145 NAE16 145 NAE17 145 NAE18 145 NAE19 145 NAE20 145 NAE21 145 NAE22 145 NAE23 145 NAE24 145 NAE25 145 NAE26 145 NAE27 145 NAE28 145 NAE29 145 NAE30 137 KI1 137 KI2 137 KI3 137 KI4 137 KI5 137 KI6 137 KI7 137 KI8 137 KI9 137 KI10 137 KI11 137 KI12 137 KI13 137 KI14 137 KI15 137 KI16 137 KI17 137 KI18 137 KI19 137 KI20 137 KI21 137 KI22 137 KI23 137 KI24 137 KI25 137 KI26 137 KI27 137 KI28 137 KI29 137 KI30 0.005 EPS -81 VA1 -81 VA2 -81 VA3 -81 VA4 -81 VA5 -81 VA6 -81 VA7 -81 VA8 -81 VA9 -81 VA10 -81 VA11 -81 VA12 -81 VA13 -81 VA14 -81 VA15 -81 VA16 -81 VA17 -81 VA18 -81 VA19 -81 VA20 -81 VA21 -81 VA22 -81 VA23 -81 VA24 -81 VA25 -81 VA26 -81 VA27 -81 VA28 -81 VA29 -81 VA30 125 KIA1 125 KIA2 125 KIA3 125 KIA4 125 KIA5 125 KIA6 125 KIA7 125 KIA8 125 KIA9 125 KIA10 125 KIA11 125 KIA12 125 KIA13 125 KIA14 125 KIA15 125 KIA16 125 KIA17 125 KIA18 125 KIA19 125 KIA20 125 KIA21 125 KIA22 125 KIA23 125 KIA24 125 KIA25 125 KIA26 125 KIA27 125 KIA28 125 KIA29 125 KIA30 5.6 NAIA1 5.6 NAIA2 5.6 NAIA3 5.6 NAIA4 5.6 NAIA5 5.6 NAIA6 5.6 NAIA7 5.6 NAIA8 5.6 NAIA9 5.6 NAIA10 5.6 NAIA11 5.6 NAIA12 5.6 NAIA13 5.6 NAIA14 5.6 NAIA15 5.6 NAIA16 5.6 NAIA17 5.6 NAIA18 5.6 NAIA19 5.6 NAIA20 5.6 NAIA21 5.6 NAIA22 5.6 NAIA23 5.6 NAIA24 5.6 NAIA25 5.6 NAIA26 5.6 NAIA27 5.6 NAIA28 5.6 NAIA29 5.6 NAIA30 # Ending ICs -76.49355504704421 V1 -76.46658416879936 V2 -76.43892419986985 V3 -76.40815387389293 V4 -76.36270769945094 V5 -76.32009728980491 V6 -76.28126637824681 V7 -76.24507501235112 V8 -76.2107103983318 V9 -76.17751827422322 V10 -76.14493119625799 V11 -76.11243398849139 V12 -76.07954548005461 V13 -76.04580857242404 V14 -76.01078532601234 V15 -75.97405592000256 V16 -75.93522181170749 V17 -75.89391490167077 V18 -75.84981634970642 V19 -75.8026910176353 V20 -75.75244609606864 V21 -75.69922414683593 V22 -75.64353829000757 V23 -75.58644373551114 V24 -75.5297063607281 V25 -75.47587797989216 V26 -75.42816253033389 V27 -75.39004453744258 V28 -75.36489806820171 V29 -75.35593735081194 V30 0.1772240411293468 N1 0.1775051402284623 N2 0.1777937855532351 N3 0.1781153220645479 N4 0.1785910472668328 N5 0.1790379918998424 N6 0.1794460555640106 N7 0.179827035977719 N8 0.1801893716285006 N9 0.1805398859795554 N10 0.1808845284065511 N11 0.1812287309496677 N12 0.1815775971722277 N13 0.1819360054340584 N14 0.1823086606712029 N15 0.1827001061985584 N16 0.1831146914713534 N17 0.183556475775267 N18 0.184029027634879 N19 0.1845350537354649 N20 0.1850757616892024 N21 0.1856498403895903 N22 0.1862519657636077 N23 0.1868708847865994 N24 0.1874874991841927 N25 0.1880739441974114 N26 0.188594967752098 N27 0.1890119882419888 N28 0.1892874838102997 N29 0.1893857284944042 N30 0.704747478804617 HP1 0.7045791600074703 HP2 0.7045573655187839 HP3 0.7046440915612748 HP4 0.704686052613842 HP5 0.7047972249862647 HP6 0.7050143492549434 HP7 0.7053099321814634 HP8 0.7056624844173727 HP9 0.7060539889958218 HP10 0.706468697698481 HP11 0.7068924464299844 HP12 0.7073122186373354 HP13 0.707715832691635 HP14 0.708091693643099 HP15 0.7084285841339617 HP16 0.7087154921356362 HP17 0.7089414948714058 HP18 0.7090957464955853 HP19 0.7091676585886147 HP20 0.7091474226276961 HP21 0.709027100607391 HP22 0.7088025771480824 HP23 0.7084766335327745 HP24 0.7080630788082415 HP25 0.7075910450776109 HP26 0.7071074178540957 HP27 0.7066750352523243 HP28 0.7063660697018448 HP29 0.7062511982821967 HP30 8.735968429481968 KE1 8.739145373932898 KE2 8.748336236470028 KE3 8.763205468019654 KE4 8.783392957782301 KE5 8.808369582878166 KE6 8.837668758253269 KE7 8.870865491427237 KE8 8.907562331709476 KE9 8.947379455477053 KE10 8.989946768114523 KE11 9.034897117960865 KE12 9.081860185296531 KE13 9.130456814536185 KE14 9.180293662449598 KE15 9.23095809785152 KE16 9.282013336024516 KE17 9.332993842651202 KE18 9.383401115731127 KE19 9.432700072745599 KE20 9.4803164631708 KE21 9.525636025349609 KE22 9.568006532391118 KE23 9.606744392766947 KE24 9.641147918763206 KE25 9.670519332457639 KE26 9.694196410359012 KE27 9.711592067525359 KE28 9.722237259781444 KE29 9.725822445515416 KE30 15.96676335119309 NAI1 15.90241689442846 NAI2 15.83353971673538 NAI3 15.75482021016129 NAI4 15.64013941019014 NAI5 15.53033384652198 NAI6 15.42765336619744 NAI7 15.32973717997048 NAI8 15.23497565770341 NAI9 15.14212507062129 NAI10 15.05014535726557 NAI11 14.9581232628623 NAI12 14.86523267849653 NAI13 14.7707137376581 NAI14 14.67386290443426 NAI15 14.57403118773247 NAI16 14.47063080999493 NAI17 14.36315368740683 NAI18 14.25120865702018 NAI19 14.13458872856375 NAI20 14.01338419123376 NAI21 13.88815983779745 NAI22 13.76020859009418 NAI23 13.63186667620081 NAI24 13.50681151026952 NAI25 13.39017005862823 NAI26 13.28822750660736 NAI27 13.20770430062067 NAI28 13.15504033102693 NAI29 13.13637471025625 NAI30 55.81245399001104 NAE1 55.79966851624788 NAE2 55.76217400715286 NAE3 55.70095637164756 NAE4 55.6171905518903 NAE5 55.5125754681265 NAE6 55.38879138050511 NAE7 55.24751450921327 NAE8 55.09044202562204 NAE9 54.91930696891241 NAE10 54.73588803016656 NAE11 54.54201631558146 NAE12 54.3395801204527 NAE13 54.1305282777593 NAE14 53.9168724100254 NAE15 53.70068827407817 NAE16 53.48411628174875 NAE17 53.26936116227447 NAE18 53.0586905609023 NAE19 52.85443208946943 NAE20 52.65896789115301 NAE21 52.47472507913976 NAE22 52.30415941140721 NAE23 52.14972834885015 NAE24 52.01384861438862 NAE25 51.89883351516576 NAE26 51.8068081056668 NAE27 51.73960647890379 NAE28 51.69866245701495 NAE29 51.68490538694089 NAE30 136.7103676878165 KI1 136.7153293412226 KI2 136.7003627677181 KI3 136.6731170431222 KI4 136.645489907782 KI5 136.6028531679111 KI6 136.5672061141418 KI7 136.536304318414 KI8 136.50922783123 KI9 136.4856502395835 KI10 136.4655629424669 KI11 136.4491450817161 KI12 136.4366930179906 KI13 136.4285778892703 KI14 136.4252175134267 KI15 136.4270553982352 KI16 136.4345418263904 KI17 136.4481120558787 KI18 136.4681553421013 KI19 136.4949661040689 KI20 136.528665746398 KI21 136.5690821374307 KI22 136.6155777837079 KI23 136.6668354411463 KI24 136.720650443678 KI25 136.7738367076874 KI26 136.8223794058454 KI27 136.861880741656 KI28 136.8881260651104 KI29 136.8974533533934 KI30 0 EPS -71.14470787402092 VA1 -71.1353581162284 VA2 -71.10833425950517 VA3 -71.06468028290804 VA4 -71.00553817722147 VA5 -70.93256651398342 VA6 -70.84722773075919 VA7 -70.75087239151657 VA8 -70.64476611722772 VA9 -70.53011510825918 VA10 -70.40808611147466 VA11 -70.27982335883915 VA12 -70.14646366342272 VA13 -70.00915028891083 VA14 -69.86904593333345 VA15 -69.72734501613641 VA16 -69.58528535154498 VA17 -69.44415917899967 VA18 -69.30532335347303 VA19 -69.17020821471368 VA20 -69.04032417948392 VA21 -68.91726435159522 VA22 -68.80270036576086 VA23 -68.69836735017637 VA24 -68.60603273665286 VA25 -68.52744372898756 VA26 -68.46425116813586 VA27 -68.41791404366165 VA28 -68.38959605227086 VA29 -68.38006530260077 VA30 127.042165688339 KIA1 127.0408323156544 KIA2 127.0367931262801 KIA3 127.0298402931793 KIA4 127.0195759746313 KIA5 127.0053746655468 KIA6 126.9910201584516 KIA7 126.9765613260008 KIA8 126.9620449748297 KIA9 126.9475172435508 KIA10 126.933024514577 KIA11 126.9186140311161 KIA12 126.9043343257698 KIA13 126.8902355278288 KIA14 126.8763695976284 KIA15 126.8627905301344 KIA16 126.849554573008 KIA17 126.8367205165961 KIA18 126.8243501347941 KIA19 126.8125088872978 KIA20 126.8012670312372 KIA21 126.7907013214595 KIA22 126.7808974695412 KIA23 126.7719534126795 KIA24 126.763983089994 KIA25 126.7571196711432 KIA26 126.7515159203029 KIA27 126.7473378095481 KIA28 126.7447466026399 KIA29 126.7438667808206 KIA30 3.558651457713748 NAIA1 3.559985605629302 NAIA2 3.564027035672119 NAIA3 3.570983488319302 NAIA4 3.581252710601262 NAIA5 3.595460070087161 NAIA6 3.609821653003561 NAIA7 3.624288474701586 NAIA8 3.638813623614221 NAIA9 3.653350861120357 NAIA10 3.667853708059424 NAIA11 3.682274826355883 NAIA12 3.696565589146974 NAIA13 3.710675772348834 NAIA14 3.724553319222054 NAIA15 3.738144135772136 NAIA16 3.75139187169566 NAIA17 3.764237629505144 NAIA18 3.776619522800027 NAIA19 3.788471973297287 NAIA20 3.79972459861593 NAIA21 3.810300511831381 NAIA22 3.82011386275985 NAIA23 3.829066570334616 NAIA24 3.837044548894101 NAIA25 3.843914483907711 NAIA26 3.849523474326392 NAIA27 3.853705427095059 NAIA28 3.856298981978166 NAIA29 3.857179594031019 NAIA30 # Parameters 8310 R 310 Temp 96485 F 922 S 2160 voln 0.1 delta 1600 Sas 2000 vola 3 gna -34 thm 5 sigm 0.05 phih 0.4 gnap 10000 taubar -40 thmp 6 sigmp -48 thhp -6 sighp -49 vt 6 sig 5 gk -55 thn 14 sgn 0.05 taun0 0.27 taun1 -40 thnt -12 sn 0.8 phin 0.3 gl -70 vl 0 iapp 10 rhon 0.02 d5 0.013 d8 4.8e-06 pka 1.5e-08 pnaa 0 sigmagap 50 gkir 0.2 gamma 10 rhoa 1 cma 0 iappa 1 ncell # Graphics 0.7071097251469876 rm -0.4999999999826661 rm 0.4999958366025516 rm 0.7071038372138505 rm 0.5000041633974482 rm -0.4999999999826661 rm 0 rm 0.7071038372138505 rm 0.7071097251469876 rm 0 1 1 1 0 0 15 1 1 1 0 30 1 1 10 0 1 1 1 0 0 1 1 1 0 0 1 1 1 0 0 1 1 1 0 0 1 1 1 0 0 1 1 1 0 0 1 1 1 0 -1000 1000 0 0 3DFlag 1 Timeflag 0 Colorflag 0 Type 1 color scale 0 minscale 50000 xmax 0 xmin 0 ymax -90 ymin 12 zmax -12 zmin 25000 4e-05 -45 0.02222222222222222 0 0.08333333333333333 45 Theta 45 Phi 0 xshft 0 yshft 0 zshft 0 xlo -90 ylo 0 -90 50000 xhi 0 yhi 50000 0 # Transpose variables etc V1 2 n columns 1 n rows 1 row skip 1 col skip 1 row 0 # Coupling stuff for H funs 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 # Array plot stuff v1 30 NCols 0 Row 1 100001 NRows 1 RowSkip -90 Zmin 0 Zmax # Torus information 0 Torus flag 1=ON 6.283185307179586 Torus period # Range information R -1 eq-range stab col 0 shoot flag 1=on 10 eq-range steps 0 eq_range low 1 eq_range high V1 V1 20 Range steps 0 Cycle color 1=on 1 Reset data 1=on 1 Use old I.C.s 1=yes 0 Par1 low 0 Par2 low 1 Par1 high 1 Par2 high R 0 BVP side 0 color cycle flag 1=on 10 BVP range steps 0 BVP range low 1 BVP range high RHS etc ... dV1/dT=(-(INA(V1,N1,VNA1)+INAP(V1,HP1,VNA1)+IK(V1,N1,VK1)+IL(V1)+IPUMP1)+IAPP) dV2/dT=(-(INA(V2,N2,VNA2)+INAP(V2,HP2,VNA2)+IK(V2,N2,VK2)+IL(V2)+IPUMP2)+IAPP) dV3/dT=(-(INA(V3,N3,VNA3)+INAP(V3,HP3,VNA3)+IK(V3,N3,VK3)+IL(V3)+IPUMP3)+IAPP) dV4/dT=(-(INA(V4,N4,VNA4)+INAP(V4,HP4,VNA4)+IK(V4,N4,VK4)+IL(V4)+IPUMP4)+IAPP) dV5/dT=(-(INA(V5,N5,VNA5)+INAP(V5,HP5,VNA5)+IK(V5,N5,VK5)+IL(V5)+IPUMP5)+IAPP) dV6/dT=(-(INA(V6,N6,VNA6)+INAP(V6,HP6,VNA6)+IK(V6,N6,VK6)+IL(V6)+IPUMP6)+IAPP) dV7/dT=(-(INA(V7,N7,VNA7)+INAP(V7,HP7,VNA7)+IK(V7,N7,VK7)+IL(V7)+IPUMP7)+IAPP) dV8/dT=(-(INA(V8,N8,VNA8)+INAP(V8,HP8,VNA8)+IK(V8,N8,VK8)+IL(V8)+IPUMP8)+IAPP) dV9/dT=(-(INA(V9,N9,VNA9)+INAP(V9,HP9,VNA9)+IK(V9,N9,VK9)+IL(V9)+IPUMP9)+IAPP) dV10/dT=(-(INA(V10,N10,VNA10)+INAP(V10,HP10,VNA10)+IK(V10,N10,VK10)+IL(V10)+IPUMP10)+IAPP) dV11/dT=(-(INA(V11,N11,VNA11)+INAP(V11,HP11,VNA11)+IK(V11,N11,VK11)+IL(V11)+IPUMP11)+IAPP) dV12/dT=(-(INA(V12,N12,VNA12)+INAP(V12,HP12,VNA12)+IK(V12,N12,VK12)+IL(V12)+IPUMP12)+IAPP) dV13/dT=(-(INA(V13,N13,VNA13)+INAP(V13,HP13,VNA13)+IK(V13,N13,VK13)+IL(V13)+IPUMP13)+IAPP) dV14/dT=(-(INA(V14,N14,VNA14)+INAP(V14,HP14,VNA14)+IK(V14,N14,VK14)+IL(V14)+IPUMP14)+IAPP) dV15/dT=(-(INA(V15,N15,VNA15)+INAP(V15,HP15,VNA15)+IK(V15,N15,VK15)+IL(V15)+IPUMP15)+IAPP) dV16/dT=(-(INA(V16,N16,VNA16)+INAP(V16,HP16,VNA16)+IK(V16,N16,VK16)+IL(V16)+IPUMP16)+IAPP) dV17/dT=(-(INA(V17,N17,VNA17)+INAP(V17,HP17,VNA17)+IK(V17,N17,VK17)+IL(V17)+IPUMP17)+IAPP) dV18/dT=(-(INA(V18,N18,VNA18)+INAP(V18,HP18,VNA18)+IK(V18,N18,VK18)+IL(V18)+IPUMP18)+IAPP) dV19/dT=(-(INA(V19,N19,VNA19)+INAP(V19,HP19,VNA19)+IK(V19,N19,VK19)+IL(V19)+IPUMP19)+IAPP) dV20/dT=(-(INA(V20,N20,VNA20)+INAP(V20,HP20,VNA20)+IK(V20,N20,VK20)+IL(V20)+IPUMP20)+IAPP) dV21/dT=(-(INA(V21,N21,VNA21)+INAP(V21,HP21,VNA21)+IK(V21,N21,VK21)+IL(V21)+IPUMP21)+IAPP) dV22/dT=(-(INA(V22,N22,VNA22)+INAP(V22,HP22,VNA22)+IK(V22,N22,VK22)+IL(V22)+IPUMP22)+IAPP) dV23/dT=(-(INA(V23,N23,VNA23)+INAP(V23,HP23,VNA23)+IK(V23,N23,VK23)+IL(V23)+IPUMP23)+IAPP) dV24/dT=(-(INA(V24,N24,VNA24)+INAP(V24,HP24,VNA24)+IK(V24,N24,VK24)+IL(V24)+IPUMP24)+IAPP) dV25/dT=(-(INA(V25,N25,VNA25)+INAP(V25,HP25,VNA25)+IK(V25,N25,VK25)+IL(V25)+IPUMP25)+IAPP) dV26/dT=(-(INA(V26,N26,VNA26)+INAP(V26,HP26,VNA26)+IK(V26,N26,VK26)+IL(V26)+IPUMP26)+IAPP) dV27/dT=(-(INA(V27,N27,VNA27)+INAP(V27,HP27,VNA27)+IK(V27,N27,VK27)+IL(V27)+IPUMP27)+IAPP) dV28/dT=(-(INA(V28,N28,VNA28)+INAP(V28,HP28,VNA28)+IK(V28,N28,VK28)+IL(V28)+IPUMP28)+IAPP) dV29/dT=(-(INA(V29,N29,VNA29)+INAP(V29,HP29,VNA29)+IK(V29,N29,VK29)+IL(V29)+IPUMP29)+IAPP) dV30/dT=(-(INA(V30,N30,VNA30)+INAP(V30,HP30,VNA30)+IK(V30,N30,VK30)+IL(V30)+IPUMP30)+IAPP) dN1/dT=PHIN*(NINF(V1)-N1)/TAUN(V1) dN2/dT=PHIN*(NINF(V2)-N2)/TAUN(V2) dN3/dT=PHIN*(NINF(V3)-N3)/TAUN(V3) dN4/dT=PHIN*(NINF(V4)-N4)/TAUN(V4) dN5/dT=PHIN*(NINF(V5)-N5)/TAUN(V5) dN6/dT=PHIN*(NINF(V6)-N6)/TAUN(V6) dN7/dT=PHIN*(NINF(V7)-N7)/TAUN(V7) dN8/dT=PHIN*(NINF(V8)-N8)/TAUN(V8) dN9/dT=PHIN*(NINF(V9)-N9)/TAUN(V9) dN10/dT=PHIN*(NINF(V10)-N10)/TAUN(V10) dN11/dT=PHIN*(NINF(V11)-N11)/TAUN(V11) dN12/dT=PHIN*(NINF(V12)-N12)/TAUN(V12) dN13/dT=PHIN*(NINF(V13)-N13)/TAUN(V13) dN14/dT=PHIN*(NINF(V14)-N14)/TAUN(V14) dN15/dT=PHIN*(NINF(V15)-N15)/TAUN(V15) dN16/dT=PHIN*(NINF(V16)-N16)/TAUN(V16) dN17/dT=PHIN*(NINF(V17)-N17)/TAUN(V17) dN18/dT=PHIN*(NINF(V18)-N18)/TAUN(V18) dN19/dT=PHIN*(NINF(V19)-N19)/TAUN(V19) dN20/dT=PHIN*(NINF(V20)-N20)/TAUN(V20) dN21/dT=PHIN*(NINF(V21)-N21)/TAUN(V21) dN22/dT=PHIN*(NINF(V22)-N22)/TAUN(V22) dN23/dT=PHIN*(NINF(V23)-N23)/TAUN(V23) dN24/dT=PHIN*(NINF(V24)-N24)/TAUN(V24) dN25/dT=PHIN*(NINF(V25)-N25)/TAUN(V25) dN26/dT=PHIN*(NINF(V26)-N26)/TAUN(V26) dN27/dT=PHIN*(NINF(V27)-N27)/TAUN(V27) dN28/dT=PHIN*(NINF(V28)-N28)/TAUN(V28) dN29/dT=PHIN*(NINF(V29)-N29)/TAUN(V29) dN30/dT=PHIN*(NINF(V30)-N30)/TAUN(V30) dHP1/dT=PHIH*(HINFP(V1)-HP1)/TAUHP(V1) dHP2/dT=PHIH*(HINFP(V2)-HP2)/TAUHP(V2) dHP3/dT=PHIH*(HINFP(V3)-HP3)/TAUHP(V3) dHP4/dT=PHIH*(HINFP(V4)-HP4)/TAUHP(V4) dHP5/dT=PHIH*(HINFP(V5)-HP5)/TAUHP(V5) dHP6/dT=PHIH*(HINFP(V6)-HP6)/TAUHP(V6) dHP7/dT=PHIH*(HINFP(V7)-HP7)/TAUHP(V7) dHP8/dT=PHIH*(HINFP(V8)-HP8)/TAUHP(V8) dHP9/dT=PHIH*(HINFP(V9)-HP9)/TAUHP(V9) dHP10/dT=PHIH*(HINFP(V10)-HP10)/TAUHP(V10) dHP11/dT=PHIH*(HINFP(V11)-HP11)/TAUHP(V11) dHP12/dT=PHIH*(HINFP(V12)-HP12)/TAUHP(V12) dHP13/dT=PHIH*(HINFP(V13)-HP13)/TAUHP(V13) dHP14/dT=PHIH*(HINFP(V14)-HP14)/TAUHP(V14) dHP15/dT=PHIH*(HINFP(V15)-HP15)/TAUHP(V15) dHP16/dT=PHIH*(HINFP(V16)-HP16)/TAUHP(V16) dHP17/dT=PHIH*(HINFP(V17)-HP17)/TAUHP(V17) dHP18/dT=PHIH*(HINFP(V18)-HP18)/TAUHP(V18) dHP19/dT=PHIH*(HINFP(V19)-HP19)/TAUHP(V19) dHP20/dT=PHIH*(HINFP(V20)-HP20)/TAUHP(V20) dHP21/dT=PHIH*(HINFP(V21)-HP21)/TAUHP(V21) dHP22/dT=PHIH*(HINFP(V22)-HP22)/TAUHP(V22) dHP23/dT=PHIH*(HINFP(V23)-HP23)/TAUHP(V23) dHP24/dT=PHIH*(HINFP(V24)-HP24)/TAUHP(V24) dHP25/dT=PHIH*(HINFP(V25)-HP25)/TAUHP(V25) dHP26/dT=PHIH*(HINFP(V26)-HP26)/TAUHP(V26) dHP27/dT=PHIH*(HINFP(V27)-HP27)/TAUHP(V27) dHP28/dT=PHIH*(HINFP(V28)-HP28)/TAUHP(V28) dHP29/dT=PHIH*(HINFP(V29)-HP29)/TAUHP(V29) dHP30/dT=PHIH*(HINFP(V30)-HP30)/TAUHP(V30) dKE1/dT=EPS+C2*(IK(V1,N1,VK1)+IKPUMP1)+KDIFF1+C3E*(IKA(VA1,KE1,KIA1)-2*IPUMPA(KE1,NAIA1)) dKE2/dT=EPS+C2*(IK(V2,N2,VK2)+IKPUMP2)+KDIFF2+C3E*(IKA(VA2,KE2,KIA2)-2*IPUMPA(KE2,NAIA2)) dKE3/dT=EPS+C2*(IK(V3,N3,VK3)+IKPUMP3)+KDIFF3+C3E*(IKA(VA3,KE3,KIA3)-2*IPUMPA(KE3,NAIA3)) dKE4/dT=EPS+C2*(IK(V4,N4,VK4)+IKPUMP4)+KDIFF4+C3E*(IKA(VA4,KE4,KIA4)-2*IPUMPA(KE4,NAIA4)) dKE5/dT=EPS+C2*(IK(V5,N5,VK5)+IKPUMP5)+KDIFF5+C3E*(IKA(VA5,KE5,KIA5)-2*IPUMPA(KE5,NAIA5)) dKE6/dT=C2*(IK(V6,N6,VK6)+IKPUMP6)+KDIFF6+C3E*(IKA(VA6,KE6,KIA6)-2*IPUMPA(KE6,NAIA6)) dKE7/dT=C2*(IK(V7,N7,VK7)+IKPUMP7)+KDIFF7+C3E*(IKA(VA7,KE7,KIA7)-2*IPUMPA(KE7,NAIA7)) dKE8/dT=C2*(IK(V8,N8,VK8)+IKPUMP8)+KDIFF8+C3E*(IKA(VA8,KE8,KIA8)-2*IPUMPA(KE8,NAIA8)) dKE9/dT=C2*(IK(V9,N9,VK9)+IKPUMP9)+KDIFF9+C3E*(IKA(VA9,KE9,KIA9)-2*IPUMPA(KE9,NAIA9)) dKE10/dT=C2*(IK(V10,N10,VK10)+IKPUMP10)+KDIFF10+C3E*(IKA(VA10,KE10,KIA10)-2*IPUMPA(KE10,NAIA10)) dKE11/dT=C2*(IK(V11,N11,VK11)+IKPUMP11)+KDIFF11+C3E*(IKA(VA11,KE11,KIA11)-2*IPUMPA(KE11,NAIA11)) dKE12/dT=C2*(IK(V12,N12,VK12)+IKPUMP12)+KDIFF12+C3E*(IKA(VA12,KE12,KIA12)-2*IPUMPA(KE12,NAIA12)) dKE13/dT=C2*(IK(V13,N13,VK13)+IKPUMP13)+KDIFF13+C3E*(IKA(VA13,KE13,KIA13)-2*IPUMPA(KE13,NAIA13)) dKE14/dT=C2*(IK(V14,N14,VK14)+IKPUMP14)+KDIFF14+C3E*(IKA(VA14,KE14,KIA14)-2*IPUMPA(KE14,NAIA14)) dKE15/dT=C2*(IK(V15,N15,VK15)+IKPUMP15)+KDIFF15+C3E*(IKA(VA15,KE15,KIA15)-2*IPUMPA(KE15,NAIA15)) dKE16/dT=C2*(IK(V16,N16,VK16)+IKPUMP16)+KDIFF16+C3E*(IKA(VA16,KE16,KIA16)-2*IPUMPA(KE16,NAIA16)) dKE17/dT=C2*(IK(V17,N17,VK17)+IKPUMP17)+KDIFF17+C3E*(IKA(VA17,KE17,KIA17)-2*IPUMPA(KE17,NAIA17)) dKE18/dT=C2*(IK(V18,N18,VK18)+IKPUMP18)+KDIFF18+C3E*(IKA(VA18,KE18,KIA18)-2*IPUMPA(KE18,NAIA18)) dKE19/dT=C2*(IK(V19,N19,VK19)+IKPUMP19)+KDIFF19+C3E*(IKA(VA19,KE19,KIA19)-2*IPUMPA(KE19,NAIA19)) dKE20/dT=C2*(IK(V20,N20,VK20)+IKPUMP20)+KDIFF20+C3E*(IKA(VA20,KE20,KIA20)-2*IPUMPA(KE20,NAIA20)) dKE21/dT=C2*(IK(V21,N21,VK21)+IKPUMP21)+KDIFF21+C3E*(IKA(VA21,KE21,KIA21)-2*IPUMPA(KE21,NAIA21)) dKE22/dT=C2*(IK(V22,N22,VK22)+IKPUMP22)+KDIFF22+C3E*(IKA(VA22,KE22,KIA22)-2*IPUMPA(KE22,NAIA22)) dKE23/dT=C2*(IK(V23,N23,VK23)+IKPUMP23)+KDIFF23+C3E*(IKA(VA23,KE23,KIA23)-2*IPUMPA(KE23,NAIA23)) dKE24/dT=C2*(IK(V24,N24,VK24)+IKPUMP24)+KDIFF24+C3E*(IKA(VA24,KE24,KIA24)-2*IPUMPA(KE24,NAIA24)) dKE25/dT=C2*(IK(V25,N25,VK25)+IKPUMP25)+KDIFF25+C3E*(IKA(VA25,KE25,KIA25)-2*IPUMPA(KE25,NAIA25)) dKE26/dT=C2*(IK(V26,N26,VK26)+IKPUMP26)+KDIFF26+C3E*(IKA(VA26,KE26,KIA26)-2*IPUMPA(KE26,NAIA26)) dKE27/dT=C2*(IK(V27,N27,VK27)+IKPUMP27)+KDIFF27+C3E*(IKA(VA27,KE27,KIA27)-2*IPUMPA(KE27,NAIA27)) dKE28/dT=C2*(IK(V28,N28,VK28)+IKPUMP28)+KDIFF28+C3E*(IKA(VA28,KE28,KIA28)-2*IPUMPA(KE28,NAIA28)) dKE29/dT=C2*(IK(V29,N29,VK29)+IKPUMP29)+KDIFF29+C3E*(IKA(VA29,KE29,KIA29)-2*IPUMPA(KE29,NAIA29)) dKE30/dT=C2*(IK(V30,N30,VK30)+IKPUMP30)+KDIFF30+C3E*(IKA(VA30,KE30,KIA30)-2*IPUMPA(KE30,NAIA30)) dNAI1/dT=-C1*(INA(V1,N1,VNA1)+INAP(V1,HP1,VNA1)+INAPUMP1) dNAI2/dT=-C1*(INA(V2,N2,VNA2)+INAP(V2,HP2,VNA2)+INAPUMP2) dNAI3/dT=-C1*(INA(V3,N3,VNA3)+INAP(V3,HP3,VNA3)+INAPUMP3) dNAI4/dT=-C1*(INA(V4,N4,VNA4)+INAP(V4,HP4,VNA4)+INAPUMP4) dNAI5/dT=-C1*(INA(V5,N5,VNA5)+INAP(V5,HP5,VNA5)+INAPUMP5) dNAI6/dT=-C1*(INA(V6,N6,VNA6)+INAP(V6,HP6,VNA6)+INAPUMP6) dNAI7/dT=-C1*(INA(V7,N7,VNA7)+INAP(V7,HP7,VNA7)+INAPUMP7) dNAI8/dT=-C1*(INA(V8,N8,VNA8)+INAP(V8,HP8,VNA8)+INAPUMP8) dNAI9/dT=-C1*(INA(V9,N9,VNA9)+INAP(V9,HP9,VNA9)+INAPUMP9) dNAI10/dT=-C1*(INA(V10,N10,VNA10)+INAP(V10,HP10,VNA10)+INAPUMP10) dNAI11/dT=-C1*(INA(V11,N11,VNA11)+INAP(V11,HP11,VNA11)+INAPUMP11) dNAI12/dT=-C1*(INA(V12,N12,VNA12)+INAP(V12,HP12,VNA12)+INAPUMP12) dNAI13/dT=-C1*(INA(V13,N13,VNA13)+INAP(V13,HP13,VNA13)+INAPUMP13) dNAI14/dT=-C1*(INA(V14,N14,VNA14)+INAP(V14,HP14,VNA14)+INAPUMP14) dNAI15/dT=-C1*(INA(V15,N15,VNA15)+INAP(V15,HP15,VNA15)+INAPUMP15) dNAI16/dT=-C1*(INA(V16,N16,VNA16)+INAP(V16,HP16,VNA16)+INAPUMP16) dNAI17/dT=-C1*(INA(V17,N17,VNA17)+INAP(V17,HP17,VNA17)+INAPUMP17) dNAI18/dT=-C1*(INA(V18,N18,VNA18)+INAP(V18,HP18,VNA18)+INAPUMP18) dNAI19/dT=-C1*(INA(V19,N19,VNA19)+INAP(V19,HP19,VNA19)+INAPUMP19) dNAI20/dT=-C1*(INA(V20,N20,VNA20)+INAP(V20,HP20,VNA20)+INAPUMP20) dNAI21/dT=-C1*(INA(V21,N21,VNA21)+INAP(V21,HP21,VNA21)+INAPUMP21) dNAI22/dT=-C1*(INA(V22,N22,VNA22)+INAP(V22,HP22,VNA22)+INAPUMP22) dNAI23/dT=-C1*(INA(V23,N23,VNA23)+INAP(V23,HP23,VNA23)+INAPUMP23) dNAI24/dT=-C1*(INA(V24,N24,VNA24)+INAP(V24,HP24,VNA24)+INAPUMP24) dNAI25/dT=-C1*(INA(V25,N25,VNA25)+INAP(V25,HP25,VNA25)+INAPUMP25) dNAI26/dT=-C1*(INA(V26,N26,VNA26)+INAP(V26,HP26,VNA26)+INAPUMP26) dNAI27/dT=-C1*(INA(V27,N27,VNA27)+INAP(V27,HP27,VNA27)+INAPUMP27) dNAI28/dT=-C1*(INA(V28,N28,VNA28)+INAP(V28,HP28,VNA28)+INAPUMP28) dNAI29/dT=-C1*(INA(V29,N29,VNA29)+INAP(V29,HP29,VNA29)+INAPUMP29) dNAI30/dT=-C1*(INA(V30,N30,VNA30)+INAP(V30,HP30,VNA30)+INAPUMP30) dNAE1/dT=C2*(INA(V1,N1,VNA1)+INAP(V1,HP1,VNA1)+INAPUMP1)+NADIFF1+C3E*(INAA(VA1,NAE1,NAIA1)+3*IPUMPA(KE1,NAIA1)) dNAE2/dT=C2*(INA(V2,N2,VNA2)+INAP(V2,HP2,VNA2)+INAPUMP2)+NADIFF2+C3E*(INAA(VA2,NAE2,NAIA2)+3*IPUMPA(KE2,NAIA2)) dNAE3/dT=C2*(INA(V3,N3,VNA3)+INAP(V3,HP3,VNA3)+INAPUMP3)+NADIFF3+C3E*(INAA(VA3,NAE3,NAIA3)+3*IPUMPA(KE3,NAIA3)) dNAE4/dT=C2*(INA(V4,N4,VNA4)+INAP(V4,HP4,VNA4)+INAPUMP4)+NADIFF4+C3E*(INAA(VA4,NAE4,NAIA4)+3*IPUMPA(KE4,NAIA4)) dNAE5/dT=C2*(INA(V5,N5,VNA5)+INAP(V5,HP5,VNA5)+INAPUMP5)+NADIFF5+C3E*(INAA(VA5,NAE5,NAIA5)+3*IPUMPA(KE5,NAIA5)) dNAE6/dT=C2*(INA(V6,N6,VNA6)+INAP(V6,HP6,VNA6)+INAPUMP6)+NADIFF6+C3E*(INAA(VA6,NAE6,NAIA6)+3*IPUMPA(KE6,NAIA6)) dNAE7/dT=C2*(INA(V7,N7,VNA7)+INAP(V7,HP7,VNA7)+INAPUMP7)+NADIFF7+C3E*(INAA(VA7,NAE7,NAIA7)+3*IPUMPA(KE7,NAIA7)) dNAE8/dT=C2*(INA(V8,N8,VNA8)+INAP(V8,HP8,VNA8)+INAPUMP8)+NADIFF8+C3E*(INAA(VA8,NAE8,NAIA8)+3*IPUMPA(KE8,NAIA8)) dNAE9/dT=C2*(INA(V9,N9,VNA9)+INAP(V9,HP9,VNA9)+INAPUMP9)+NADIFF9+C3E*(INAA(VA9,NAE9,NAIA9)+3*IPUMPA(KE9,NAIA9)) dNAE10/dT=C2*(INA(V10,N10,VNA10)+INAP(V10,HP10,VNA10)+INAPUMP10)+NADIFF10+C3E*(INAA(VA10,NAE10,NAIA10)+3*IPUMPA(KE10,NAIA10)) dNAE11/dT=C2*(INA(V11,N11,VNA11)+INAP(V11,HP11,VNA11)+INAPUMP11)+NADIFF11+C3E*(INAA(VA11,NAE11,NAIA11)+3*IPUMPA(KE11,NAIA11)) dNAE12/dT=C2*(INA(V12,N12,VNA12)+INAP(V12,HP12,VNA12)+INAPUMP12)+NADIFF12+C3E*(INAA(VA12,NAE12,NAIA12)+3*IPUMPA(KE12,NAIA12)) dNAE13/dT=C2*(INA(V13,N13,VNA13)+INAP(V13,HP13,VNA13)+INAPUMP13)+NADIFF13+C3E*(INAA(VA13,NAE13,NAIA13)+3*IPUMPA(KE13,NAIA13)) dNAE14/dT=C2*(INA(V14,N14,VNA14)+INAP(V14,HP14,VNA14)+INAPUMP14)+NADIFF14+C3E*(INAA(VA14,NAE14,NAIA14)+3*IPUMPA(KE14,NAIA14)) dNAE15/dT=C2*(INA(V15,N15,VNA15)+INAP(V15,HP15,VNA15)+INAPUMP15)+NADIFF15+C3E*(INAA(VA15,NAE15,NAIA15)+3*IPUMPA(KE15,NAIA15)) dNAE16/dT=C2*(INA(V16,N16,VNA16)+INAP(V16,HP16,VNA16)+INAPUMP16)+NADIFF16+C3E*(INAA(VA16,NAE16,NAIA16)+3*IPUMPA(KE16,NAIA16)) dNAE17/dT=C2*(INA(V17,N17,VNA17)+INAP(V17,HP17,VNA17)+INAPUMP17)+NADIFF17+C3E*(INAA(VA17,NAE17,NAIA17)+3*IPUMPA(KE17,NAIA17)) dNAE18/dT=C2*(INA(V18,N18,VNA18)+INAP(V18,HP18,VNA18)+INAPUMP18)+NADIFF18+C3E*(INAA(VA18,NAE18,NAIA18)+3*IPUMPA(KE18,NAIA18)) dNAE19/dT=C2*(INA(V19,N19,VNA19)+INAP(V19,HP19,VNA19)+INAPUMP19)+NADIFF19+C3E*(INAA(VA19,NAE19,NAIA19)+3*IPUMPA(KE19,NAIA19)) dNAE20/dT=C2*(INA(V20,N20,VNA20)+INAP(V20,HP20,VNA20)+INAPUMP20)+NADIFF20+C3E*(INAA(VA20,NAE20,NAIA20)+3*IPUMPA(KE20,NAIA20)) dNAE21/dT=C2*(INA(V21,N21,VNA21)+INAP(V21,HP21,VNA21)+INAPUMP21)+NADIFF21+C3E*(INAA(VA21,NAE21,NAIA21)+3*IPUMPA(KE21,NAIA21)) dNAE22/dT=C2*(INA(V22,N22,VNA22)+INAP(V22,HP22,VNA22)+INAPUMP22)+NADIFF22+C3E*(INAA(VA22,NAE22,NAIA22)+3*IPUMPA(KE22,NAIA22)) dNAE23/dT=C2*(INA(V23,N23,VNA23)+INAP(V23,HP23,VNA23)+INAPUMP23)+NADIFF23+C3E*(INAA(VA23,NAE23,NAIA23)+3*IPUMPA(KE23,NAIA23)) dNAE24/dT=C2*(INA(V24,N24,VNA24)+INAP(V24,HP24,VNA24)+INAPUMP24)+NADIFF24+C3E*(INAA(VA24,NAE24,NAIA24)+3*IPUMPA(KE24,NAIA24)) dNAE25/dT=C2*(INA(V25,N25,VNA25)+INAP(V25,HP25,VNA25)+INAPUMP25)+NADIFF25+C3E*(INAA(VA25,NAE25,NAIA25)+3*IPUMPA(KE25,NAIA25)) dNAE26/dT=C2*(INA(V26,N26,VNA26)+INAP(V26,HP26,VNA26)+INAPUMP26)+NADIFF26+C3E*(INAA(VA26,NAE26,NAIA26)+3*IPUMPA(KE26,NAIA26)) dNAE27/dT=C2*(INA(V27,N27,VNA27)+INAP(V27,HP27,VNA27)+INAPUMP27)+NADIFF27+C3E*(INAA(VA27,NAE27,NAIA27)+3*IPUMPA(KE27,NAIA27)) dNAE28/dT=C2*(INA(V28,N28,VNA28)+INAP(V28,HP28,VNA28)+INAPUMP28)+NADIFF28+C3E*(INAA(VA28,NAE28,NAIA28)+3*IPUMPA(KE28,NAIA28)) dNAE29/dT=C2*(INA(V29,N29,VNA29)+INAP(V29,HP29,VNA29)+INAPUMP29)+NADIFF29+C3E*(INAA(VA29,NAE29,NAIA29)+3*IPUMPA(KE29,NAIA29)) dNAE30/dT=C2*(INA(V30,N30,VNA30)+INAP(V30,HP30,VNA30)+INAPUMP30)+NADIFF30+C3E*(INAA(VA30,NAE30,NAIA30)+3*IPUMPA(KE30,NAIA30)) dKI1/dT=-C1*(IK(V1,N1,VK1)+IKPUMP1) dKI2/dT=-C1*(IK(V2,N2,VK2)+IKPUMP2) dKI3/dT=-C1*(IK(V3,N3,VK3)+IKPUMP3) dKI4/dT=-C1*(IK(V4,N4,VK4)+IKPUMP4) dKI5/dT=-C1*(IK(V5,N5,VK5)+IKPUMP5) dKI6/dT=-C1*(IK(V6,N6,VK6)+IKPUMP6) dKI7/dT=-C1*(IK(V7,N7,VK7)+IKPUMP7) dKI8/dT=-C1*(IK(V8,N8,VK8)+IKPUMP8) dKI9/dT=-C1*(IK(V9,N9,VK9)+IKPUMP9) dKI10/dT=-C1*(IK(V10,N10,VK10)+IKPUMP10) dKI11/dT=-C1*(IK(V11,N11,VK11)+IKPUMP11) dKI12/dT=-C1*(IK(V12,N12,VK12)+IKPUMP12) dKI13/dT=-C1*(IK(V13,N13,VK13)+IKPUMP13) dKI14/dT=-C1*(IK(V14,N14,VK14)+IKPUMP14) dKI15/dT=-C1*(IK(V15,N15,VK15)+IKPUMP15) dKI16/dT=-C1*(IK(V16,N16,VK16)+IKPUMP16) dKI17/dT=-C1*(IK(V17,N17,VK17)+IKPUMP17) dKI18/dT=-C1*(IK(V18,N18,VK18)+IKPUMP18) dKI19/dT=-C1*(IK(V19,N19,VK19)+IKPUMP19) dKI20/dT=-C1*(IK(V20,N20,VK20)+IKPUMP20) dKI21/dT=-C1*(IK(V21,N21,VK21)+IKPUMP21) dKI22/dT=-C1*(IK(V22,N22,VK22)+IKPUMP22) dKI23/dT=-C1*(IK(V23,N23,VK23)+IKPUMP23) dKI24/dT=-C1*(IK(V24,N24,VK24)+IKPUMP24) dKI25/dT=-C1*(IK(V25,N25,VK25)+IKPUMP25) dKI26/dT=-C1*(IK(V26,N26,VK26)+IKPUMP26) dKI27/dT=-C1*(IK(V27,N27,VK27)+IKPUMP27) dKI28/dT=-C1*(IK(V28,N28,VK28)+IKPUMP28) dKI29/dT=-C1*(IK(V29,N29,VK29)+IKPUMP29) dKI30/dT=-C1*(IK(V30,N30,VK30)+IKPUMP30) dEPS/dT=0 dVA1/dT=-(IKA(VA1,KE1,KIA1)+INAA(VA1,NAE1,NAIA1)+IPUMPA(KE1,NAIA1)+IGAP1-IAPPA)/CMA dVA2/dT=-(IKA(VA2,KE2,KIA2)+INAA(VA2,NAE2,NAIA2)+IPUMPA(KE2,NAIA2)+IGAP2-IAPPA)/CMA dVA3/dT=-(IKA(VA3,KE3,KIA3)+INAA(VA3,NAE3,NAIA3)+IPUMPA(KE3,NAIA3)+IGAP3-IAPPA)/CMA dVA4/dT=-(IKA(VA4,KE4,KIA4)+INAA(VA4,NAE4,NAIA4)+IPUMPA(KE4,NAIA4)+IGAP4-IAPPA)/CMA dVA5/dT=-(IKA(VA5,KE5,KIA5)+INAA(VA5,NAE5,NAIA5)+IPUMPA(KE5,NAIA5)+IGAP5-IAPPA)/CMA dVA6/dT=-(IKA(VA6,KE6,KIA6)+INAA(VA6,NAE6,NAIA6)+IPUMPA(KE6,NAIA6)+IGAP6-IAPPA)/CMA dVA7/dT=-(IKA(VA7,KE7,KIA7)+INAA(VA7,NAE7,NAIA7)+IPUMPA(KE7,NAIA7)+IGAP7-IAPPA)/CMA dVA8/dT=-(IKA(VA8,KE8,KIA8)+INAA(VA8,NAE8,NAIA8)+IPUMPA(KE8,NAIA8)+IGAP8-IAPPA)/CMA dVA9/dT=-(IKA(VA9,KE9,KIA9)+INAA(VA9,NAE9,NAIA9)+IPUMPA(KE9,NAIA9)+IGAP9-IAPPA)/CMA dVA10/dT=-(IKA(VA10,KE10,KIA10)+INAA(VA10,NAE10,NAIA10)+IPUMPA(KE10,NAIA10)+IGAP10-IAPPA)/CMA dVA11/dT=-(IKA(VA11,KE11,KIA11)+INAA(VA11,NAE11,NAIA11)+IPUMPA(KE11,NAIA11)+IGAP11-IAPPA)/CMA dVA12/dT=-(IKA(VA12,KE12,KIA12)+INAA(VA12,NAE12,NAIA12)+IPUMPA(KE12,NAIA12)+IGAP12-IAPPA)/CMA dVA13/dT=-(IKA(VA13,KE13,KIA13)+INAA(VA13,NAE13,NAIA13)+IPUMPA(KE13,NAIA13)+IGAP13-IAPPA)/CMA dVA14/dT=-(IKA(VA14,KE14,KIA14)+INAA(VA14,NAE14,NAIA14)+IPUMPA(KE14,NAIA14)+IGAP14-IAPPA)/CMA dVA15/dT=-(IKA(VA15,KE15,KIA15)+INAA(VA15,NAE15,NAIA15)+IPUMPA(KE15,NAIA15)+IGAP15-IAPPA)/CMA dVA16/dT=-(IKA(VA16,KE16,KIA16)+INAA(VA16,NAE16,NAIA16)+IPUMPA(KE16,NAIA16)+IGAP16-IAPPA)/CMA dVA17/dT=-(IKA(VA17,KE17,KIA17)+INAA(VA17,NAE17,NAIA17)+IPUMPA(KE17,NAIA17)+IGAP17-IAPPA)/CMA dVA18/dT=-(IKA(VA18,KE18,KIA18)+INAA(VA18,NAE18,NAIA18)+IPUMPA(KE18,NAIA18)+IGAP18-IAPPA)/CMA dVA19/dT=-(IKA(VA19,KE19,KIA19)+INAA(VA19,NAE19,NAIA19)+IPUMPA(KE19,NAIA19)+IGAP19-IAPPA)/CMA dVA20/dT=-(IKA(VA20,KE20,KIA20)+INAA(VA20,NAE20,NAIA20)+IPUMPA(KE20,NAIA20)+IGAP20-IAPPA)/CMA dVA21/dT=-(IKA(VA21,KE21,KIA21)+INAA(VA21,NAE21,NAIA21)+IPUMPA(KE21,NAIA21)+IGAP21-IAPPA)/CMA dVA22/dT=-(IKA(VA22,KE22,KIA22)+INAA(VA22,NAE22,NAIA22)+IPUMPA(KE22,NAIA22)+IGAP22-IAPPA)/CMA dVA23/dT=-(IKA(VA23,KE23,KIA23)+INAA(VA23,NAE23,NAIA23)+IPUMPA(KE23,NAIA23)+IGAP23-IAPPA)/CMA dVA24/dT=-(IKA(VA24,KE24,KIA24)+INAA(VA24,NAE24,NAIA24)+IPUMPA(KE24,NAIA24)+IGAP24-IAPPA)/CMA dVA25/dT=-(IKA(VA25,KE25,KIA25)+INAA(VA25,NAE25,NAIA25)+IPUMPA(KE25,NAIA25)+IGAP25-IAPPA)/CMA dVA26/dT=-(IKA(VA26,KE26,KIA26)+INAA(VA26,NAE26,NAIA26)+IPUMPA(KE26,NAIA26)+IGAP26-IAPPA)/CMA dVA27/dT=-(IKA(VA27,KE27,KIA27)+INAA(VA27,NAE27,NAIA27)+IPUMPA(KE27,NAIA27)+IGAP27-IAPPA)/CMA dVA28/dT=-(IKA(VA28,KE28,KIA28)+INAA(VA28,NAE28,NAIA28)+IPUMPA(KE28,NAIA28)+IGAP28-IAPPA)/CMA dVA29/dT=-(IKA(VA29,KE29,KIA29)+INAA(VA29,NAE29,NAIA29)+IPUMPA(KE29,NAIA29)+IGAP29-IAPPA)/CMA dVA30/dT=-(IKA(VA30,KE30,KIA30)+INAA(VA30,NAE30,NAIA30)+IPUMPA(KE30,NAIA30)+IGAP30-IAPPA)/CMA dKIA1/dT=-C3A*(IKA(VA1,KE1,KIA1)-2*IPUMPA(KE1,NAIA1)+IKGAP1) dKIA2/dT=-C3A*(IKA(VA2,KE2,KIA2)-2*IPUMPA(KE2,NAIA2)+IKGAP2) dKIA3/dT=-C3A*(IKA(VA3,KE3,KIA3)-2*IPUMPA(KE3,NAIA3)+IKGAP3) dKIA4/dT=-C3A*(IKA(VA4,KE4,KIA4)-2*IPUMPA(KE4,NAIA4)+IKGAP4) dKIA5/dT=-C3A*(IKA(VA5,KE5,KIA5)-2*IPUMPA(KE5,NAIA5)+IKGAP5) dKIA6/dT=-C3A*(IKA(VA6,KE6,KIA6)-2*IPUMPA(KE6,NAIA6)+IKGAP6) dKIA7/dT=-C3A*(IKA(VA7,KE7,KIA7)-2*IPUMPA(KE7,NAIA7)+IKGAP7) dKIA8/dT=-C3A*(IKA(VA8,KE8,KIA8)-2*IPUMPA(KE8,NAIA8)+IKGAP8) dKIA9/dT=-C3A*(IKA(VA9,KE9,KIA9)-2*IPUMPA(KE9,NAIA9)+IKGAP9) dKIA10/dT=-C3A*(IKA(VA10,KE10,KIA10)-2*IPUMPA(KE10,NAIA10)+IKGAP10) dKIA11/dT=-C3A*(IKA(VA11,KE11,KIA11)-2*IPUMPA(KE11,NAIA11)+IKGAP11) dKIA12/dT=-C3A*(IKA(VA12,KE12,KIA12)-2*IPUMPA(KE12,NAIA12)+IKGAP12) dKIA13/dT=-C3A*(IKA(VA13,KE13,KIA13)-2*IPUMPA(KE13,NAIA13)+IKGAP13) dKIA14/dT=-C3A*(IKA(VA14,KE14,KIA14)-2*IPUMPA(KE14,NAIA14)+IKGAP14) dKIA15/dT=-C3A*(IKA(VA15,KE15,KIA15)-2*IPUMPA(KE15,NAIA15)+IKGAP15) dKIA16/dT=-C3A*(IKA(VA16,KE16,KIA16)-2*IPUMPA(KE16,NAIA16)+IKGAP16) dKIA17/dT=-C3A*(IKA(VA17,KE17,KIA17)-2*IPUMPA(KE17,NAIA17)+IKGAP17) dKIA18/dT=-C3A*(IKA(VA18,KE18,KIA18)-2*IPUMPA(KE18,NAIA18)+IKGAP18) dKIA19/dT=-C3A*(IKA(VA19,KE19,KIA19)-2*IPUMPA(KE19,NAIA19)+IKGAP19) dKIA20/dT=-C3A*(IKA(VA20,KE20,KIA20)-2*IPUMPA(KE20,NAIA20)+IKGAP20) dKIA21/dT=-C3A*(IKA(VA21,KE21,KIA21)-2*IPUMPA(KE21,NAIA21)+IKGAP21) dKIA22/dT=-C3A*(IKA(VA22,KE22,KIA22)-2*IPUMPA(KE22,NAIA22)+IKGAP22) dKIA23/dT=-C3A*(IKA(VA23,KE23,KIA23)-2*IPUMPA(KE23,NAIA23)+IKGAP23) dKIA24/dT=-C3A*(IKA(VA24,KE24,KIA24)-2*IPUMPA(KE24,NAIA24)+IKGAP24) dKIA25/dT=-C3A*(IKA(VA25,KE25,KIA25)-2*IPUMPA(KE25,NAIA25)+IKGAP25) dKIA26/dT=-C3A*(IKA(VA26,KE26,KIA26)-2*IPUMPA(KE26,NAIA26)+IKGAP26) dKIA27/dT=-C3A*(IKA(VA27,KE27,KIA27)-2*IPUMPA(KE27,NAIA27)+IKGAP27) dKIA28/dT=-C3A*(IKA(VA28,KE28,KIA28)-2*IPUMPA(KE28,NAIA28)+IKGAP28) dKIA29/dT=-C3A*(IKA(VA29,KE29,KIA29)-2*IPUMPA(KE29,NAIA29)+IKGAP29) dKIA30/dT=-C3A*(IKA(VA30,KE30,KIA30)-2*IPUMPA(KE30,NAIA30)+IKGAP30) dNAIA1/dT=-C3A*(INAA(VA1,NAE1,NAIA1)+3*IPUMPA(KE1,NAIA1)+INAGAP1) dNAIA2/dT=-C3A*(INAA(VA2,NAE2,NAIA2)+3*IPUMPA(KE2,NAIA2)+INAGAP2) dNAIA3/dT=-C3A*(INAA(VA3,NAE3,NAIA3)+3*IPUMPA(KE3,NAIA3)+INAGAP3) dNAIA4/dT=-C3A*(INAA(VA4,NAE4,NAIA4)+3*IPUMPA(KE4,NAIA4)+INAGAP4) dNAIA5/dT=-C3A*(INAA(VA5,NAE5,NAIA5)+3*IPUMPA(KE5,NAIA5)+INAGAP5) dNAIA6/dT=-C3A*(INAA(VA6,NAE6,NAIA6)+3*IPUMPA(KE6,NAIA6)+INAGAP6) dNAIA7/dT=-C3A*(INAA(VA7,NAE7,NAIA7)+3*IPUMPA(KE7,NAIA7)+INAGAP7) dNAIA8/dT=-C3A*(INAA(VA8,NAE8,NAIA8)+3*IPUMPA(KE8,NAIA8)+INAGAP8) dNAIA9/dT=-C3A*(INAA(VA9,NAE9,NAIA9)+3*IPUMPA(KE9,NAIA9)+INAGAP9) dNAIA10/dT=-C3A*(INAA(VA10,NAE10,NAIA10)+3*IPUMPA(KE10,NAIA10)+INAGAP10) dNAIA11/dT=-C3A*(INAA(VA11,NAE11,NAIA11)+3*IPUMPA(KE11,NAIA11)+INAGAP11) dNAIA12/dT=-C3A*(INAA(VA12,NAE12,NAIA12)+3*IPUMPA(KE12,NAIA12)+INAGAP12) dNAIA13/dT=-C3A*(INAA(VA13,NAE13,NAIA13)+3*IPUMPA(KE13,NAIA13)+INAGAP13) dNAIA14/dT=-C3A*(INAA(VA14,NAE14,NAIA14)+3*IPUMPA(KE14,NAIA14)+INAGAP14) dNAIA15/dT=-C3A*(INAA(VA15,NAE15,NAIA15)+3*IPUMPA(KE15,NAIA15)+INAGAP15) dNAIA16/dT=-C3A*(INAA(VA16,NAE16,NAIA16)+3*IPUMPA(KE16,NAIA16)+INAGAP16) dNAIA17/dT=-C3A*(INAA(VA17,NAE17,NAIA17)+3*IPUMPA(KE17,NAIA17)+INAGAP17) dNAIA18/dT=-C3A*(INAA(VA18,NAE18,NAIA18)+3*IPUMPA(KE18,NAIA18)+INAGAP18) dNAIA19/dT=-C3A*(INAA(VA19,NAE19,NAIA19)+3*IPUMPA(KE19,NAIA19)+INAGAP19) dNAIA20/dT=-C3A*(INAA(VA20,NAE20,NAIA20)+3*IPUMPA(KE20,NAIA20)+INAGAP20) dNAIA21/dT=-C3A*(INAA(VA21,NAE21,NAIA21)+3*IPUMPA(KE21,NAIA21)+INAGAP21) dNAIA22/dT=-C3A*(INAA(VA22,NAE22,NAIA22)+3*IPUMPA(KE22,NAIA22)+INAGAP22) dNAIA23/dT=-C3A*(INAA(VA23,NAE23,NAIA23)+3*IPUMPA(KE23,NAIA23)+INAGAP23) dNAIA24/dT=-C3A*(INAA(VA24,NAE24,NAIA24)+3*IPUMPA(KE24,NAIA24)+INAGAP24) dNAIA25/dT=-C3A*(INAA(VA25,NAE25,NAIA25)+3*IPUMPA(KE25,NAIA25)+INAGAP25) dNAIA26/dT=-C3A*(INAA(VA26,NAE26,NAIA26)+3*IPUMPA(KE26,NAIA26)+INAGAP26) dNAIA27/dT=-C3A*(INAA(VA27,NAE27,NAIA27)+3*IPUMPA(KE27,NAIA27)+INAGAP27) dNAIA28/dT=-C3A*(INAA(VA28,NAE28,NAIA28)+3*IPUMPA(KE28,NAIA28)+INAGAP28) dNAIA29/dT=-C3A*(INAA(VA29,NAE29,NAIA29)+3*IPUMPA(KE29,NAIA29)+INAGAP29) dNAIA30/dT=-C3A*(INAA(VA30,NAE30,NAIA30)+3*IPUMPA(KE30,NAIA30)+INAGAP30) IPUMP1=IPUMP1 IPUMP2=IPUMP2 IPUMP3=IPUMP3 IPUMP4=IPUMP4 IPUMP5=IPUMP5 IPUMP6=IPUMP6 IPUMP7=IPUMP7 IPUMP8=IPUMP8 IPUMP9=IPUMP9 IPUMP10=IPUMP10 IPUMP11=IPUMP11 IPUMP12=IPUMP12 IPUMP13=IPUMP13 IPUMP14=IPUMP14 IPUMP15=IPUMP15 IPUMP16=IPUMP16 IPUMP17=IPUMP17 IPUMP18=IPUMP18 IPUMP19=IPUMP19 IPUMP20=IPUMP20 IPUMP21=IPUMP21 IPUMP22=IPUMP22 IPUMP23=IPUMP23 IPUMP24=IPUMP24 IPUMP25=IPUMP25 IPUMP26=IPUMP26 IPUMP27=IPUMP27 IPUMP28=IPUMP28 IPUMP29=IPUMP29 IPUMP30=IPUMP30 EK1=EK(KE1,KIA1) EK2=EK(KE2,KIA2) EK3=EK(KE3,KIA3) EK4=EK(KE4,KIA4) EK5=EK(KE5,KIA5) EK6=EK(KE6,KIA6) EK7=EK(KE7,KIA7) EK8=EK(KE8,KIA8) EK9=EK(KE9,KIA9) EK10=EK(KE10,KIA10) EK11=EK(KE11,KIA11) EK12=EK(KE12,KIA12) EK13=EK(KE13,KIA13) EK14=EK(KE14,KIA14) EK15=EK(KE15,KIA15) EK16=EK(KE16,KIA16) EK17=EK(KE17,KIA17) EK18=EK(KE18,KIA18) EK19=EK(KE19,KIA19) EK20=EK(KE20,KIA20) EK21=EK(KE21,KIA21) EK22=EK(KE22,KIA22) EK23=EK(KE23,KIA23) EK24=EK(KE24,KIA24) EK25=EK(KE25,KIA25) EK26=EK(KE26,KIA26) EK27=EK(KE27,KIA27) EK28=EK(KE28,KIA28) EK29=EK(KE29,KIA29) EK30=EK(KE30,KIA30) IKA1=IKA(VA1,KE1,KIA1) IKA2=IKA(VA2,KE2,KIA2) IKA3=IKA(VA3,KE3,KIA3) IKA4=IKA(VA4,KE4,KIA4) IKA5=IKA(VA5,KE5,KIA5) IKA6=IKA(VA6,KE6,KIA6) IKA7=IKA(VA7,KE7,KIA7) IKA8=IKA(VA8,KE8,KIA8) IKA9=IKA(VA9,KE9,KIA9) IKA10=IKA(VA10,KE10,KIA10) IKA11=IKA(VA11,KE11,KIA11) IKA12=IKA(VA12,KE12,KIA12) IKA13=IKA(VA13,KE13,KIA13) IKA14=IKA(VA14,KE14,KIA14) IKA15=IKA(VA15,KE15,KIA15) IKA16=IKA(VA16,KE16,KIA16) IKA17=IKA(VA17,KE17,KIA17) IKA18=IKA(VA18,KE18,KIA18) IKA19=IKA(VA19,KE19,KIA19) IKA20=IKA(VA20,KE20,KIA20) IKA21=IKA(VA21,KE21,KIA21) IKA22=IKA(VA22,KE22,KIA22) IKA23=IKA(VA23,KE23,KIA23) IKA24=IKA(VA24,KE24,KIA24) IKA25=IKA(VA25,KE25,KIA25) IKA26=IKA(VA26,KE26,KIA26) IKA27=IKA(VA27,KE27,KIA27) IKA28=IKA(VA28,KE28,KIA28) IKA29=IKA(VA29,KE29,KIA29) IKA30=IKA(VA30,KE30,KIA30) where ... FRT = R*TEMP/F VOLE = DELTA*VOLN C1 = 10.0*S/(F*VOLN) C2 = 10.0*S/(F*VOLE) C3A = 10*SAS/(F*VOLA) C3E = 10*SAS/(F*VOLE) IPUMP1 = RHON/(((1+2/KE1)^2)*(1+7.7/NAI1)^3) IPUMP2 = RHON/(((1+2/KE2)^2)*(1+7.7/NAI2)^3) IPUMP3 = RHON/(((1+2/KE3)^2)*(1+7.7/NAI3)^3) IPUMP4 = RHON/(((1+2/KE4)^2)*(1+7.7/NAI4)^3) IPUMP5 = RHON/(((1+2/KE5)^2)*(1+7.7/NAI5)^3) IPUMP6 = RHON/(((1+2/KE6)^2)*(1+7.7/NAI6)^3) IPUMP7 = RHON/(((1+2/KE7)^2)*(1+7.7/NAI7)^3) IPUMP8 = RHON/(((1+2/KE8)^2)*(1+7.7/NAI8)^3) IPUMP9 = RHON/(((1+2/KE9)^2)*(1+7.7/NAI9)^3) IPUMP10 = RHON/(((1+2/KE10)^2)*(1+7.7/NAI10)^3) IPUMP11 = RHON/(((1+2/KE11)^2)*(1+7.7/NAI11)^3) IPUMP12 = RHON/(((1+2/KE12)^2)*(1+7.7/NAI12)^3) IPUMP13 = RHON/(((1+2/KE13)^2)*(1+7.7/NAI13)^3) IPUMP14 = RHON/(((1+2/KE14)^2)*(1+7.7/NAI14)^3) IPUMP15 = RHON/(((1+2/KE15)^2)*(1+7.7/NAI15)^3) IPUMP16 = RHON/(((1+2/KE16)^2)*(1+7.7/NAI16)^3) IPUMP17 = RHON/(((1+2/KE17)^2)*(1+7.7/NAI17)^3) IPUMP18 = RHON/(((1+2/KE18)^2)*(1+7.7/NAI18)^3) IPUMP19 = RHON/(((1+2/KE19)^2)*(1+7.7/NAI19)^3) IPUMP20 = RHON/(((1+2/KE20)^2)*(1+7.7/NAI20)^3) IPUMP21 = RHON/(((1+2/KE21)^2)*(1+7.7/NAI21)^3) IPUMP22 = RHON/(((1+2/KE22)^2)*(1+7.7/NAI22)^3) IPUMP23 = RHON/(((1+2/KE23)^2)*(1+7.7/NAI23)^3) IPUMP24 = RHON/(((1+2/KE24)^2)*(1+7.7/NAI24)^3) IPUMP25 = RHON/(((1+2/KE25)^2)*(1+7.7/NAI25)^3) IPUMP26 = RHON/(((1+2/KE26)^2)*(1+7.7/NAI26)^3) IPUMP27 = RHON/(((1+2/KE27)^2)*(1+7.7/NAI27)^3) IPUMP28 = RHON/(((1+2/KE28)^2)*(1+7.7/NAI28)^3) IPUMP29 = RHON/(((1+2/KE29)^2)*(1+7.7/NAI29)^3) IPUMP30 = RHON/(((1+2/KE30)^2)*(1+7.7/NAI30)^3) INAPUMP1 = 3*IPUMP1 INAPUMP2 = 3*IPUMP2 INAPUMP3 = 3*IPUMP3 INAPUMP4 = 3*IPUMP4 INAPUMP5 = 3*IPUMP5 INAPUMP6 = 3*IPUMP6 INAPUMP7 = 3*IPUMP7 INAPUMP8 = 3*IPUMP8 INAPUMP9 = 3*IPUMP9 INAPUMP10 = 3*IPUMP10 INAPUMP11 = 3*IPUMP11 INAPUMP12 = 3*IPUMP12 INAPUMP13 = 3*IPUMP13 INAPUMP14 = 3*IPUMP14 INAPUMP15 = 3*IPUMP15 INAPUMP16 = 3*IPUMP16 INAPUMP17 = 3*IPUMP17 INAPUMP18 = 3*IPUMP18 INAPUMP19 = 3*IPUMP19 INAPUMP20 = 3*IPUMP20 INAPUMP21 = 3*IPUMP21 INAPUMP22 = 3*IPUMP22 INAPUMP23 = 3*IPUMP23 INAPUMP24 = 3*IPUMP24 INAPUMP25 = 3*IPUMP25 INAPUMP26 = 3*IPUMP26 INAPUMP27 = 3*IPUMP27 INAPUMP28 = 3*IPUMP28 INAPUMP29 = 3*IPUMP29 INAPUMP30 = 3*IPUMP30 IKPUMP1 = -2*IPUMP1 IKPUMP2 = -2*IPUMP2 IKPUMP3 = -2*IPUMP3 IKPUMP4 = -2*IPUMP4 IKPUMP5 = -2*IPUMP5 IKPUMP6 = -2*IPUMP6 IKPUMP7 = -2*IPUMP7 IKPUMP8 = -2*IPUMP8 IKPUMP9 = -2*IPUMP9 IKPUMP10 = -2*IPUMP10 IKPUMP11 = -2*IPUMP11 IKPUMP12 = -2*IPUMP12 IKPUMP13 = -2*IPUMP13 IKPUMP14 = -2*IPUMP14 IKPUMP15 = -2*IPUMP15 IKPUMP16 = -2*IPUMP16 IKPUMP17 = -2*IPUMP17 IKPUMP18 = -2*IPUMP18 IKPUMP19 = -2*IPUMP19 IKPUMP20 = -2*IPUMP20 IKPUMP21 = -2*IPUMP21 IKPUMP22 = -2*IPUMP22 IKPUMP23 = -2*IPUMP23 IKPUMP24 = -2*IPUMP24 IKPUMP25 = -2*IPUMP25 IKPUMP26 = -2*IPUMP26 IKPUMP27 = -2*IPUMP27 IKPUMP28 = -2*IPUMP28 IKPUMP29 = -2*IPUMP29 IKPUMP30 = -2*IPUMP30 VK1 = FRT*LN(KE1/KI1) VK2 = FRT*LN(KE2/KI2) VK3 = FRT*LN(KE3/KI3) VK4 = FRT*LN(KE4/KI4) VK5 = FRT*LN(KE5/KI5) VK6 = FRT*LN(KE6/KI6) VK7 = FRT*LN(KE7/KI7) VK8 = FRT*LN(KE8/KI8) VK9 = FRT*LN(KE9/KI9) VK10 = FRT*LN(KE10/KI10) VK11 = FRT*LN(KE11/KI11) VK12 = FRT*LN(KE12/KI12) VK13 = FRT*LN(KE13/KI13) VK14 = FRT*LN(KE14/KI14) VK15 = FRT*LN(KE15/KI15) VK16 = FRT*LN(KE16/KI16) VK17 = FRT*LN(KE17/KI17) VK18 = FRT*LN(KE18/KI18) VK19 = FRT*LN(KE19/KI19) VK20 = FRT*LN(KE20/KI20) VK21 = FRT*LN(KE21/KI21) VK22 = FRT*LN(KE22/KI22) VK23 = FRT*LN(KE23/KI23) VK24 = FRT*LN(KE24/KI24) VK25 = FRT*LN(KE25/KI25) VK26 = FRT*LN(KE26/KI26) VK27 = FRT*LN(KE27/KI27) VK28 = FRT*LN(KE28/KI28) VK29 = FRT*LN(KE29/KI29) VK30 = FRT*LN(KE30/KI30) VNA1 = FRT*LN(NAE1/NAI1) VNA2 = FRT*LN(NAE2/NAI2) VNA3 = FRT*LN(NAE3/NAI3) VNA4 = FRT*LN(NAE4/NAI4) VNA5 = FRT*LN(NAE5/NAI5) VNA6 = FRT*LN(NAE6/NAI6) VNA7 = FRT*LN(NAE7/NAI7) VNA8 = FRT*LN(NAE8/NAI8) VNA9 = FRT*LN(NAE9/NAI9) VNA10 = FRT*LN(NAE10/NAI10) VNA11 = FRT*LN(NAE11/NAI11) VNA12 = FRT*LN(NAE12/NAI12) VNA13 = FRT*LN(NAE13/NAI13) VNA14 = FRT*LN(NAE14/NAI14) VNA15 = FRT*LN(NAE15/NAI15) VNA16 = FRT*LN(NAE16/NAI16) VNA17 = FRT*LN(NAE17/NAI17) VNA18 = FRT*LN(NAE18/NAI18) VNA19 = FRT*LN(NAE19/NAI19) VNA20 = FRT*LN(NAE20/NAI20) VNA21 = FRT*LN(NAE21/NAI21) VNA22 = FRT*LN(NAE22/NAI22) VNA23 = FRT*LN(NAE23/NAI23) VNA24 = FRT*LN(NAE24/NAI24) VNA25 = FRT*LN(NAE25/NAI25) VNA26 = FRT*LN(NAE26/NAI26) VNA27 = FRT*LN(NAE27/NAI27) VNA28 = FRT*LN(NAE28/NAI28) VNA29 = FRT*LN(NAE29/NAI29) VNA30 = FRT*LN(NAE30/NAI30) KDIFF1 = D5*2*(KE2-KE1) KDIFF2 = D5*(KE3+KE1-2*KE2) KDIFF3 = D5*(KE4+KE2-2*KE3) KDIFF4 = D5*(KE5+KE3-2*KE4) KDIFF5 = D5*(KE6+KE4-2*KE5) KDIFF6 = D5*(KE7+KE5-2*KE6) KDIFF7 = D5*(KE8+KE6-2*KE7) KDIFF8 = D5*(KE9+KE7-2*KE8) KDIFF9 = D5*(KE10+KE8-2*KE9) KDIFF10 = D5*(KE11+KE9-2*KE10) KDIFF11 = D5*(KE12+KE10-2*KE11) KDIFF12 = D5*(KE13+KE11-2*KE12) KDIFF13 = D5*(KE14+KE12-2*KE13) KDIFF14 = D5*(KE15+KE13-2*KE14) KDIFF15 = D5*(KE16+KE14-2*KE15) KDIFF16 = D5*(KE17+KE15-2*KE16) KDIFF17 = D5*(KE18+KE16-2*KE17) KDIFF18 = D5*(KE19+KE17-2*KE18) KDIFF19 = D5*(KE20+KE18-2*KE19) KDIFF20 = D5*(KE21+KE19-2*KE20) KDIFF21 = D5*(KE22+KE20-2*KE21) KDIFF22 = D5*(KE23+KE21-2*KE22) KDIFF23 = D5*(KE24+KE22-2*KE23) KDIFF24 = D5*(KE25+KE23-2*KE24) KDIFF25 = D5*(KE26+KE24-2*KE25) KDIFF26 = D5*(KE27+KE25-2*KE26) KDIFF27 = D5*(KE28+KE26-2*KE27) KDIFF28 = D5*(KE29+KE27-2*KE28) KDIFF29 = D5*(KE30+KE28-2*KE29) KDIFF30 = D5*2*(KE29-KE30) NADIFF1 = D8*2*(NAE2-NAE1) NADIFF2 = D8*(NAE3+NAE1-2*NAE2) NADIFF3 = D8*(NAE4+NAE2-2*NAE3) NADIFF4 = D8*(NAE5+NAE3-2*NAE4) NADIFF5 = D8*(NAE6+NAE4-2*NAE5) NADIFF6 = D8*(NAE7+NAE5-2*NAE6) NADIFF7 = D8*(NAE8+NAE6-2*NAE7) NADIFF8 = D8*(NAE9+NAE7-2*NAE8) NADIFF9 = D8*(NAE10+NAE8-2*NAE9) NADIFF10 = D8*(NAE11+NAE9-2*NAE10) NADIFF11 = D8*(NAE12+NAE10-2*NAE11) NADIFF12 = D8*(NAE13+NAE11-2*NAE12) NADIFF13 = D8*(NAE14+NAE12-2*NAE13) NADIFF14 = D8*(NAE15+NAE13-2*NAE14) NADIFF15 = D8*(NAE16+NAE14-2*NAE15) NADIFF16 = D8*(NAE17+NAE15-2*NAE16) NADIFF17 = D8*(NAE18+NAE16-2*NAE17) NADIFF18 = D8*(NAE19+NAE17-2*NAE18) NADIFF19 = D8*(NAE20+NAE18-2*NAE19) NADIFF20 = D8*(NAE21+NAE19-2*NAE20) NADIFF21 = D8*(NAE22+NAE20-2*NAE21) NADIFF22 = D8*(NAE23+NAE21-2*NAE22) NADIFF23 = D8*(NAE24+NAE22-2*NAE23) NADIFF24 = D8*(NAE25+NAE23-2*NAE24) NADIFF25 = D8*(NAE26+NAE24-2*NAE25) NADIFF26 = D8*(NAE27+NAE25-2*NAE26) NADIFF27 = D8*(NAE28+NAE26-2*NAE27) NADIFF28 = D8*(NAE29+NAE27-2*NAE28) NADIFF29 = D8*(NAE30+NAE28-2*NAE29) NADIFF30 = D8*2*(NAE29-NAE30) PKGAP = SIGMAGAP*PKA PNAGAP = .8*PKGAP IKGA1 = KG(VA1,VA2,KIA1,KIA2)+KG(VA1,VA3,KIA1,KIA3)+KG(VA1,VA4,KIA1,KIA4)+KG(VA1,VA5,KIA1,KIA5)+KG(VA1,VA6,KIA1,KIA6) IKGA2 = KG(VA2,VA3,KIA2,KIA3)+KG(VA2,VA4,KIA2,KIA4)+KG(VA2,VA5,KIA2,KIA5)+KG(VA2,VA6,KIA2,KIA6)+KG(VA2,VA7,KIA2,KIA7) IKGA3 = KG(VA3,VA4,KIA3,KIA4)+KG(VA3,VA5,KIA3,KIA5)+KG(VA3,VA6,KIA3,KIA6)+KG(VA3,VA7,KIA3,KIA7)+KG(VA3,VA8,KIA3,KIA8) IKGA4 = KG(VA4,VA5,KIA4,KIA5)+KG(VA4,VA6,KIA4,KIA6)+KG(VA4,VA7,KIA4,KIA7)+KG(VA4,VA8,KIA4,KIA8)+KG(VA4,VA9,KIA4,KIA9) IKGA5 = KG(VA5,VA6,KIA5,KIA6)+KG(VA5,VA7,KIA5,KIA7)+KG(VA5,VA8,KIA5,KIA8)+KG(VA5,VA9,KIA5,KIA9)+KG(VA5,VA10,KIA5,KIA10) IKGA6 = KG(VA6,VA7,KIA6,KIA7)+KG(VA6,VA8,KIA6,KIA8)+KG(VA6,VA9,KIA6,KIA9)+KG(VA6,VA10,KIA6,KIA10)+KG(VA6,VA11,KIA6,KIA11) IKGA7 = KG(VA7,VA8,KIA7,KIA8)+KG(VA7,VA9,KIA7,KIA9)+KG(VA7,VA10,KIA7,KIA10)+KG(VA7,VA11,KIA7,KIA11)+KG(VA7,VA12,KIA7,KIA12) IKGA8 = KG(VA8,VA9,KIA8,KIA9)+KG(VA8,VA10,KIA8,KIA10)+KG(VA8,VA11,KIA8,KIA11)+KG(VA8,VA12,KIA8,KIA12)+KG(VA8,VA13,KIA8,KIA13) IKGA9 = KG(VA9,VA10,KIA9,KIA10)+KG(VA9,VA11,KIA9,KIA11)+KG(VA9,VA12,KIA9,KIA12)+KG(VA9,VA13,KIA9,KIA13)+KG(VA9,VA14,KIA9,KIA14) IKGA10 = KG(VA10,VA11,KIA10,KIA11)+KG(VA10,VA12,KIA10,KIA12)+KG(VA10,VA13,KIA10,KIA13)+KG(VA10,VA14,KIA10,KIA14)+KG(VA10,VA15,KIA10,KIA15) IKGA11 = KG(VA11,VA12,KIA11,KIA12)+KG(VA11,VA13,KIA11,KIA13)+KG(VA11,VA14,KIA11,KIA14)+KG(VA11,VA15,KIA11,KIA15)+KG(VA11,VA16,KIA11,KIA16) IKGA12 = KG(VA12,VA13,KIA12,KIA13)+KG(VA12,VA14,KIA12,KIA14)+KG(VA12,VA15,KIA12,KIA15)+KG(VA12,VA16,KIA12,KIA16)+KG(VA12,VA17,KIA12,KIA17) IKGA13 = KG(VA13,VA14,KIA13,KIA14)+KG(VA13,VA15,KIA13,KIA15)+KG(VA13,VA16,KIA13,KIA16)+KG(VA13,VA17,KIA13,KIA17)+KG(VA13,VA18,KIA13,KIA18) IKGA14 = KG(VA14,VA15,KIA14,KIA15)+KG(VA14,VA16,KIA14,KIA16)+KG(VA14,VA17,KIA14,KIA17)+KG(VA14,VA18,KIA14,KIA18)+KG(VA14,VA19,KIA14,KIA19) IKGA15 = KG(VA15,VA16,KIA15,KIA16)+KG(VA15,VA17,KIA15,KIA17)+KG(VA15,VA18,KIA15,KIA18)+KG(VA15,VA19,KIA15,KIA19)+KG(VA15,VA20,KIA15,KIA20) IKGA16 = KG(VA16,VA17,KIA16,KIA17)+KG(VA16,VA18,KIA16,KIA18)+KG(VA16,VA19,KIA16,KIA19)+KG(VA16,VA20,KIA16,KIA20)+KG(VA16,VA21,KIA16,KIA21) IKGA17 = KG(VA17,VA18,KIA17,KIA18)+KG(VA17,VA19,KIA17,KIA19)+KG(VA17,VA20,KIA17,KIA20)+KG(VA17,VA21,KIA17,KIA21)+KG(VA17,VA22,KIA17,KIA22) IKGA18 = KG(VA18,VA19,KIA18,KIA19)+KG(VA18,VA20,KIA18,KIA20)+KG(VA18,VA21,KIA18,KIA21)+KG(VA18,VA22,KIA18,KIA22)+KG(VA18,VA23,KIA18,KIA23) IKGA19 = KG(VA19,VA20,KIA19,KIA20)+KG(VA19,VA21,KIA19,KIA21)+KG(VA19,VA22,KIA19,KIA22)+KG(VA19,VA23,KIA19,KIA23)+KG(VA19,VA24,KIA19,KIA24) IKGA20 = KG(VA20,VA21,KIA20,KIA21)+KG(VA20,VA22,KIA20,KIA22)+KG(VA20,VA23,KIA20,KIA23)+KG(VA20,VA24,KIA20,KIA24)+KG(VA20,VA25,KIA20,KIA25) IKGA21 = KG(VA21,VA22,KIA21,KIA22)+KG(VA21,VA23,KIA21,KIA23)+KG(VA21,VA24,KIA21,KIA24)+KG(VA21,VA25,KIA21,KIA25)+KG(VA21,VA26,KIA21,KIA26) IKGA22 = KG(VA22,VA23,KIA22,KIA23)+KG(VA22,VA24,KIA22,KIA24)+KG(VA22,VA25,KIA22,KIA25)+KG(VA22,VA26,KIA22,KIA26)+KG(VA22,VA27,KIA22,KIA27) IKGA23 = KG(VA23,VA24,KIA23,KIA24)+KG(VA23,VA25,KIA23,KIA25)+KG(VA23,VA26,KIA23,KIA26)+KG(VA23,VA27,KIA23,KIA27)+KG(VA23,VA28,KIA23,KIA28) IKGA24 = KG(VA24,VA25,KIA24,KIA25)+KG(VA24,VA26,KIA24,KIA26)+KG(VA24,VA27,KIA24,KIA27)+KG(VA24,VA28,KIA24,KIA28)+KG(VA24,VA29,KIA24,KIA29) IKGA25 = KG(VA25,VA26,KIA25,KIA26)+KG(VA25,VA27,KIA25,KIA27)+KG(VA25,VA28,KIA25,KIA28)+KG(VA25,VA29,KIA25,KIA29)+KG(VA25,VA30,KIA25,KIA30) IKGA26 = KG(VA26,VA27,KIA26,KIA27)+KG(VA26,VA28,KIA26,KIA28)+KG(VA26,VA29,KIA26,KIA29)+KG(VA26,VA30,KIA26,KIA30) IKGA27 = KG(VA27,VA28,KIA27,KIA28)+KG(VA27,VA29,KIA27,KIA29)+KG(VA27,VA30,KIA27,KIA30) IKGA28 = KG(VA28,VA29,KIA28,KIA29)+KG(VA28,VA30,KIA28,KIA30) IKGA29 = KG(VA29,VA30,KIA29,KIA30) IKGA30 = 0 IKGB6 = KG(VA6,VA5,KIA6,KIA5)+KG(VA6,VA4,KIA6,KIA4)+KG(VA6,VA3,KIA6,KIA3)+KG(VA6,VA2,KIA6,KIA2)+KG(VA6,VA1,KIA6,KIA1) IKGB7 = KG(VA7,VA6,KIA7,KIA6)+KG(VA7,VA5,KIA7,KIA5)+KG(VA7,VA4,KIA7,KIA4)+KG(VA7,VA3,KIA7,KIA3)+KG(VA7,VA2,KIA7,KIA2) IKGB8 = KG(VA8,VA7,KIA8,KIA7)+KG(VA8,VA6,KIA8,KIA6)+KG(VA8,VA5,KIA8,KIA5)+KG(VA8,VA4,KIA8,KIA4)+KG(VA8,VA3,KIA8,KIA3) IKGB9 = KG(VA9,VA8,KIA9,KIA8)+KG(VA9,VA7,KIA9,KIA7)+KG(VA9,VA6,KIA9,KIA6)+KG(VA9,VA5,KIA9,KIA5)+KG(VA9,VA4,KIA9,KIA4) IKGB10 = KG(VA10,VA9,KIA10,KIA9)+KG(VA10,VA8,KIA10,KIA8)+KG(VA10,VA7,KIA10,KIA7)+KG(VA10,VA6,KIA10,KIA6)+KG(VA10,VA5,KIA10,KIA5) IKGB11 = KG(VA11,VA10,KIA11,KIA10)+KG(VA11,VA9,KIA11,KIA9)+KG(VA11,VA8,KIA11,KIA8)+KG(VA11,VA7,KIA11,KIA7)+KG(VA11,VA6,KIA11,KIA6) IKGB12 = KG(VA12,VA11,KIA12,KIA11)+KG(VA12,VA10,KIA12,KIA10)+KG(VA12,VA9,KIA12,KIA9)+KG(VA12,VA8,KIA12,KIA8)+KG(VA12,VA7,KIA12,KIA7) IKGB13 = KG(VA13,VA12,KIA13,KIA12)+KG(VA13,VA11,KIA13,KIA11)+KG(VA13,VA10,KIA13,KIA10)+KG(VA13,VA9,KIA13,KIA9)+KG(VA13,VA8,KIA13,KIA8) IKGB14 = KG(VA14,VA13,KIA14,KIA13)+KG(VA14,VA12,KIA14,KIA12)+KG(VA14,VA11,KIA14,KIA11)+KG(VA14,VA10,KIA14,KIA10)+KG(VA14,VA9,KIA14,KIA9) IKGB15 = KG(VA15,VA14,KIA15,KIA14)+KG(VA15,VA13,KIA15,KIA13)+KG(VA15,VA12,KIA15,KIA12)+KG(VA15,VA11,KIA15,KIA11)+KG(VA15,VA10,KIA15,KIA10) IKGB16 = KG(VA16,VA15,KIA16,KIA15)+KG(VA16,VA14,KIA16,KIA14)+KG(VA16,VA13,KIA16,KIA13)+KG(VA16,VA12,KIA16,KIA12)+KG(VA16,VA11,KIA16,KIA11) IKGB17 = KG(VA17,VA16,KIA17,KIA16)+KG(VA17,VA15,KIA17,KIA15)+KG(VA17,VA14,KIA17,KIA14)+KG(VA17,VA13,KIA17,KIA13)+KG(VA17,VA12,KIA17,KIA12) IKGB18 = KG(VA18,VA17,KIA18,KIA17)+KG(VA18,VA16,KIA18,KIA16)+KG(VA18,VA15,KIA18,KIA15)+KG(VA18,VA14,KIA18,KIA14)+KG(VA18,VA13,KIA18,KIA13) IKGB19 = KG(VA19,VA18,KIA19,KIA18)+KG(VA19,VA17,KIA19,KIA17)+KG(VA19,VA16,KIA19,KIA16)+KG(VA19,VA15,KIA19,KIA15)+KG(VA19,VA14,KIA19,KIA14) IKGB20 = KG(VA20,VA19,KIA20,KIA19)+KG(VA20,VA18,KIA20,KIA18)+KG(VA20,VA17,KIA20,KIA17)+KG(VA20,VA16,KIA20,KIA16)+KG(VA20,VA15,KIA20,KIA15) IKGB21 = KG(VA21,VA20,KIA21,KIA20)+KG(VA21,VA19,KIA21,KIA19)+KG(VA21,VA18,KIA21,KIA18)+KG(VA21,VA17,KIA21,KIA17)+KG(VA21,VA16,KIA21,KIA16) IKGB22 = KG(VA22,VA21,KIA22,KIA21)+KG(VA22,VA20,KIA22,KIA20)+KG(VA22,VA19,KIA22,KIA19)+KG(VA22,VA18,KIA22,KIA18)+KG(VA22,VA17,KIA22,KIA17) IKGB23 = KG(VA23,VA22,KIA23,KIA22)+KG(VA23,VA21,KIA23,KIA21)+KG(VA23,VA20,KIA23,KIA20)+KG(VA23,VA19,KIA23,KIA19)+KG(VA23,VA18,KIA23,KIA18) IKGB24 = KG(VA24,VA23,KIA24,KIA23)+KG(VA24,VA22,KIA24,KIA22)+KG(VA24,VA21,KIA24,KIA21)+KG(VA24,VA20,KIA24,KIA20)+KG(VA24,VA19,KIA24,KIA19) IKGB25 = KG(VA25,VA24,KIA25,KIA24)+KG(VA25,VA23,KIA25,KIA23)+KG(VA25,VA22,KIA25,KIA22)+KG(VA25,VA21,KIA25,KIA21)+KG(VA25,VA20,KIA25,KIA20) IKGB26 = KG(VA26,VA25,KIA26,KIA25)+KG(VA26,VA24,KIA26,KIA24)+KG(VA26,VA23,KIA26,KIA23)+KG(VA26,VA22,KIA26,KIA22)+KG(VA26,VA21,KIA26,KIA21) IKGB27 = KG(VA27,VA26,KIA27,KIA26)+KG(VA27,VA25,KIA27,KIA25)+KG(VA27,VA24,KIA27,KIA24)+KG(VA27,VA23,KIA27,KIA23)+KG(VA27,VA22,KIA27,KIA22) IKGB28 = KG(VA28,VA27,KIA28,KIA27)+KG(VA28,VA26,KIA28,KIA26)+KG(VA28,VA25,KIA28,KIA25)+KG(VA28,VA24,KIA28,KIA24)+KG(VA28,VA23,KIA28,KIA23) IKGB29 = KG(VA29,VA28,KIA29,KIA28)+KG(VA29,VA27,KIA29,KIA27)+KG(VA29,VA26,KIA29,KIA26)+KG(VA29,VA25,KIA29,KIA25)+KG(VA29,VA24,KIA29,KIA24) IKGB30 = KG(VA30,VA29,KIA30,KIA29)+KG(VA30,VA28,KIA30,KIA28)+KG(VA30,VA27,KIA30,KIA27)+KG(VA30,VA26,KIA30,KIA26)+KG(VA30,VA25,KIA30,KIA25) IKGB5 = KG(VA5,VA4,KIA5,KIA4)+KG(VA5,VA3,KIA5,KIA3)+KG(VA5,VA2,KIA5,KIA2)+KG(VA5,VA1,KIA5,KIA1) IKGB4 = KG(VA4,VA3,KIA4,KIA3)+KG(VA4,VA2,KIA4,KIA2)+KG(VA4,VA1,KIA4,KIA1) IKGB3 = KG(VA3,VA2,KIA3,KIA2)+KG(VA3,VA1,KIA3,KIA1) IKGB2 = KG(VA2,VA1,KIA2,KIA1) IKGB1 = 0 IKGAP1 = IKGA1+IKGB1 IKGAP2 = IKGA2+IKGB2 IKGAP3 = IKGA3+IKGB3 IKGAP4 = IKGA4+IKGB4 IKGAP5 = IKGA5+IKGB5 IKGAP6 = IKGA6+IKGB6 IKGAP7 = IKGA7+IKGB7 IKGAP8 = IKGA8+IKGB8 IKGAP9 = IKGA9+IKGB9 IKGAP10 = IKGA10+IKGB10 IKGAP11 = IKGA11+IKGB11 IKGAP12 = IKGA12+IKGB12 IKGAP13 = IKGA13+IKGB13 IKGAP14 = IKGA14+IKGB14 IKGAP15 = IKGA15+IKGB15 IKGAP16 = IKGA16+IKGB16 IKGAP17 = IKGA17+IKGB17 IKGAP18 = IKGA18+IKGB18 IKGAP19 = IKGA19+IKGB19 IKGAP20 = IKGA20+IKGB20 IKGAP21 = IKGA21+IKGB21 IKGAP22 = IKGA22+IKGB22 IKGAP23 = IKGA23+IKGB23 IKGAP24 = IKGA24+IKGB24 IKGAP25 = IKGA25+IKGB25 IKGAP26 = IKGA26+IKGB26 IKGAP27 = IKGA27+IKGB27 IKGAP28 = IKGA28+IKGB28 IKGAP29 = IKGA29+IKGB29 IKGAP30 = IKGA30+IKGB30 INAGA1 = NAG(VA1,VA2,NAIA1,NAIA2)+NAG(VA1,VA3,NAIA1,NAIA3)+NAG(VA1,VA4,NAIA1,NAIA4)+NAG(VA1,VA5,NAIA1,NAIA5)+NAG(VA1,VA6,NAIA1,NAIA6) INAGA2 = NAG(VA2,VA3,NAIA2,NAIA3)+NAG(VA2,VA4,NAIA2,NAIA4)+NAG(VA2,VA5,NAIA2,NAIA5)+NAG(VA2,VA6,NAIA2,NAIA6)+NAG(VA2,VA7,NAIA2,NAIA7) INAGA3 = NAG(VA3,VA4,NAIA3,NAIA4)+NAG(VA3,VA5,NAIA3,NAIA5)+NAG(VA3,VA6,NAIA3,NAIA6)+NAG(VA3,VA7,NAIA3,NAIA7)+NAG(VA3,VA8,NAIA3,NAIA8) INAGA4 = NAG(VA4,VA5,NAIA4,NAIA5)+NAG(VA4,VA6,NAIA4,NAIA6)+NAG(VA4,VA7,NAIA4,NAIA7)+NAG(VA4,VA8,NAIA4,NAIA8)+NAG(VA4,VA9,NAIA4,NAIA9) INAGA5 = NAG(VA5,VA6,NAIA5,NAIA6)+NAG(VA5,VA7,NAIA5,NAIA7)+NAG(VA5,VA8,NAIA5,NAIA8)+NAG(VA5,VA9,NAIA5,NAIA9)+NAG(VA5,VA10,NAIA5,NAIA10) INAGA6 = NAG(VA6,VA7,NAIA6,NAIA7)+NAG(VA6,VA8,NAIA6,NAIA8)+NAG(VA6,VA9,NAIA6,NAIA9)+NAG(VA6,VA10,NAIA6,NAIA10)+NAG(VA6,VA11,NAIA6,NAIA11) INAGA7 = NAG(VA7,VA8,NAIA7,NAIA8)+NAG(VA7,VA9,NAIA7,NAIA9)+NAG(VA7,VA10,NAIA7,NAIA10)+NAG(VA7,VA11,NAIA7,NAIA11)+NAG(VA7,VA12,NAIA7,NAIA12) INAGA8 = NAG(VA8,VA9,NAIA8,NAIA9)+NAG(VA8,VA10,NAIA8,NAIA10)+NAG(VA8,VA11,NAIA8,NAIA11)+NAG(VA8,VA12,NAIA8,NAIA12)+NAG(VA8,VA13,NAIA8,NAIA13) INAGA9 = NAG(VA9,VA10,NAIA9,NAIA10)+NAG(VA9,VA11,NAIA9,NAIA11)+NAG(VA9,VA12,NAIA9,NAIA12)+NAG(VA9,VA13,NAIA9,NAIA13)+NAG(VA9,VA14,NAIA9,NAIA14) INAGA10 = NAG(VA10,VA11,NAIA10,NAIA11)+NAG(VA10,VA12,NAIA10,NAIA12)+NAG(VA10,VA13,NAIA10,NAIA13)+NAG(VA10,VA14,NAIA10,NAIA14)+NAG(VA10,VA15,NAIA10,NAIA15) INAGA11 = NAG(VA11,VA12,NAIA11,NAIA12)+NAG(VA11,VA13,NAIA11,NAIA13)+NAG(VA11,VA14,NAIA11,NAIA14)+NAG(VA11,VA15,NAIA11,NAIA15)+NAG(VA11,VA16,NAIA11,NAIA16) INAGA12 = NAG(VA12,VA13,NAIA12,NAIA13)+NAG(VA12,VA14,NAIA12,NAIA14)+NAG(VA12,VA15,NAIA12,NAIA15)+NAG(VA12,VA16,NAIA12,NAIA16)+NAG(VA12,VA17,NAIA12,NAIA17) INAGA13 = NAG(VA13,VA14,NAIA13,NAIA14)+NAG(VA13,VA15,NAIA13,NAIA15)+NAG(VA13,VA16,NAIA13,NAIA16)+NAG(VA13,VA17,NAIA13,NAIA17)+NAG(VA13,VA18,NAIA13,NAIA18) INAGA14 = NAG(VA14,VA15,NAIA14,NAIA15)+NAG(VA14,VA16,NAIA14,NAIA16)+NAG(VA14,VA17,NAIA14,NAIA17)+NAG(VA14,VA18,NAIA14,NAIA18)+NAG(VA14,VA19,NAIA14,NAIA19) INAGA15 = NAG(VA15,VA16,NAIA15,NAIA16)+NAG(VA15,VA17,NAIA15,NAIA17)+NAG(VA15,VA18,NAIA15,NAIA18)+NAG(VA15,VA19,NAIA15,NAIA19)+NAG(VA15,VA20,NAIA15,NAIA20) INAGA16 = NAG(VA16,VA17,NAIA16,NAIA17)+NAG(VA16,VA18,NAIA16,NAIA18)+NAG(VA16,VA19,NAIA16,NAIA19)+NAG(VA16,VA20,NAIA16,NAIA20)+NAG(VA16,VA21,NAIA16,NAIA21) INAGA17 = NAG(VA17,VA18,NAIA17,NAIA18)+NAG(VA17,VA19,NAIA17,NAIA19)+NAG(VA17,VA20,NAIA17,NAIA20)+NAG(VA17,VA21,NAIA17,NAIA21)+NAG(VA17,VA22,NAIA17,NAIA22) INAGA18 = NAG(VA18,VA19,NAIA18,NAIA19)+NAG(VA18,VA20,NAIA18,NAIA20)+NAG(VA18,VA21,NAIA18,NAIA21)+NAG(VA18,VA22,NAIA18,NAIA22)+NAG(VA18,VA23,NAIA18,NAIA23) INAGA19 = NAG(VA19,VA20,NAIA19,NAIA20)+NAG(VA19,VA21,NAIA19,NAIA21)+NAG(VA19,VA22,NAIA19,NAIA22)+NAG(VA19,VA23,NAIA19,NAIA23)+NAG(VA19,VA24,NAIA19,NAIA24) INAGA20 = NAG(VA20,VA21,NAIA20,NAIA21)+NAG(VA20,VA22,NAIA20,NAIA22)+NAG(VA20,VA23,NAIA20,NAIA23)+NAG(VA20,VA24,NAIA20,NAIA24)+NAG(VA20,VA25,NAIA20,NAIA25) INAGA21 = NAG(VA21,VA22,NAIA21,NAIA22)+NAG(VA21,VA23,NAIA21,NAIA23)+NAG(VA21,VA24,NAIA21,NAIA24)+NAG(VA21,VA25,NAIA21,NAIA25)+NAG(VA21,VA26,NAIA21,NAIA26) INAGA22 = NAG(VA22,VA23,NAIA22,NAIA23)+NAG(VA22,VA24,NAIA22,NAIA24)+NAG(VA22,VA25,NAIA22,NAIA25)+NAG(VA22,VA26,NAIA22,NAIA26)+NAG(VA22,VA27,NAIA22,NAIA27) INAGA23 = NAG(VA23,VA24,NAIA23,NAIA24)+NAG(VA23,VA25,NAIA23,NAIA25)+NAG(VA23,VA26,NAIA23,NAIA26)+NAG(VA23,VA27,NAIA23,NAIA27)+NAG(VA23,VA28,NAIA23,NAIA28) INAGA24 = NAG(VA24,VA25,NAIA24,NAIA25)+NAG(VA24,VA26,NAIA24,NAIA26)+NAG(VA24,VA27,NAIA24,NAIA27)+NAG(VA24,VA28,NAIA24,NAIA28)+NAG(VA24,VA29,NAIA24,NAIA29) INAGA25 = NAG(VA25,VA26,NAIA25,NAIA26)+NAG(VA25,VA27,NAIA25,NAIA27)+NAG(VA25,VA28,NAIA25,NAIA28)+NAG(VA25,VA29,NAIA25,NAIA29)+NAG(VA25,VA30,NAIA25,NAIA30) INAGA26 = NAG(VA26,VA27,NAIA26,NAIA27)+NAG(VA26,VA28,NAIA26,NAIA28)+NAG(VA26,VA29,NAIA26,NAIA29)+NAG(VA26,VA30,NAIA26,NAIA30) INAGA27 = NAG(VA27,VA28,NAIA27,NAIA28)+NAG(VA27,VA29,NAIA27,NAIA29)+NAG(VA27,VA30,NAIA27,NAIA30) INAGA28 = NAG(VA28,VA29,NAIA28,NAIA29)+NAG(VA28,VA30,NAIA28,NAIA30) INAGA29 = NAG(VA29,VA30,NAIA29,NAIA30) INAGA30 = 0 INAGB6 = NAG(VA6,VA5,NAIA6,NAIA5)+NAG(VA6,VA4,NAIA6,NAIA4)+NAG(VA6,VA3,NAIA6,NAIA3)+NAG(VA6,VA2,NAIA6,NAIA2)+NAG(VA6,VA1,NAIA6,NAIA1) INAGB7 = NAG(VA7,VA6,NAIA7,NAIA6)+NAG(VA7,VA5,NAIA7,NAIA5)+NAG(VA7,VA4,NAIA7,NAIA4)+NAG(VA7,VA3,NAIA7,NAIA3)+NAG(VA7,VA2,NAIA7,NAIA2) INAGB8 = NAG(VA8,VA7,NAIA8,NAIA7)+NAG(VA8,VA6,NAIA8,NAIA6)+NAG(VA8,VA5,NAIA8,NAIA5)+NAG(VA8,VA4,NAIA8,NAIA4)+NAG(VA8,VA3,NAIA8,NAIA3) INAGB9 = NAG(VA9,VA8,NAIA9,NAIA8)+NAG(VA9,VA7,NAIA9,NAIA7)+NAG(VA9,VA6,NAIA9,NAIA6)+NAG(VA9,VA5,NAIA9,NAIA5)+NAG(VA9,VA4,NAIA9,NAIA4) INAGB10 = NAG(VA10,VA9,NAIA10,NAIA9)+NAG(VA10,VA8,NAIA10,NAIA8)+NAG(VA10,VA7,NAIA10,NAIA7)+NAG(VA10,VA6,NAIA10,NAIA6)+NAG(VA10,VA5,NAIA10,NAIA5) INAGB11 = NAG(VA11,VA10,NAIA11,NAIA10)+NAG(VA11,VA9,NAIA11,NAIA9)+NAG(VA11,VA8,NAIA11,NAIA8)+NAG(VA11,VA7,NAIA11,NAIA7)+NAG(VA11,VA6,NAIA11,NAIA6) INAGB12 = NAG(VA12,VA11,NAIA12,NAIA11)+NAG(VA12,VA10,NAIA12,NAIA10)+NAG(VA12,VA9,NAIA12,NAIA9)+NAG(VA12,VA8,NAIA12,NAIA8)+NAG(VA12,VA7,NAIA12,NAIA7) INAGB13 = NAG(VA13,VA12,NAIA13,NAIA12)+NAG(VA13,VA11,NAIA13,NAIA11)+NAG(VA13,VA10,NAIA13,NAIA10)+NAG(VA13,VA9,NAIA13,NAIA9)+NAG(VA13,VA8,NAIA13,NAIA8) INAGB14 = NAG(VA14,VA13,NAIA14,NAIA13)+NAG(VA14,VA12,NAIA14,NAIA12)+NAG(VA14,VA11,NAIA14,NAIA11)+NAG(VA14,VA10,NAIA14,NAIA10)+NAG(VA14,VA9,NAIA14,NAIA9) INAGB15 = NAG(VA15,VA14,NAIA15,NAIA14)+NAG(VA15,VA13,NAIA15,NAIA13)+NAG(VA15,VA12,NAIA15,NAIA12)+NAG(VA15,VA11,NAIA15,NAIA11)+NAG(VA15,VA10,NAIA15,NAIA10) INAGB16 = NAG(VA16,VA15,NAIA16,NAIA15)+NAG(VA16,VA14,NAIA16,NAIA14)+NAG(VA16,VA13,NAIA16,NAIA13)+NAG(VA16,VA12,NAIA16,NAIA12)+NAG(VA16,VA11,NAIA16,NAIA11) INAGB17 = NAG(VA17,VA16,NAIA17,NAIA16)+NAG(VA17,VA15,NAIA17,NAIA15)+NAG(VA17,VA14,NAIA17,NAIA14)+NAG(VA17,VA13,NAIA17,NAIA13)+NAG(VA17,VA12,NAIA17,NAIA12) INAGB18 = NAG(VA18,VA17,NAIA18,NAIA17)+NAG(VA18,VA16,NAIA18,NAIA16)+NAG(VA18,VA15,NAIA18,NAIA15)+NAG(VA18,VA14,NAIA18,NAIA14)+NAG(VA18,VA13,NAIA18,NAIA13) INAGB19 = NAG(VA19,VA18,NAIA19,NAIA18)+NAG(VA19,VA17,NAIA19,NAIA17)+NAG(VA19,VA16,NAIA19,NAIA16)+NAG(VA19,VA15,NAIA19,NAIA15)+NAG(VA19,VA14,NAIA19,NAIA14) INAGB20 = NAG(VA20,VA19,NAIA20,NAIA19)+NAG(VA20,VA18,NAIA20,NAIA18)+NAG(VA20,VA17,NAIA20,NAIA17)+NAG(VA20,VA16,NAIA20,NAIA16)+NAG(VA20,VA15,NAIA20,NAIA15) INAGB21 = NAG(VA21,VA20,NAIA21,NAIA20)+NAG(VA21,VA19,NAIA21,NAIA19)+NAG(VA21,VA18,NAIA21,NAIA18)+NAG(VA21,VA17,NAIA21,NAIA17)+NAG(VA21,VA16,NAIA21,NAIA16) INAGB22 = NAG(VA22,VA21,NAIA22,NAIA21)+NAG(VA22,VA20,NAIA22,NAIA20)+NAG(VA22,VA19,NAIA22,NAIA19)+NAG(VA22,VA18,NAIA22,NAIA18)+NAG(VA22,VA17,NAIA22,NAIA17) INAGB23 = NAG(VA23,VA22,NAIA23,NAIA22)+NAG(VA23,VA21,NAIA23,NAIA21)+NAG(VA23,VA20,NAIA23,NAIA20)+NAG(VA23,VA19,NAIA23,NAIA19)+NAG(VA23,VA18,NAIA23,NAIA18) INAGB24 = NAG(VA24,VA23,NAIA24,NAIA23)+NAG(VA24,VA22,NAIA24,NAIA22)+NAG(VA24,VA21,NAIA24,NAIA21)+NAG(VA24,VA20,NAIA24,NAIA20)+NAG(VA24,VA19,NAIA24,NAIA19) INAGB25 = NAG(VA25,VA24,NAIA25,NAIA24)+NAG(VA25,VA23,NAIA25,NAIA23)+NAG(VA25,VA22,NAIA25,NAIA22)+NAG(VA25,VA21,NAIA25,NAIA21)+NAG(VA25,VA20,NAIA25,NAIA20) INAGB26 = NAG(VA26,VA25,NAIA26,NAIA25)+NAG(VA26,VA24,NAIA26,NAIA24)+NAG(VA26,VA23,NAIA26,NAIA23)+NAG(VA26,VA22,NAIA26,NAIA22)+NAG(VA26,VA21,NAIA26,NAIA21) INAGB27 = NAG(VA27,VA26,NAIA27,NAIA26)+NAG(VA27,VA25,NAIA27,NAIA25)+NAG(VA27,VA24,NAIA27,NAIA24)+NAG(VA27,VA23,NAIA27,NAIA23)+NAG(VA27,VA22,NAIA27,NAIA22) INAGB28 = NAG(VA28,VA27,NAIA28,NAIA27)+NAG(VA28,VA26,NAIA28,NAIA26)+NAG(VA28,VA25,NAIA28,NAIA25)+NAG(VA28,VA24,NAIA28,NAIA24)+NAG(VA28,VA23,NAIA28,NAIA23) INAGB29 = NAG(VA29,VA28,NAIA29,NAIA28)+NAG(VA29,VA27,NAIA29,NAIA27)+NAG(VA29,VA26,NAIA29,NAIA26)+NAG(VA29,VA25,NAIA29,NAIA25)+NAG(VA29,VA24,NAIA29,NAIA24) INAGB30 = NAG(VA30,VA29,NAIA30,NAIA29)+NAG(VA30,VA28,NAIA30,NAIA28)+NAG(VA30,VA27,NAIA30,NAIA27)+NAG(VA30,VA26,NAIA30,NAIA26)+NAG(VA30,VA25,NAIA30,NAIA25) INAGB5 = NAG(VA5,VA4,NAIA5,NAIA4)+NAG(VA5,VA3,NAIA5,NAIA3)+NAG(VA5,VA2,NAIA5,NAIA2)+NAG(VA5,VA1,NAIA5,NAIA1) INAGB4 = NAG(VA4,VA3,NAIA4,NAIA3)+NAG(VA4,VA2,NAIA4,NAIA2)+NAG(VA4,VA1,NAIA4,NAIA1) INAGB3 = NAG(VA3,VA2,NAIA3,NAIA2)+NAG(VA3,VA1,NAIA3,NAIA1) INAGB2 = NAG(VA2,VA1,NAIA2,NAIA1) INAGB1 = 0 INAGAP1 = INAGA1+INAGB1 INAGAP2 = INAGA2+INAGB2 INAGAP3 = INAGA3+INAGB3 INAGAP4 = INAGA4+INAGB4 INAGAP5 = INAGA5+INAGB5 INAGAP6 = INAGA6+INAGB6 INAGAP7 = INAGA7+INAGB7 INAGAP8 = INAGA8+INAGB8 INAGAP9 = INAGA9+INAGB9 INAGAP10 = INAGA10+INAGB10 INAGAP11 = INAGA11+INAGB11 INAGAP12 = INAGA12+INAGB12 INAGAP13 = INAGA13+INAGB13 INAGAP14 = INAGA14+INAGB14 INAGAP15 = INAGA15+INAGB15 INAGAP16 = INAGA16+INAGB16 INAGAP17 = INAGA17+INAGB17 INAGAP18 = INAGA18+INAGB18 INAGAP19 = INAGA19+INAGB19 INAGAP20 = INAGA20+INAGB20 INAGAP21 = INAGA21+INAGB21 INAGAP22 = INAGA22+INAGB22 INAGAP23 = INAGA23+INAGB23 INAGAP24 = INAGA24+INAGB24 INAGAP25 = INAGA25+INAGB25 INAGAP26 = INAGA26+INAGB26 INAGAP27 = INAGA27+INAGB27 INAGAP28 = INAGA28+INAGB28 INAGAP29 = INAGA29+INAGB29 INAGAP30 = INAGA30+INAGB30 IGAP1 = IKGAP1+INAGAP1 IGAP2 = IKGAP2+INAGAP2 IGAP3 = IKGAP3+INAGAP3 IGAP4 = IKGAP4+INAGAP4 IGAP5 = IKGAP5+INAGAP5 IGAP6 = IKGAP6+INAGAP6 IGAP7 = IKGAP7+INAGAP7 IGAP8 = IKGAP8+INAGAP8 IGAP9 = IKGAP9+INAGAP9 IGAP10 = IKGAP10+INAGAP10 IGAP11 = IKGAP11+INAGAP11 IGAP12 = IKGAP12+INAGAP12 IGAP13 = IKGAP13+INAGAP13 IGAP14 = IKGAP14+INAGAP14 IGAP15 = IKGAP15+INAGAP15 IGAP16 = IKGAP16+INAGAP16 IGAP17 = IKGAP17+INAGAP17 IGAP18 = IKGAP18+INAGAP18 IGAP19 = IKGAP19+INAGAP19 IGAP20 = IKGAP20+INAGAP20 IGAP21 = IKGAP21+INAGAP21 IGAP22 = IKGAP22+INAGAP22 IGAP23 = IKGAP23+INAGAP23 IGAP24 = IKGAP24+INAGAP24 IGAP25 = IKGAP25+INAGAP25 IGAP26 = IKGAP26+INAGAP26 IGAP27 = IKGAP27+INAGAP27 IGAP28 = IKGAP28+INAGAP28 IGAP29 = IKGAP29+INAGAP29 IGAP30 = IKGAP30+INAGAP30 User-defined functions: MINF(V) = 1./(1.+EXP(-(V-THM)/SIGM)) INA(V,N,VNA) = GNA*MINF(V)^3*(1-N)*(V-VNA) MINFP(V) = 1./(1.+EXP(-(V-THMP)/SIGMP)) HINFP(V) = 1./(1.+EXP(-(V-THHP)/SIGHP)) TAUHP(V) = TAUBAR/COSH((V-VT)/(2*SIG)) INAP(V,HP,VNA) = GNAP*MINFP(V)*HP*(V-VNA) NINF(V) = 1./(1.+EXP(-(V-THN)/SGN)) TAUN(V) = TAUN0+TAUN1/(1+EXP(-(V-THNT)/SN)) IK(V,N,VK) = GK*(N^4)*(V-VK) IL(V) = GL*(V-VL) PHIA(V) = V/FRT PHIGAP(V1,V2) = (V1-V2)/FRT EK(KE,KI) = FRT*LOG(KE/KI) IKA(V,KE,KI) = (1-GAMMA)*PKA*F*PHIA(V)*(KE*EXP(-PHIA(V))-KI)/(EXP(-PHIA(V))-1) INAA(V,NAE,NAI) = PNAA*F*PHIA(V)*(NAE*EXP(-PHIA(V))-NAI)/(EXP(-PHIA(V))-1) IPUMPA(KE,NAI) = RHOA/(((1+2/KE)^2)*(1+10/NAI)^3) GAP(X,Y,A,B) = F*PHIGAP(X,Y)*((B*EXP(-PHIGAP(X,Y))-A)/(EXP(-PHIGAP(X,Y))-1)) KG(V1,V2,KI1,KI2) = PKGAP*GAP(V1,V2,KI1,KI2) NAG(V1,V2,NAI1,NAI2) = PNAGAP*GAP(V1,V2,NAI1,NAI2)