/**************************************************/
/*********************H-H?model********************/
/**************************************************/
/****************equation?form???f=dx/dt***********/
/**************************************************/

#include?
#include?
#include?

double?xx[3];
double?f[3];
int?neq=2;
long?int?nstep?=?5000;
double?dt=0.01;

double?eps=20.0;?
???
double??V1=-1.2;
double??V2=18;
double??V3=2;
double??V4=30;

double?Vk=-84;??????
double?Vca=120;
double?Vl=-60;????

double??gl=2.0;
double??gk=8.0;????????
double??gca=4.4;???????

double?I;??????//parameter?bifurcation?（left?subcritical?Hoft?:stable?and?unstable?cycle-------[88.293193.85761]）
double??fai=0.04;

double?minf（double?x）;
double?winf（double?x）;
double?tao（double?x）;



void??ftn（）;

main（）
{
?FILE?*fi;
?FILE?*ff;
?FILE?*fp;
?long??int?j;
?double?f1[3]f2[3]f3[3]f4[3]x0[3];
?int?iisioldisinewkk;
?double?n;
?double?x1x2x3;
//?int?T;





??ff=fopen（“ML.dat““w“）;
??fp=fopen（“ML_ISI.dat““w“）;
??fi=fopen（“ML_initial.dat““w“）;


??for（n=1;n<=15000;n=n+1）
?{
	??I=n;
?printf（“%lf\n“n）;
??

??xx[1]?=?-30;
??xx[2]?=?0.11;


??x1=xx[1];
??x2=xx[1];
??x3=xx[1];

??isiold=0;
??isinew=0;
??kk=0;

???????				
????

??xx[0]?=?0.0;
??x0[0]?=?0.0;


???printf（“please?wait?%lf\n“n）;
?

?for（j=1;j<=500*nstep;j++）
?{
	/*??
????if（j>=20*nstep&&j<=40*nstep）
		I=50;
	else
		I=0;
????*/

??for（i=1;i<=neq;i++）
?????x0[i]?=?xx[i];

?????ftn（）;

??for（i=1;i<=neq;i++）
????{
?????f1[i]=f[i]*dt;
?????xx[i]=x0[i]+f1[i]/2.0;
????}

?????ftn（）;

???for（i=1;i<=