Mathc initiation/Fichiers c : c73c01

Un livre de Wikilivres.


Sommaire


Installer et compiler ces fichiers dans votre répertoire de travail.

c01a.c
/* ---------------------------------- */
/* save as c1a.c                      */
/* ---------------------------------- */
#include "x_hfile.h"
#include      "fa.h"
/* ---------------------------------- */
int main(void)
{
int      n =  2*250;
double   a =  0.;
double   b =  3.;

 clrscrn();

 printf(" With the Simpson's rule.    (n = %d)\n\n"
        "    (%.3f\n"
        " int(      (%s)  dx = %.6f\n"
        "    (%.3f\n\n\n\n",n,  b, feq, simpson(f,a,b,n), a);

 printf(" With the antiderivative of f.\n\n"
        " F(x) = %s \n\n\n" 
        " F(%.3f) -  F(%.3f)  = %.6f \n\n\n", Feq, b,a, F(b)-F(a));
 
 stop();

 return 0;
}
/* ---------------------------------- */


Calculons l'intégrale avec la fonction simpson(f,a,b,n); puis avec sa primitive F(x).


Exemple de sortie écran :
 With the Simpson's rule.    (n = 500)

    (3.000
 int(      (x**3 * exp(2*x))  dx = 3479.948346
    (0.000



 With the antiderivative of f.

 F(x) = x**3 exp(2*x)/2- 3 x**2 exp(2*x)/4+ 6 x exp(2*x)/8- 6 exp(2*x)/16 


 F(3.000) -  F(0.000)  = 3479.948344 


 Press return to continue.



Calculons la primitive :
             
Calculer la primitive de 

       
        /               
       | x**3 * exp(2*x)  dx = 
       /               
    __________________        _______________
    |Signe|  Diff     |       |   Int       |
    |_____|___________|       |  exp(2*x)   |--> Attention petit décalage 
    |  +  |     x**3  |------>|  exp(2*x)/2 |    d'une ligne vers le haut   
    |  -  |   3 x**2  |------>|  exp(2*x)/4 |    pour éviter les fèches       
    |  +  |   6 x     |------>|  exp(2*x)/8 |    en diagonales.          
    |  -  |   6       |------>|  exp(2*x)/16|  
    |  +  |   0       |------>|             |       
    |_____|___________|       |_____________|         
            
            
            
        /               
       | x**3 * exp(2*x)  dx = +   x**3 exp(2*x)/2 
       /                       - 3 x**2 exp(2*x)/4
                               + 6 x    exp(2*x)/8 
                               - 6      exp(2*x)/16