lprint(`Úloha 109:`);lprint(`Funkci f aproximujte jen s pomocí tabulky hodnot a čtyř základních`);
lprint(`aritmetických operací (eventuálně umocňování s celočíselným`);
lprint(`exponentem) tak, abyste dosáhli relativní chybu nejvýše epsilon na`);
lprint(`intervalu nejméně (a,b), kde`);

f := erf; erf('xx') = Int(2/sqrt(Pi) * e^(-'t'^2), 't'=0..'xx');
a := -30; 
b := 30;
epsilon := 10^(-6);