#include "modelQuadraticNeuron.h"
#include "math.h"
#include "iostream.h"
void
modelQuadraticNeuron (double t, double *x, double *dx, double *parameter,
		 double *extra)
{
  double v,  b, vreset,vpeak,k,c,vrest;
  double type;
  b = parameter[0];
  vreset = parameter[1];
  vpeak = parameter[2];
 k=parameter[3];
 vrest=parameter[4];
 c=parameter[5];
x[1]=0; 
dx[0] =k*pow((x[0]-parameter[4]),2.0)+parameter[5]*b+extra[1];

  if (x[0]>vpeak)
  {
  x[0]=vreset;
  x[1]=1;
  }

}