home *** CD-ROM | disk | FTP | other *** search
/ Power-Programmierung / CD2.mdf / c / library / dos / math / anum / raphson1.c < prev    next >
Encoding:
C/C++ Source or Header  |  1991-08-24  |  872 b   |  57 lines
  1. /* OK */
  2.  
  3. #include <math.h>
  4. #include <conio.h>
  5. #include <stdio.h>
  6.  
  7. #include "anum.h"
  8. #include "sysio.h"
  9.  
  10.  
  11.  
  12. double f1(double x)
  13.  
  14. {    return(cos(x)-x);
  15. }
  16.  
  17.  
  18. double f2(double x)
  19. {    return(-sin(x)-1);
  20. }
  21.  
  22.  
  23.  
  24.  
  25.  
  26.  
  27. void main(void)
  28.  
  29.  
  30. {    double pt0,rac,tol,val,racderiv;
  31.     int iter,maxiter,errcode;
  32.  
  33.     pt0=1e2;
  34.     tol=1.e-10;
  35.     maxiter=1000;
  36.  
  37.  
  38.     puts("Essai de la procedure Newton-Raphson");
  39.     printf("A partir de %lf et  avec une \n"
  40.            "tolérance de %le et un nb max d'itérations de %d\n",
  41.            pt0,tol,maxiter);
  42.  
  43.  
  44.  
  45.  
  46.     newton_raphson(pt0,tol,maxiter,&rac,&val,&racderiv,&iter,
  47.                 &errcode,*f1,*f2);
  48.         SYSMSG(errcode,stderr);
  49.  
  50.     printf("\n\nRésultats:\n----------\n"
  51.            "\n\nracine : %16.12lf\n valeur: %12.8le\n"
  52.            "derivée : %16.12lf\niterations : %d\nerrcode : %d\n",
  53.            rac,val,racderiv,iter,errcode);
  54.  
  55.  
  56. }
  57.