function s_p = Hausaufgabe3( xe )
%HAUSAUFGABE3 Berechnung der Sinus-Funktion durch
%   schrittweise Reduktion auf das Intervall [0, pi/6]
x=xe;
s_v=1;

% Reduktion auf [0, 2*pi]
if (x>=2*pi)
    x=x-2*pi*floor(x/(2*pi));
end

% Reduktion auf [0, pi]
if (pi<=x)
    x=x-pi;
    s_v=-1*s_v;
end

% Reduktion auf [0, pi/2]
if (pi/2<=x)
    x=pi-x;    
end

% Reduktion auf [0, pi/6]
x=0.333333333333333*x;
sindrittel=x-0.166666666666667*x.^3+0.008333333333333*x.^5;

s_p=3*sindrittel-4*sindrittel.^3;
s_p=s_v*s_p;
    
end
