xdnewton.h

c源码

/* x.h                    freeware                  [e-mail]  */

/* ------------------------------------- FUNCTION --------------------------  */
/* Do    :                                                                    */
/*                                                                            */
/* -------------------------------------------------------------------------- */
double Newton_s_Method(
double x,
   int imax,
double (*P_f)(double x),
double (*PDf)(double x)
)
{
   int i;

       for(i=1; i<imax; i++)
          {
           x = x - ((*P_f)(x)) / ((*PDf)(x));
          }
 return(x);
}

/* ------------------------------------- FUNCTION --------------------------  */
/* Do    :                                                                    */
/*                                                                            */
/* -------------------------------------------------------------------------- */
double p_Newton_s_Method(
double x,
   int imax,
double (*P_f)(double x),
double (*PDf)(double x)
)
{
   int i;

       for(i=1; i<(imax+1); i++)
          {
           printf(" x[%d] = %.15f\n",i,x);
           x = x - ((*P_f)(x)) / ((*PDf)(x));
          }
 return(x);
}