Предыдущая страница реферата | 1  2  3  4  5 | Следующая страница реферата

/* Bracketing function''s root. The function is supposed to have unlimited

   domain and be continuous.

   int BracketRoot(double x0,double *a,double *b,double d0,

                    double di, double dmax,double (*fun)(double));

   Parameters:

    x0 - initial guess on input;

    a  - left bound on output;

    b  - right bound on output;

    d0 - initial interval of hunting;

    di - interval increment (geometric progression multiplier);

    dmax - maximal interval;

    fun - pointer to the function.

   Returns:

    1 - if a root is bracketed;

    0 - on failure

*/

int BracketRoot(double x0,double *a,double *b,double d0,

                 double di, double dmax,double (*fun)(double)) {

  double fa,fb,f0;

  /* get initial function guess, initial a,b,fa,fb */

  f0=(*fun)(x0); *a=x0-d0; *b=x0+d0; fa=(*fun)(*a); fb=(*fun)(*b);

  /* while the increased search interval is less than maximal,

     process cycle */

  while((d0*=di)<dmax) {

    /* check up the bracketing success. Case f0>0. */

Рекомендуем скачать другие рефераты по теме: курсовая работа производство, зимнее сочинение.



Предыдущая страница реферата | 1  2  3  4  5 | Следующая страница реферата

Поделитесь этой записью или добавьте в закладки