1

シンプソンズ法と台形法を使って関数の積分を計算するプログラムを作らなければなりません。MinGW を使用しているコンピューターでは問題なく動作しますが、大学のコンピューターでコンパイルしようとすると、次のようになります。

It = 0 //should be 0.954499
*pn = 0 //should be 18
Is = 0 //should be 0.954500
*pn = 0 //should be 6

これが私が思いついたものです（変数とコメントはポルトガル語で申し訳ありません。家に帰ってから修正します）：

積分.h:

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

#define eps 1.e-6
#define kmax 20

double trapezio(double (*f)(double x), double a, double b, int *n);
double simpson(double (*f)(double x), double a, double b, int *n);

積分.c:

#include "integral.h"
#define xmin -2
#define xmax 2

double f(double x);

int main(){
    double It,Is;
    int n = 0;
    int *pn = NULL;

    pn = &n;
    It = trapezio(f,xmin,xmax,pn)/sqrt(2*M_PI);
    printf("Pelo metodo dos trapezios a integral vale aproximadamente %lf\n", It);
    printf("O numero de iteracoes usadas foi %d\n\n",*pn);

    *pn = 0;
    Is = simpson(f,xmin,xmax,pn)/sqrt(2*M_PI);
    printf("Pelo metodo de simpson a integral vale aproximadamente %lf\n", Is);
    printf("O numero de iteracoes usadas foi %d\n",*pn);

    return 0;
}

double f(double x){
    return exp(-0.5*x*x); // Funcao que sera integrada
}

トラペジオ.c:

#include "integral.h"

double trapezio(double (*f)(double x), double a, double b, int *n){
    double To, Tk;
    double soma;
    int i, k = 1;
    Tk = 0.5*(f(a) - f(b))*(b - a);

    while (fabs((Tk-To)/To) > eps && k < kmax){
        soma = 0; // Resetando variavel soma
        To = Tk; // To e' T(k - 1), caso o loop se repita o ultimo Tk vira To
        for (i = 1 ; i <= (pow(2,k)-1) ; i += 2) soma += f(a + i*(b - a)/pow(2.,k));
        Tk = 0.5*To + soma*(b - a)/pow(2.,k);
        k++;
        *n += 1;
    }

    return Tk;
}

simpson.c:

#include "integral.h"

double simpson(double (*f)(double x), double a, double b, int *n){
    double So, Sk = 0;
    double somaimp, somapar;
    int i, k = 1;

    while (fabs((Sk-So)/So) > eps && k < kmax){
        somaimp = 0;
        somapar = 0;
        So = Sk; // So e' S(k - 1)
        for (i = 1; i <= (pow(2,k)-1); i += 2) somaimp += f(a + i*(b - a)/pow(2.,k));
        for (i = 2; i <= (pow(2,k)-2); i += 2) somapar += f(a + i*(b - a)/pow(2.,k));
        Sk = (b - a)*(f(a) + 4*somaimp + 2*somapar + f(b))/(3*pow(2.,k));
        k++;
        *n += 1;
    }

    return Sk;
}

編集: ポインターを取り出すと、trapezio は機能しますが、simpson は引き続き 0 を返すことを忘れていました。

4

1 に答える 1

6

では、ループで使用する前にtrapezio()初期化することはありません。Towhile

double To, Tk;
/* ... no assignment to To ... */
while (fabs((Tk-To)/To) > eps && k < kmax){

これは、定義されていない方法で動作し、おそらくまったくループに入らないことを意味します。

于 2013-06-19T12:27:10.927 に答える