(PTA)习题5-7 使用函数求余弦函数的近似值【C语言实现】

本题要求实现一个函数，用下列公式求cos(x)的近似值，精确到最后一项的绝对值小于e：

cos(x)=x0/0!−x2/2!+x4/4!−x6/6!+⋯

函数接口定义：

double funcos( double e, double x );

其中用户传入的参数为误差上限e和自变量x；函数funcos应返回用给定公式计算出来、并且满足误差要求的cos(x)的近似值。输入输出均在双精度范围内。

裁判测试程序样例：

#include <stdio.h>
#include <math.h>

double funcos( double e, double x );

int main()
{    
    double e, x;

    scanf("%lf %lf", &e, &x);
    printf("cos(%.2f) = %.6f\n", x, funcos(e, x));
    
    return 0;
}

/* 你的代码将被嵌在这里 */

输入样例：

0.01 -3.14

输出样例：

cos(-3.14) = -0.999899

解题思路：它必须先进行一次计算，然后依次判断最后一项的绝对值是否小于e。如果不先计算一次的话，就没有项与e比较。所以这里使用的是do...while()循环实现。期间需要注意符号的变换和幂方是每项+2的递增，而不是+1，需要注意；然后我分母阶乘实现的方法是重写了一个计算阶乘的函数st_mul()。具体代码实现如下：

double st_mul(int x)
{
    if (x == 1 || x == 0)
        return 1;
    else
        return x * st_mul(x - 1);
}
double funcos(double e, double x)
{
    int flag = 1;    //用于改变加减符号
    double item = 0; //判断最后一项的绝对值小于e
    double value = 0; //用于存放最后的值
    int i = 0;
    //循环实现
    do
    {
        //单独写一个计算阶乘的函数st_mul()
        value = value + pow(x, i) * flag / st_mul(i);
        //pow(x, i)是math库里中计算幂函数的方法
        item = pow(x, i) / st_mul(i);
        flag = -flag;    //改变符号
        i += 2;
    } while (item >= e);

    return value;
}