function RungehoPriklad()

N = 10;

syms x
f = 1/(1+x^2);

xx = -5:0.1:5;
% omezenost derivaci
fakt = 1;
ff = eval(subs(f,x,xx));
plot(xx,ff/fakt); pause
df = f;
for n=1:10
    n
    df = diff(df,x);
    fakt = fakt*n;
    dff = eval(subs(df,x,xx));
    plot(xx,dff/fakt); pause
end
    