--description--

Euler's method numerically approximates solutions of first-order ordinary differential equations (ODEs) with a given initial value. It is an explicit method for solving initial value problems (IVPs), as described in <a href="https://www.freecodecamp.org/news/eulers-method-explained-with-examples/" title="Euler's Method Explained with Examples" target="_blank" rel="noopener noreferrer nofollow">this article</a>.

The ODE has to be provided in the following form:

with an initial value

To get a numeric solution, we replace the derivative on the LHS with a finite difference approximation:

<ul style='list-style: none;'> <li><big>$\frac{dy(t)}{dt} \approx \frac{y(t+h)-y(t)}{h}$</big></li> </ul>

then solve for $y(t+h)$:

<ul style='list-style: none;'> <li><big>$y(t+h) \approx y(t) + h \, \frac{dy(t)}{dt}$</big></li> </ul>

which is the same as

<ul style='list-style: none;'> <li><big>$y(t+h) \approx y(t) + h \, f(t,y(t))$</big></li> </ul>

The iterative solution rule is then:

where $h$ is the step size, the most relevant parameter for accuracy of the solution. A smaller step size increases accuracy but also the computation cost, so it has always has to be hand-picked according to the problem at hand.

Example: Newton's Cooling Law

Newton's cooling law describes how an object of initial temperature $T(t_0) = T_0$ cools down in an environment of temperature $T_R$:

<ul style='list-style: none;'> <li><big>$\frac{dT(t)}{dt} = -k \, \Delta T$</big></li> </ul>

It says that the cooling rate $\frac{dT(t)}{dt}$ of the object is proportional to the current temperature difference $\Delta T = (T(t) - T_R)$ to the surrounding environment.

The analytical solution, which we will compare to the numerical approximation, is

--instructions--

Implement a routine of Euler's method and then use it to solve the given example of Newton's cooling law for three different step sizes of:

and compare with the analytical solution.

Initial values:

<ul> <li>initial temperature <big>$T_0$</big> shall be <code>100 °C</code></li> <li>room temperature <big>$T_R$</big> shall be <code>20 °C</code></li> <li>cooling constant <big>$k$</big> shall be <code>0.07</code></li> <li>time interval to calculate shall be from <code>0 s</code> to <code>100 s</code></li> </ul>

First parameter to the function is initial time, second parameter is initial temperature, third parameter is elapsed time and fourth parameter is step size.

--hints--

eulersMethod should be a function.

assert(typeof eulersMethod === 'function');

eulersMethod(0, 100, 100, 2) should return a number.

assert(typeof eulersMethod(0, 100, 100, 2) === 'number');

eulersMethod(0, 100, 100, 2) should return 20.0424631833732.

assert.equal(eulersMethod(0, 100, 100, 2), 20.0424631833732);

eulersMethod(0, 100, 100, 5) should return 20.01449963666907.

assert.equal(eulersMethod(0, 100, 100, 5), 20.01449963666907);

eulersMethod(0, 100, 100, 10) should return 20.000472392.

assert.equal(eulersMethod(0, 100, 100, 10), 20.000472392);

--seed--

--seed-contents--

function eulersMethod(x1, y1, x2, h) {

}

--solutions--

function eulersMethod(x1, y1, x2, h) {
  let x = x1;
  let y = y1;

  while ((x < x2 && x1 < x2) || (x > x2 && x1 > x2)) {
    y += h * (-0.07 * (y - 20));
    x += h;
  }

  return y;
}