用泰勒级数求C中的sin（x）

Question

我有一个作业，教授要我们用泰勒级数计算sin（x）。他希望我们在两个连续分数之间的差小于10^-6时停止迭代。

最后，我说，例如x^5/5！与（x^3/3！）*相同（x^2/4*5），所有分数都是如此。所以我可以保留之前计算的分数，并在下一次迭代中使用。问题是，我最终得到的数字与它的实际罪过有点偏差，我不知道为什么。提前谢谢。这是我的代码：

#include <stdio.h>
#include <Math.h>
#define pi 3.14159265358979323846

int main(int argc, int **argv){
    
    int sign = -1, pwr = 3;
    double previous, current, rad,sum, degr;
    printf("Calculating sin using Taylor Series\n\n");
    printf("Give degrees: ");
    scanf("%lf", &degr);
    
    // translate to rads
    rad = degr*(pi/180);
    
    sum = rad;
    previous = rad;
    do{
        current = (previous * pow(rad, 2))/(pwr* pwr-1);
        
        sum += sign*current;
        pwr += 2;
        sign *= -1;
    }
    while(abs(current - previous) > pow(10, -6));
    
    printf("The sin of %lf degrees is ", degr);
    printf("%.6f\n", sum);
    printf("%.6f", sin(rad));
    return 0;   
}

Answer 1

您正在使用< code>abs函数，该函数需要一个< code>int并返回一个< code>int。如果当前项和先前项之间的差小于1，这会导致循环存在，因为它会将< code>diff设置为0。

相反，你想要期待并返回双精度的晶圆厂。