神经元 ML模型

/**************************************************/
/*********************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<=