第4週目プロジェクト2 sinx値を求める

/*copyright(c)2016.         
* All rights reserved,
*     ：my dream，Cpp
*   ：   
*     ：2016 3 19 
*    ：vc++6.0
*     ：
            sin    ，  sin(π/2) sin(56°)  ，            
           (            0.00001 ，    ，          
                 myabs  )。（        ）
*/
#include<iostream>
using namespace std;
const double pi=3.1415926;
int main()
{
    double mysin(double m);
    cout<<"sin(π/2)   "<<mysin(pi/2)<<endl;
    cout<<"sin(56°)   "<<mysin((pi/180.0)*56.0)<<endl;
    return 0;
}
double mysin(double m)
{
    int di(int x);
    double sq=m*m,s=0.0,j,d;
    int i=1,t=1,R;
    R=di(i);
    d=m;
    j=d/(R*1.0);
    s=j;
    while(j>1e-6)
    {
        i+=2;
        R=di(i);
        d=d*sq;
        j=d/(R*1.0);
        t=-t;
        s=s+t*j;
    }
    return s;
}
int di(int x)
{
    if(x==1)
        return 1;
    else
        return x*di(x-1);
}

以下に図面を示す.