@Axoren An ODE system is written as follows: Let $m \in \mathbb{N}, f:[a,b] \times \mathbb{R}^m \to \mathbb{R}^m$ and $y_0 \in \mathbb{R}^m$. We are looking for a function $y:[a,b] \to \mathbb{R}^m$ such that $\left\{\begin{matrix}
y'(t)=f(t,y(t)) &, t \in [a,b] \\
y(a)=y_0 &
\end{matrix}\right.$. For example in this case Euler's method is $y_k^{n+1}=y_k^n+hf(t^n,y_k^n)$