#include "Euler.h"
Euler::Euler()
{
}
vec Euler::solve(std::vector< functPtr > functions)
{
vec solution;
functPtr temp_ptr;
for (int i = 0; i < functions.size(); i++)
{
temp_ptr = functions[i];
solution.push_back(this->initialConditions[i] + this->timeStep*temp_ptr(this->timeStep, this->initialConditions));
}
return solution;
}
void Euler::setTimeStep(double time)
{
this->timeStep = time;
}
void Euler::setInitial(std::vector< double > initial)
{
this->initialConditions.clear();
for(int i = 0; i < initial.size(); i++)
this->initialConditions.push_back(initial[i]);
}