function y = sinus(x)
% ukazka vypoctu sumace a 
% implicitniho vypoctu mocnin a faktorialu
diff = 1e-16;
y = x;
t = x;
n = 2;
x2 = -x*x;   
while abs(t) > diff
    t = t*x2/n/(n+1);
    y = y + t;
    n = n + 2;
end
