0

単一の関数ポインターを使用して Newton-Raphson メソッドを実装しようとしています。関数には、方程式とその導関数の両方が含まれている必要があります。テスト関数内でこれら 2 つの別個の関数を通過するのが困難です。

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

double NR(double, double(*)(double, double*), double);

void test_function( double x, double * f, double * f_prime )
{
  *f = (x-2) * (x-2);
    *f_prime = 2*x - 4;
}

double NR( double x0, double (*test_function)(double x, double *f, double *f_prime), double precision )
{
  int i;
    while(!isnan(x0)){
      i = x0;
      x0 = (x0 - (test_function(x0, f, 0)/test_function(x0, 0, f_prime)));
      if(!isnan(x0))
    printf("%f\n",x0);
      if ( i - x0 < 0 )
    printf("NO ROOT FOUND");
      return -1;
      else if ( i - x0 > 0 && i - x0 < precision )
    break;
    }

}


int main(void)
{
  double x0 = 300;
  double precision = .0000001;

  double root = NR( x0, test_function, precision);

  printf("%f\n",root);

  return 0;

}

ありがとうございました

4

1 に答える 1

0

NR3 つのパラメーターを受け取る関数ポインターを受け取るように宣言を更新します。

double NR(double, double(*)(double, double*, double *), double);

それ以外の

double NR(double, double(*)(double, double*), double);

0(またはNULL) をprimeまたはに渡すとf_prime、テスト関数を次のように更新できます。

void test_function( double x, double * f, double * f_prime )
{
    if(f)
      *f = (x-2) * (x-2);
    if(f_prime)
      *f_prime = 2*x - 4;
}

NRまた、以下のように更新して定義する必要がありfます。f_prime

double NR( double x0, double (*test_function)(double x, double *f, double *f_prime), double precision )
{
  int i;
  double f, f_prime;
  ...
  //your code
  ...
  //----------------------------v pass address -------------v
  x0 = (x0 - (test_function(x0, &f, 0)/test_function(x0, 0, &f_prime)));
  ....
  //your code

}
于 2013-11-11T18:55:07.517 に答える