Синус и косинус, ряд Тейлора

#include <iostream>
#define _USE_MATH_DEFINES
#include <math.h>
using namespace std;
int main()
{
    int i,s;
    double z,n,x0,x,n2,y1,y2;
    cout.precision(18);
    cout<<"vvedite 4islo  ";
    cin>>x;
    x0=x;y2=0;y1=0;n=x;n2=1;s=0;z=0;
    if ((x<M_PI_2)&&(x>M_PI_4))
      {
          x=fmod(x,M_PI_4);
           n=1;n2=x;
    }
    else 
        {
            if ((x>M_PI_2)&&(x<M_PI))
            {
                x=fmod(x,M_PI_4);
                n=-x;n2=1;
            }
            else if((x>M_PI)&&(x<3*M_PI_2))
                {
                    x=fmod(x,M_PI_4);
                n=1;n2=x; 
            }
            else x=fmod(x,M_PI);
    }
    cout<<"\nzna4enie M_PI="<<M_PI;
    cout<<"\nzna4enie x="<<x;
    for (i=1;;i++)
    {
        y1+=n;
        y2+=n2;
        n*=-x*x/((2*i)*(2*i+1));
        n2*=-x*x/((2*i)*(2*i-1));
        s++;
        if ((y1==y1+n)&&(y2==y2+n2))
        break;
    }
    if ((x0>M_PI_2)&&(x0<M_PI))
        {
        z=y1;y1=y2;y2=z;
    }
       cout<<"     summ bilo  "<<s;
       cout<<"\nsin 4erez funkjiu="<<sin(x0);
       cout<<"\nsin 4erez pi formulu="<<y1;
       cout<<"\ncos 4erez funkjiu="<<cos(x0);
       cout<<"\ncos 4erez pi formulu="<<y2;
    cin.get();
    cin.get();
    return 0;
}

while(x>2*M_PI) x-=2*M_PI;