第12周项目2--太乐了

/*
 * Copyright (c) 2013, 烟台大学计算机学院
 * All rights reserved.
 * 作    者：沈远宏
 * 完成日期：2013 年 11 月 13 日
 * 版 本 号：v1.0
 * 问题描述：用in和cos的泰勒展开式编程序，求出下面四个式子的值，当最后一项的绝对值小于0.00001时，累加结束。
 * 问题分析：要求自定义四个函数，分别是mysin,mycos,myabs,fac求sin的泰勒展开式，cos的泰勒展开式，最后一项的绝对值和分母的绝对值。
 */#include <iostream>
using namespace std;
const double PI=3.1415926;
double mysin(double);
double mycos(double);
double myabs(double);
int fac(int);
int main()
{
    cout<<"sin(π/2)的值为"<<mysin(PI/2)<<endl;
    cout<<"sin(56°)的值为"<<mysin((56.0/180)*PI)<<endl;
    cout<<"cos(87°)的值为"<<mycos((87.0/180)*PI)<<endl;
    cout<<"cos(π/3)的值为"<<mycos(PI/3)<<endl;
    return 0;
}
double mysin(double x)
{
    int i=2;
    double sinx=x,xs=x,absxxs=1,a,xxs;//sinx为每一个单项式的加和，xs 为每一项的分子，a 为每一项的分母也就是阶乘，xxs为分子分母的比值。absxxs为xxs的绝对值
    while(absxxs>=0.00001)
    {

        a =fac(2*i-1);
        ++i;
        xs=(xs*x*x*(-1));
        xxs=xs/a;
        sinx=sinx+xxs;
        absxxs=myabs(xxs);
    }
    return sinx;
}
double mycos(double x)
{
    int i=1;
    double cosx=1,xs=1,absxxs=1,a,xxs;
    while(absxxs>=0.00001)
    {
        a=fac(2*i);
        ++i;
        xs=xs*x*x*(-1);
        xxs=xs/a;
        cosx=cosx+xxs;
        absxxs=myabs(xxs);
    }
    return cosx;
}
double myabs(double a)
{
    double n;
    if(a<=0)
        n=(-a);
    else
        n=a;
    return n;
}
int fac(int a)
{
    int fac=1, i=1;
    for(i=1; i<=a; ++i)
    {
        fac=fac*i;
    }
    return fac;
}

运行结果：

心得体会：

首先absxxs是>=0.0001,我写成了<0.0001,导致纠结了很久，其次，这里面的几个循环没看清楚有点晕，所以耽误了很多时间。这说明要保持清醒啊！！！