泰勒公式

/*
*Copyright (c) 2016,烟台大学计算机学院
*All rights reserved.
*文件名称：text.cpp
*作者：汤善晔
*完成日期：2016年8月23日
*版本号：v1.0
 */

#include <iostream>     

using namespace std;    
const double pi=3.1415926;    
double mysin(double x);    
double myabs(double x);    
int main( )    
{    
    cout<<"sin(π/2)的值为"<<mysin(pi/2)<<endl;    
    cout<<"sin(56°)的值为"<<mysin((56.0/180)*pi)<<endl;    
   return 0;    
}    
double mysin(double x)    
{    
    double sum=x,x_pow=x,item;    
    int n=1,fact=1,sign=1;      
    do    
    {    
        fact=fact*(n+1)*(n+2);      
        x_pow*=x*x;                 
        sign=-sign;                
        item =x_pow/fact*sign;    
        sum+=item;                   
        n+=2;    
    }while(myabs(item)>1e-5);    
    return sum;    
}        
double myabs(double x)    
{    
    return ((x>=0)?x:-x);    

}