因此,我正在尝试创建一个使用泰勒近似计算 cos(x) 的程序。
该程序非常简单:用户输入一个参数 x(x 是一个以弧度为单位的角度)和一个浮点 ε,它是 cos(x) 值的精度。
基本上,程序唯一要做的就是计算这个总和:x^0/0!- x^2/2!+ x^4/4!- x^6!+ x^8/8!- ...,直到项小于 ε,即 cos(x) 的值,它将在我们的精度范围内。
所以,这里是代码:
#include <stdio.h>
/* Calculates cos(x) by using a Taylor approximation:
cos(x) = x^0/(0!) - x^2/(2!) + x^4/(4!) - x^6/(6!) + x^8/(8!) - ... */
int main(void)
{
int k; // dummy variable k
float x, // parameter of cos(x), in radians
epsilon; // precision of cos(x) (cos = sum ± epsilon)
sum, // sum of the terms of the polynomial series
term; // variable that stores each term of the summation
scanf("%f %f", &x, &epsilon);
sum = term = 1, k = 0;
while (term >= epsilon && -term <= epsilon)
// while abs(term) is smaller than epsilon
{
k += 2;
term *= -(x*x)/(k*(k-1));
sum += term;
}
printf("cos(%f) = %f\n", x, sum);
return 0;
}
起初,我试图通过计算单独变量“fact”的阶乘来解决它,尽管即使 ε 具有合理的大值也会导致溢出。
为了解决这个问题,我注意到我可以将前一项乘以 -x² / (k(k-1)),在每次迭代中将 k 增加 2,以获得下一项。我以为这会解决我的问题,但话又说回来,它不起作用。
该程序编译得很好,但是例如,如果我输入:
3.141593 0.001
输出是:
cos(3.141593) = -3.934803
……这显然是错误的。有人能帮我吗?