#include "modelIh.h" #include "math.h" void modelIh(double t,double *x,double *dx,double *parameter,double *extra) { double v,hinf,tinf; v=extra[1]; hinf=1./(1+exp(-(v+35.)/10.0)); tinf=3900.0/(exp(-7.68-.086*v)+exp(5.04+.0701*v)); dx[0]=(hinf-x[0])/tinf; } /*m1=1./(1+exp(-(v+60)/8.5)); m2=1./(1+exp(-(v+36)/20)); tma=(1./(exp((v+35.82)/19.69)+exp(-(v+79.69)/12.7))+.37); h1=1./(1+exp((v+78)/6)); tadef=1./(exp((v+46.05)/5.)+exp(-(v+238.4)/37.45)); if(v<-63) ta1=tadef; else ta1=19; if(v<-73) ta2=tadef; else ta2=60; dx[0]=(m1-x[0])/tma; dx[1]=(m2-x[1])/tma; dx[2]=(h1-x[2])/ta1; dx[3]=(h1-x[3])/ta2;*/