Take the 2-minute tour ×
Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. It's 100% free, no registration required.

Consider a system of the form:

(1) $x' = Ax + g$

For appropriate matrices $x'$, $A$, $x$, and $g$.

If we let $y_n$ be the approximation to the solution of (1) at time step $t_n$, what matrix equation will give $y_{n+1}$ for Euler's method? How about the improved Euler method (I understand this case may involve two equations)?

share|improve this question

1 Answer 1

up vote 1 down vote accepted

Euler's Method works as follows: Given an ODE $y'(t) = f(t,y(t)$ we can compute $y_{n+1} = y_n +h f(t_n,y(t_n)) $ with some given stepsize $h$. In your case ( even though you are speaking of $y_{n+1}$, but your ODE reads $x'=Ax+g$, I assume you then have the ODE $y'=Ay+g$ ) $y_{n+1} = y_n +h (A y_n +g) = (A+h)y_n+ hg $.

Improved Euler method can be split into two equation, but it should rather be called two steps. You basically perform one Euler step, but this time you call the solution $y^*_{n+1}$. You then use the average of the slope at $y_n$ and $y^*_{n+1}$ to approximate the slope for your step $y_{n+1} = y_{n} +h\cdot s$ with $s = (f(y_n)+ f(y_{n+1}^*)/2$. Applied to your problem this would result in the following term. \begin{align}y_{n+1} &=y_n + h\cdot(f(y_n)+ f(y_{n+1}^*)/2 \\ &= y_n+ \frac{h}{2}( A y_n+g+(A+h)y_n+hg) \text{ by Euler method for }y^*_{n+1} \\ &= (1+{Ah}+ \frac{ h^2}{2})y_n +(\frac{h}{2}+\frac{h^2}{2})g\end{align}

See also Wiki for improved Euler Method aka Heun's Method.

Thomas

share|improve this answer

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.