#include "modelHenrysynapse.h"
#include "math.h"


void modelHenrysynapse(double t,double *x,double *dx,double *parameter,double
*extra)
{
double So,tau,s1;
 tau =parameter[3]-parameter[2]; 
  s1 = parameter[3]/tau;
 So = 1. / 2 * (1 + tanh (120 * (extra[1] - .1)));
//So = 1. / 2 * (1 + tanh (120 * (extra[1] +25)));

dx[0]=(1./tau)*(So-x[0])/(s1-So);
}


//the 'extra' parameter signifies the the ode of the object needs input from other
//object at everytime tims step.