Mathc initiation/Fichiers c : c73c01
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