第4周项目2—太乐了

/*copyright(c)2016.烟台大学计算机学院

* All rights reserved

* 文件名称：my dream，Cpp

* 作者：杨驰

* 完成日期：2016年3月20日

* 问题描述： 用sin泰勒展式，求出sin(π/2)和sin(56°)的值，精确度要达到小数点后六位
           (即当最后一项的绝对值小于0.00001时，累加结束，求绝对值的函数也可以
            自定义函数myabs实现)。（泰勒展式见附录）
*/

#include<iostream>
using namespace std;
const double pi=3.1415926;
double mysin(double x);
int factor(int x);
int main()
{
   cout<<"sin(π/2)的值为"<<mysin(pi/2)<<endl;
   cout<<"sin(56°)的值为"<<mysin(56.0/180.0*pi)<<endl;
   return 0;
}
double mysin(double x)
{
    double r=0.0;
    double e=x;
    double sqr=x*x;
    int n=1;
    while(e/factor(n)>=1e-6)
    {
        double f=e/factor(n);
        r=(n%4==1)?r+f:r-f;
        e=e*sqr;
        n+=2;
    }
    return r;
}
int factor(int x)
{
   if(x==1)
    return 1;
   else
    return x*factor(x-1);
}

附录：

运行结果：