0
$\begingroup$

Consider the IVP $$y' = \frac{1}{t^2} - \frac{y}{t} = y^2$$ with $t$ an element of $[1,2]$ and $y(1) = -1$. The exact solution is $y(t) = -1/t$

(a) How can we use Euler's method with $h = 0.05$ to approximate the solution & and compare it to actual values of $y$?

(b) Use answer in part (a), and linear interpolation to approx the following values of $y$ and compare also with actual values:

$y(1.052)$

$y(1.555)$

$y(1.978)$

I'm stuck at the euler setup in part (a), and also, with the portion of comparing errors. Can you work out for one of the values in part (b), so I can get it?

1 Answers 1

1

Here is some pseudo code for (a):

h := 0.05;
// The number of xs minus 1 is the number of steps you'll do:
number_of_xs := length([1: h : 2]) - 1;

y_Euler(0) := -1; // Initialise your array of ys.
for i = 0 to i = (number_of_xs-1) do
  y_Euler(i+1) = y_Euler(i) + h*y_Euler(i)^2;

// Then you want to compute the exact values:
xs := [1: h : 2];
y_exact := -1/xs; // Alternatively, write a for loop.

// Then you want to compute the (local) error:
err := abs(y_Euler - y_exact);

You might want to read about local and global error in the second paragraph here. Hope this helps.

  • 0
    Any ideas on how I can do this on Maple?2013-10-31