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

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