function f=faktorial_r(n)

if n<0 || fix(n)~=n
    error ('Vstupem musi byt neyaporne cele cislo!')
end

if n==0
    f=1;
    return
end

f=n;
if n==1
    return
else
f=f*faktorial_r(n-1);
end