Mathc initiation/Fichiers c : c79c04
Installer et compiler ces fichiers dans votre répertoire de travail.
|
c01d.c |
|---|
/* ---------------------------------- */
/* save as c1d.c */
/* ---------------------------------- */
#include "x_hfile.h"
#include "fd.h"
/* ---------------------------------- */
int main(void)
{
int n = 2*50;
double a = .5;
double b = 1.;
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)
(1.000
int( (1/(1-sin(x)) dx = 1.722427
(0.500
With the antiderivative of f.
F(x) = tan(x) + sec(x)
F(1.000) - F(0.500) = 1.722427
Press return to continue.
Calculons la primitive :
/
| 1
1) Calculer la primitive de | -------- dx
| 1-sin(x)
/
/ /
| 1 | 1 1+sin(x)
| -------- dx = | -------- -------- dx
| 1-sin(x) | 1-sin(x) 1+sin(x)
/ /
/
| 1+sin(x)
= | ------------ dx
| 1-sin(x)**2
/
sin(x)**2+cos(x)**2 = 1
/
| 1+sin(x)
= | --------- dx
| cos(x)**2
/
/ /
| 1 | sin(x)
= | --------- dx + | --------- dx
| cos(x)**2 | cos(x)**2
/ /
/ /
= | sec(x)**2 dx + | tan(x) sec(x) dx
/ /
/
| 1/(1-sin(x)) dx = tan(x) + sec(x) + c
/