4Misc_Start<4Platform<¸@ þÊROGIXX¨ý¨ýp”þþ4<èèèè    4!ø*TÐ<winspoolHP LaserJet Professional P1102w (Copy 1)USB001HP LaserJet Professional P1102œ4Ûê oXLetterSDDMHP LaserJet Professional P1102w-€€€ÿÿÿÿÿÿÿÿÿ(dÞÞ€R²¼4„ þÊROGIXX¨ý¨ýp”þþ4<èèèè    4!ø*TÐ<winspoolHP LaserJet Professional P1102w (Copy 1)USB001HP LaserJet Professional P1102œ4Ûê oXLetterSDDMHP LaserJet Professional P1102w-€€€ÿÿÿÿÿÿÿÿÿ(dÞÞ€R²¼4„ þÊROGIXX¨ý¨ýp”þþ4<èèèè    4!ø*TÐ<winspoolHP LaserJet Professional P1102w (Copy 1)USB001HP LaserJet Professional P1102œ4Ûê oXLetterSDDMHP LaserJet Professional P1102w-€€€ÿÿÿÿÿÿÿÿÿ(dÞÞ€R²¼4„ þÊ ROGIXX¨ý¨ýp”þþ4<èèèè    4!ø*TÐ<winspoolHP LaserJet Professional P1102w (Copy 1)USB001HP LaserJet Professional P1102œ4Ûê oXLetterSDDMHP LaserJet Professional P1102w-€€€ÿÿÿÿÿÿÿÿÿ(dÞÞ€R²¼4„^Graph*@@ð?ð? WDashSettings#úÿ  !¼ˆºÀüNormalÿÿÿÿ@ ÿArialÿÿÿÿÿÿ<ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿHHHH$$¼ˆºÀüNormalÿÿÿÿ@ ÿArialÿÿÿÿÿÿ<ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿHHHH$$444444¼ˆ +QùüNormalÿÿÿÿ@ ÿArialÿÿÿÿÿÿ<ÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿÿHHHH$$4 4 4 4 4 4 homeMEdÿÿÿÿû¨™·´C:Users:Brian:Desktop:ionchannelmodel:C:Users:Brian:Desktop:ionchannelmodelÜç 2ñ—þP]× p]× €çcÞwaøwuIußé°ëÐí$ÿÿÿÿ„çˆç€ü°^uâÄœþÿÿÿIu„uÎë-= ctŠHéä1ûwuÎë;C{²ãéÙ‰ionchannelmodelMEdÿÿÿÿû¨™·´C:Users:Brian:Desktop:ionchannelmodel:C:Users:Brian:Desktop:ionchannelmodelÜç 3ñ—þP]× p]× €çcÞwaøwuIußé°ëÐí$ÿÿÿÿ„çˆç€ü°^uâÄœþÿÿÿIu„uÎë-= ctŠHéä1ûwuÎë;C{²ãéÙ‰?RecentWindows HH (4).ipf 4Misc_End<”XOPState_Start`<Data Browser°aä HH (4).ipfGizmondowPeakFunctions28;4XOPState_End<\Ú•RunSimulation() Vm= mV -65 t= ms 0 n 0.317677 m 8.67996e-12 h 0.792612 Vm= mV -65.1191 t= ms 0.101 n 0.31766 m 8.65832e-12 h 0.792625 Vm= mV -65.2293 t= ms 0.201 n 0.317611 m 8.60627e-12 h 0.792663 Vm= mV -65.331 t= ms 0.3 n 0.317535 m 8.53736e-12 h 0.792722 Vm= mV -65.4267 t= ms 0.4 n 0.317433 m 8.45866e-12 h 0.792803 Vm= mV -65.5156 t= ms 0.5 n 0.317307 m 8.37637e-12 h 0.792902 Vm= mV 30.8514 t= ms 0.6 n 0.34297 m -1.39339e-08 h 0.0312071 Vm= mV 21.5173 t= ms 0.7 n 0.392933 m -1.03521e-07 h 7.86194e-06 Vm= mV 9.64498 t= ms 0.8 n 0.432592 m -1.29707e-07 h 3.59339e-05 Vm= mV -3.6904 t= ms 0.9 n 0.462321 m -1.344e-07 h 0.000151829 Vm= mV -17.0276 t= ms 1 n 0.483081 m -1.32829e-07 h 0.000499094 Vm= mV -29.1378 t= ms 1.1 n 0.496349 m -1.27768e-07 h 0.00128441 Vm= mV -39.4494 t= ms 1.201 n 0.503921 m -1.18515e-07 h 0.00276386 Vm= mV -47.628 t= ms 1.301 n 0.507326 m -1.04698e-07 h 0.00513108 Vm= mV -53.9787 t= ms 1.401 n 0.507989 m -8.7041e-08 h 0.00854397 Vm= mV -58.8176 t= ms 1.501 n 0.506845 m -6.77159e-08 h 0.0130446 Vm= mV -62.4668 t= ms 1.601 n 0.504511 m -4.93224e-08 h 0.0185783 Vm= mV -65.2056 t= ms 1.701 n 0.501389 m -3.38046e-08 h 0.0250254 Vm= mV -67.2575 t= ms 1.801 n 0.497743 m -2.19677e-08 h 0.032236 Vm= mV -68.7944 t= ms 1.901 n 0.493749 m -1.36499e-08 h 0.0400563 Vm= mV -69.946 t= ms 2.001 n 0.489528 m -8.17614e-09 h 0.0483442 Vm= mV -70.8086 t= ms 2.101 n 0.485164 m -4.75545e-09 h 0.0569771 Vm= mV -71.4539 t= ms 2.201 n 0.480718 m -2.70226e-09 h 0.0658526 Vm= mV -71.9346 t= ms 2.301 n 0.476232 m -1.50771e-09 h 0.0748882 Vm= mV -72.2902 t= ms 2.401 n 0.471738 m -8.29137e-10 h 0.0840182 Vm= mV -72.5477 t= ms 2.5 n 0.467304 m -4.53405e-10 h 0.093099 Vm= mV -72.734 t= ms 2.6 n 0.462858 m -2.43887e-10 h 0.102274 Vm= mV -72.8627 t= ms 2.7 n 0.458456 m -1.29771e-10 h 0.11142 Vm= mV -72.9463 t= ms 2.8 n 0.454109 m -6.8081e-11 h 0.120514 Vm= mV -72.9942 t= ms 2.9 n 0.449824 m -3.49074e-11 h 0.129536 Vm= mV -73.0139 t= ms 3 n 0.445606 m -1.71311e-11 h 0.138473 Vm= mV -73.0108 t= ms 3.1 n 0.441459 m -7.62406e-12 h 0.147312 Vm= mV -72.9896 t= ms 3.2 n 0.437387 m -2.54186e-12 h 0.156048 Vm= mV -72.9534 t= ms 3.3 n 0.43339 m 1.78685e-13 h 0.164672 Vm= mV -72.9051 t= ms 3.4 n 0.429471 m 1.64095e-12 h 0.17318 Vm= mV -72.8469 t= ms 3.5 n 0.42563 m 2.4337e-12 h 0.181571 Vm= mV -72.7803 t= ms 3.6 n 0.421868 m 2.8707e-12 h 0.18984 Vm= mV -72.7068 t= ms 3.7 n 0.418186 m 3.11904e-12 h 0.197987 Vm= mV -72.6276 t= ms 3.8 n 0.414582 m 3.26772e-12 h 0.206012 Vm= mV -72.5435 t= ms 3.9 n 0.411058 m 3.3641e-12 h 0.213913 Vm= mV -72.4554 t= ms 4 n 0.407613 m 3.43339e-12 h 0.221691 Vm= mV -72.3637 t= ms 4.1 n 0.404245 m 3.48902e-12 h 0.229347 Vm= mV -72.2691 t= ms 4.2 n 0.400955 m 3.53808e-12 h 0.236881 Vm= mV -72.1721 t= ms 4.3 n 0.397742 m 3.58434e-12 h 0.244294 Vm= mV -72.0729 t= ms 4.4 n 0.394606 m 3.62975e-12 h 0.251588 Vm= mV -71.972 t= ms 4.5 n 0.391544 m 3.67533e-12 h 0.258764 Vm= mV -71.8696 t= ms 4.6 n 0.388557 m 3.72157e-12 h 0.265823 Vm= mV -71.7659 t= ms 4.7 n 0.385644 m 3.7687e-12 h 0.272766 Vm= mV -71.6612 t= ms 4.8 n 0.382803 m 3.81682e-12 h 0.279595 Vm= mV -71.5557 t= ms 4.9 n 0.380035 m 3.86593e-12 h 0.286313 !À·í*þd TXET????dHH (4).ipfHH (4).ipf|========================================================= | | Hodgkin-Huxley action potential simulation | | | by Jiahua Li 12.2.2010 |========================================================= Menu "Macros" "Run Simulation" End Function RunSimulation() variable Vrest=-65, I_inj=1 Prompt Vrest, "Resting potential (mV)?" Prompt I_inj, "the stimulus curreut is ? (mA)" // all values are SI units // wave and variable declaration variable Cmemb=1e-6, Ttot=5e-3, vNa=0.05, vK=-0.077, vL=0.0106, Vm variable gbarNa=0.12, gbarK=0.036, gbarL=0.0003, I_inject variable/D INa, IK, IL, Imem variable/D dT= 1e-6,t=0,tprint=0,tpstep=1e-4 variable/D dndt, dmdt, dhdt, dvdt variable numticks=floor(Ttot/dT) variable Ttotms variable n,m,h, an, am, ah, bn, bm, bh Vrest*=0.001 I_inj*=0.001 // initial a, b constant Vm=Vrest an= 10*(-1000*(Vm-Vrest)+10)/(exp((-1000*(Vm-Vrest)+10)/10)-1) // Vm and Vrest are in mV, an in s-1 bn=125*exp(-1000*(Vm-Vrest)/80) am=100*(25-1000*(Vm-Vrest))/(exp(25-1000*(Vm-Vrest)/10)-1) bm=4000*exp(-1000*(Vm-Vrest)/18) ah=70*exp(-1000*(Vm-Vrest)/20) bh=1000/(exp((30-1000*(Vm-Vrest))/10+1)) n=an/(an+bn) m=am/(am+bm) h=ah/(ah+bh) do if ((t>=5e-4)&(t<6e-4)) // inject current happen at 0.5-0.6ms I_inject=I_inj else I_inject=0 endif if (t>=tprint) //// print the result every 0.1ms print "Vm= mV ", Vm*1000 ,"t= ms ", t*1000, "n ", n, "m ",m,"h ",h tprint+=tpstep Endif INa=gbarNa*m*m*m*h*(Vm-vNa) IK=gbarK*n*n*n*n*(Vm-vK) IL=gbarL*(Vm-Vrest-vL) Imem=INa+IK+IL dvdt=(-Imem+I_inject)/Cmemb // use Imem to calculate dvdt Vm+=dvdt*dT // calculate new Vm an= 10*(-1000*(Vm-Vrest)+10)/(exp((-1000*(Vm-Vrest)+10)/10)-1) // Vm and Vrest are in mV, an in s-1 bn=125*exp(-1000*(Vm-Vrest)/80) am=100*(25-1000*(Vm-Vrest))/(exp(25-1000*(Vm-Vrest)/10)-1) bm=4000*exp(-1000*(Vm-Vrest)/18) ah=70*exp(-1000*(Vm-Vrest)/20) bh=1000/(exp((30-1000*(Vm-Vrest))/10+1)) dndt=an*(1-n)-bn*n dmdt=am*(1-m)-bm*m dhdt=ah*(1-h)-bh*h n+=dndt*dT // calculate new n, m, h m+=dmdt*dT h+=dhdt*dT t+=dT while (t<=Ttot) EndJ// Platform=WindowsNT, IGORVersion=6.204, architecture=Intel Silent 101 // use | as bitwise or -- not comment. NewPath/Z ionchannelmodel "C:Users:Brian:Desktop:ionchannelmodel:" DefaultFont "Arial" MoveWindow/P 5.25,42.5,504.75,337.25 OpenProc/W=(198,42.5,697.5,515.75)/J=15936 "HH (4).ipf" MoveWindow/C 8.25,442.25,960,536.75 :#pragma rtGlobals=1 // Use modern global access method.