function ov07_5c(h)
%
% Exercise 7, Task 5 c) 
%
y0 = [5210 / 249401; 20259 / 249401];
tend = 10;

[t, y] = impEuler(@ov07, 0, tend, y0, h);
plot(t, y);
xaxis([0, tend]);

function yp = ov07(t, y)
yp = [-10 * y(1) - 10 * y(2) + sin(t);
       40 * y(1) - 10 * y(2) + cos(t)];

