第十二周上机实践 泰勒公式

#include <iostream>
#include<Cmath>  //为便于对比结果，main函数中调用了Cmath中的库函数sin和cos
using namespace std;
const double pi=3.1415926;
double mysin(double);
double mycos(double);
double myabs(double);   //程序中需要求精度的绝对值，也用自定义函数完成吧
int main( )
{
    cout<<"sin(π/2)的值为"<<mysin(pi/2);


    cout<<"sin(56°)的值为"<<mysin((56.0/180)*pi);

    cout<<"cos(87°)的值为"<<mycos((87.0/180)*pi);

    cout<<"cos(π/3)的值为"<<mycos(pi/3);


    return 0;
}
//下面定义mysin函数
double mysin(double x)
{
    double sum=x,x_pow=x,item;
    int n=1,fact=1,sign=1;     //定义变量时赋初值，已经将第一项考虑到累加和sum中
    do
    {
        fact=fact*(n+1)*(n+2);  //fact用于表示阶乘，在公式中作分母
        x_pow*=x*x;             //x_pow是分子中用于表示阶乘，在公式中作分母
        sign=-sign;             //确定即将要累加的这一项的符号
        item =x_pow/fact*sign; //计算出要累加的项
        sum+=item;              //将该项累加上去
        n+=2;
    }while(myabs(item)>1e-5);
    return sum;
}

//下面定义mycos函数
double mycos(double x)
{
    double sum=1,x_pow=1,item;
    int n=0,fact=1,sign=-1;
    do
    {
        fact=fact*(n+1)*(n+2);
        x_pow*=x*x;
        item =x_pow/fact*sign;
        sum+=item;
        sign=-sign;
        n+=2;
    }while(myabs(item)>0.00001);
    return sum;
}

//下面定义myabs函数
double myabs(double x)
{
    return ((x>=0)?x:-x);
}