SIN_function_taylor_series(15)

function[y]=SIN(x)
%%First step is to bring x between -2pi and 2pi

x=x-2*pi*fix((x/(2*pi)));

%%Second step is to bring x between -pi and ppi
if(x<-pi)
    x=x+2*pi;
elseif(x>pi)
    x=x-2*pi;
end

    factorijel=1.;
    predznak=1;
    accX=1.;
    y=0.;
    N=50;
    for i=1:N
        factorijel=factorijel*i;
        accX=accX*x;
        if(mod(i,2)~=0)
            y=y+accX*predznak/double(factorijel);
            predznak=predznak*(-1);
        end
    end
   y=round(y*10000)/10000.;
end