diff --git a/ODElibrary/ForwardEuler.cpp b/ODElibrary/ForwardEuler.cpp index 14bae88..bda6609 100644 --- a/ODElibrary/ForwardEuler.cpp +++ b/ODElibrary/ForwardEuler.cpp @@ -1,3 +1,14 @@ #include +#include "ForwardEuler.h" + +void ForwardEuler::solve(){ + auto y = initialValue; + double t = initialTime; + while(t< finalTime){ + y = y + stepSize * (A*y + g(t)); // the simple Euler update + t += stepSize; + std::cout << t << " " << y << std::endl; + } +}; diff --git a/ODElibrary/ForwardEuler.h b/ODElibrary/ForwardEuler.h index 6127454..29b771c 100644 --- a/ODElibrary/ForwardEuler.h +++ b/ODElibrary/ForwardEuler.h @@ -1,4 +1,10 @@ -#include +#include "AbstractSolver.h" + + +class ForwardEuler : public AbstractSolver{ +public: + void solve() override; +};