lprint(`Uloha 604:`);
lprint(`Numericky reste diferencialni rovnici  y'(x) = f(x,y(x))                    `);
lprint(`s pocatecni podminkou  y(x0) = y0  na intervalu  [x0,xk]                    `);
lprint(`s presností alespon epsilon.                                                `);
lprint(`Najdete nejmensi pocet kroku pro dosazeni dane presnosti.                   `);
lprint(`Zduvodnete dosazenou presnost a predvedte pouziti Richardsonovy extrapolace.`);

f := (x,y) -> 8*exp(-x^2) - y;
x0:=0.;
y0:=1.;
xk:=4.;
epsilon:=1e-5;