Aller au contenu

Mathc initiation/Fichiers c : c71c03

Un livre de Wikilivres.


Sommaire


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

c01c.c
/* ---------------------------------- */
/* save as c1c.c                      */
/* ---------------------------------- */
#include "x_hfile.h"
#include      "fc.h"
/* ---------------------------------- */
int main(void)
{
int      n =  2*50;
double   a =  0.;
double   b =  .5;

 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 = 100)

    (0.500
 int(      (asin(x)) * dx = 0.127825
    (0.000



 With the antiderivative of f.

 F(x) = asin(x)*x + sqrt(1 - x**2) 


 F(0.500) -  F(0.000)  = 0.127825 


 Press return to continue.



Calculons la primitive :
                           /
Calculer la primitive de  | asin(x) dx 
                          /
       
       Utilisons l'intégration par partie
       
        u = ...    dv = ...
       du = ...     v = ...
       
        /              /
       | u dv = u v - | v du
       /              /
       
       Nous savons que la dérivé de asin(x) est 1/sqrt(1-x**2).
       Nous ne connaissons pas sa primitive. 
       Nous allons multiplier asin(x) par 1. 
       Nous allons poser que u = asin(x) et dv = 1 dx.
       
       
        /               
       | (asin(x)*1) dx = 
       /               
       

        u = asin(x)               dv = 1 dx
       du = 1/sqrt(1-x**2) dx      v = x   
       

        /                        (u*v)    /(v                   du)
       | (asin(x)*1) dx =  asin(x) * x - | (x * 1/sqrt(1-x**2)) dx
       /                                 /
 

        /                          /
       | asin(x) dx = asin(x) x - | 1/sqrt(1-x**2)  x dx
       /                          /   
      
                         _______________________
                        |         u = (1-x**2)  |
                        |        du =   -2*x dx |
                        | (-1/2) du =      x dx |
                        |_______________________|
          

        /                           /
       | asin(x) dx = asin(x) x -  | 1/sqrt(u)  (-1/2) du
       /                           /   
       

        /                                /
       | asin(x) dx = asin(x) x + (1/2) | u**(-1/2)  du
       /                                /


        /                                   
       | asin(x) dx = asin(x) x + (1/2) (2/1) u**(1/2) + C
       /                                   
       
        /                                   
       | asin(x) dx = asin(x) x +  sqrt(1-x**2) + C
       /