c - 泰勒定理 C

Question

我有泰勒定理的简单代码cosh()。

我试图抓住一个错误 - 结果有时接近真正的答案。

如何正确执行？

当我的开始为 0，结束为 5，细分为 5 时，它给出了很好的结果，但是当我将 5 作为开始，将 10 作为结束时，结果离期望值更远。

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

int poww( float number, int a )
{
    float result = 1.0;
    int i;
    if( a != 0 );
    {
        for( i = 0; i < a; i++ ) {
            result = result * number;
        }
    }
    return result;
}

int factorial(int n)
{
    switch (n) {
    case 0:
        return 1;
        break;
    default:
        return n * factorial(n-1);
    }
}

void main()
{
    puts("Enter start: ");
    float start;
    scanf("%f", &start);

    puts("\nEnter end: ");
    float end;
    scanf("%f", &end);

    puts("\nSubintervals:");
    int subinterval;
    scanf("%d", &subinterval);

    float h = (end - start) / (float)subinterval;
    printf("h is : %3.2f \n", h);

    double x, result, temp;
    int n;

    for( x = start; x <= end; x += h) {
        result = 0;
        for(n = 0 ;  ; n++) {
            temp = poww(x, 2 * n) / (factorial( 2 * n ) * 1.0);
            if(temp < 0.00001) {
                break;
            } else {
                result = result + temp;
                printf("X = %f temp = %f, result = %f\n", x, temp, result);
            }
        }
        printf("X = %f, result = %3.2f, cosH = %3.2f\n\n", x, result, cosh(x) );
    }
    puts("Press any key...");
    getchar();
}

问题解决：

函数返回一个整数而不是双精度，我也将每个浮点数更改为双精度。

Answer 1

也将所有float类型更改为double并用作and函数double的返回类型。在这种情况下，最后两个是最重要的。factorial()poww()

此外， on 的返回类型main()应该是int，而不是void。

[我刚刚删除了 poww() 中的死 if 语句，并注意到该函数仅“加速”了 pow() 计算。如果您担心性能，请担心计算每个项的阶乘和幂，而不是将前一项乘以 x^2 并除以 (2*n)*(2*n-1)。]

在您的代码的这个小修复中，我在 4 到 10 之间得到了很好的结果：

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

double poww( float number, int a )
{
    float result = 1.0;
    int i;
    for( i = 0; i < a; i++ )
    {
        result = result * number;
    }
    return result;
}

double factorial(int n)
{
     switch (n)
      {
            case 0: return 1;
                    break;
            default: return n * factorial(n-1);
      }
}

int main(){

    puts("Enter start: ");
    float start;
    scanf("%f", &start);

    puts("\nEnter end: ");
    float end;
    scanf("%f", &end);

    puts("\nSubintervals:");
    int subinterval;
    scanf("%d", &subinterval);

    float h = (end - start) / (float)subinterval;
    printf("h is : %3.2f \n", h);

    double x, result, temp;
    int n;

    for( x = start; x <= end; x += h){
        result = 0;
        for(n = 0 ;  ; n++){
            temp = poww(x, 2 * n) / (factorial( 2 * n ) * 1.0);
            if(temp < 0.00001){
                break; }
            else{
                result = result + temp;
                printf("X = %f temp = %f, result = %f\n", x, temp, result);
            }
        }
                printf("X = %f, result = %3.2f, cosH = %3.2f\n\n", x, result, cosh(x) );
    }

    puts("Press any key...");
    getchar();
    return 0;
}