35#include "FitzHughNagumo1961OdeSystem.hpp"
36#include "OdeSystemInformation.hpp"
48 boost::shared_ptr<AbstractIvpOdeSolver> pOdeSolver,
49 boost::shared_ptr<AbstractStimulusFunction> pIntracellularStimulus)
63 double membrane_V = rY[0];
64 double recovery_variable = rY[1];
69 double membrane_V_prime = 0;
73 membrane_V_prime = membrane_V*(membrane_V-
mAlpha)*(1-membrane_V)-recovery_variable+i_stim;
77 double recovery_variable_prime =
mEpsilon*(membrane_V-
mGamma*recovery_variable);
79 rDY[0] = membrane_V_prime;
80 rDY[1] = recovery_variable_prime;
87 double recovery_variable = (*pStateVariables)[1];
88 double fake_ionic_current = membrane_V*(membrane_V-
mAlpha)*(1-membrane_V)-recovery_variable;
89 return fake_ionic_current;
98 this->mVariableNames.push_back(
"V");
99 this->mVariableUnits.push_back(
"mV");
100 this->mInitialConditions.push_back(0.0);
102 this->mVariableNames.push_back(
"w");
103 this->mVariableUnits.push_back(
"");
104 this->mInitialConditions.push_back(0.0);
106 this->mInitialised =
true;
bool mSetVoltageDerivativeToZero
double GetIntracellularAreaStimulus(double time)
std::vector< double > mStateVariables
boost::shared_ptr< AbstractOdeSystemInformation > mpSystemInfo
void EvaluateYDerivatives(double time, const std::vector< double > &rY, std::vector< double > &rDY)
double GetIIonic(const std::vector< double > *pStateVariables=NULL)
~FitzHughNagumo1961OdeSystem()
FitzHughNagumo1961OdeSystem(boost::shared_ptr< AbstractIvpOdeSolver > pOdeSolver, boost::shared_ptr< AbstractStimulusFunction > pIntracellularStimulus)
static const double mGamma
static const double mEpsilon
static const double mAlpha