Accueil | Cast3M

   1 : $$$$ FZERO    NOTICE  BP208322  19/02/22    21:15:01     10118          
   2 :                                              DATE     19/02/22
   3 :                                              
   4 :   Procedure FZERO                     Voir aussi : ZLEG RACP 
   5 :     ---------------  
   6 :     XSOL = FZERO FONCTION XA XB (XTOL);
   7 : 
   8 : 
   9 :     Objet :
  10 :     _______
  11 : 
  12 :     Cette procedure calcule le zero d'une fonction x->f(x) elle-meme
  13 :     definie par une procedure.
  14 :     La methode de Brent est utilisee.
  15 : 
  16 :     Commentaire :
  17 :     _____________
  18 : 
  19 :     FONCTION : PROCEDURE definisant la fonction x -> y=f(x) telle que :
  20 :                y = FONCTION x;
  21 :                
  22 :     [XA XB]  : intervalle de recherche (FLOTTANTs)
  23 :     
  24 :     XTOL (facultatif) : precision sur la solution (FLOTTANT)
  25 :     
  26 :     Exemple :
  27 :     _________
  28 :     debp f6 x;
  29 :       y  = (exp  x) - 2. - (1. / ((10.*x)**2)) + (2./((100.*x)**3));
  30 :     finp y;
  31 :     
  32 :     x6 = FZERO f6 -4. 2.;
  33 :     y6 = f6 x6;
  34 :     * the result is : x6=0.70320 et f6=-2.7E-16 which is correct !
  35 :     
  36 :     * graphical verification :
  37 :     x6_p = prog -4. PAS 0.01 2.;
  38 :     f6_p = f6 x6_p;
  39 :     ev6 = evol bleu manu 'x' x6_p 'f(x)' f6_p;
  40 :     dess ev6 'YBOR' -5 5 'YGRA' 1 'AXES';
  41 : 
  42 :     
  43 :