Как реализовать метод деления пополам, используя этот прототип функции в C ++ - PullRequest
Новая
       8

Как реализовать метод деления пополам, используя этот прототип функции в C ++

1 голос
/ 20 марта 2011 
double p1::root(double (*pf)(double k), int a, int b, double e)

Я не уверен, как это сделать, я понимаю, что мне нужно цикл, который точно определяет середину и все такое 

double p1::root(double (*pf)(double k), int a, int b, double e) {


  // void nrerror(char error_text[]);                                           
  int j;
  float dx, f, fmid, xmid, rtb;

  f = (*pf)(a);
  fmid = (*pf)(b);
  //if (f*fmid >= 0.0) nrerror("root must be bracketed for bisection in rtbis")\
;                                                                               
      rtb = f < 0.0 ? (dx=b-a,a) : (dx=a-b,b);
      for(j = 1;j <40; j++) {
        fmid = (*pf)(xmid = rtb+(dx *= .5));
      if (fmid <= 0.0) rtb = xmid;
      if (fabs(dx) < e || fmid == 0.0) return rtb;
      }
      // nrerror("too many bisections in rtbis");                               


  return 0.0;
}




double p1::test_function(double k) {
  return (pow(k, 3) -2);
}

тогда в основном у меня есть это 

 double (*pf)(double k);
pf = &p1::test_function;

//double result = p1::root(pf, a, b, e);

1 Ответ

0 голосов
/ 20 марта 2011

Может быть Числовые рецепты даст вам идею.Подсказка: это рекурсивно.

...