lprint(`Uloha 603:`);
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) -> x^2 + y;
x0:=0.;
y0:=-1.;
xk:=3.;
epsilon:=1e-6;