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