Exercice D'intégration Numérique
Exercice D'intégration Numérique
Exercice D'intégration Numérique
a- On a:
Z 1
I= f (x) dx
0
1
Z
2
Z 1
= f (x) dx + f (x) dx
1
0 2
1
Z
2
Z 1
2
= x 1 dx + ex sin (x) dx
1
0 2
12 1
x3
1 x
= x + e (sin (x) cos (x))
3 0 2 1
2
0.279222
b- Simpson composee:
Puisque h = 0.05 alors n = 20
Utilisons python pour effectuer calcule:
def f(x):
import math as m
if 0<=x<=0.5:
return x*x-1
else:
return m.exp(x)*m.sin(x)
def Simpson(a,b,n):
h=(b-a)/float(n)
I=(f(a)+f(b))/6
for i in range(1,n) :
I=I+f(a+i*h)/3
for i in range(n) :
I=I+f(a+(2*i+1)*h/2)*2/3
return h*I
>>> Simpson(0,1,20)
0.26638531822189726
>>> Simpson(0,1,25)
0.258600640962109
>>> Simpson(0,0.5,20)+Simpson(0.5,1,20)
0.2728038204032314
Qui est un valeure mieux approchee de lintegrale que celle en question b
e-La meilleure methode des 3 decrites precedement est la derni`ere car son erreur est de lordre de 103
tandis que les erreur des deux autres est de lordre de 102
a- Dapr`es les hypoth`eses la famille des Polynomes est orthogonale, alors reste `a montrer que i {0, .., n}
deg (Li ) = i
On a deg (L0 ) = 0 et deg (L1 ) = 1. Soit i {2, .., n} et supposons que deg (Li1 ) = i 1 et que
deg (Pi2 ) = i 2.
Puisque:
li (x) = (2i x 1) li1 (x) (i 1)2 li2 (x)
Alors:
deg (li ) = deg (2i X 1) li1 (i 1)2 li2
= deg ((2i X 1) li1 )
= deg ((2i X 1)) + deg (li1 )
=1+i1
=i
b- La quadrature en question est une quadrature de Gauss, ainsi le degre de precision est :
2deg (Ln ) 1 = 2n 1
c- Cas n=2
les racines:
On a
L2 (x) = (3 x) l1 (x) l0 (x) = x2 4x + 2
Alors: x1 = 2 2 et x2 = 2 + 2
Les coefficients:
On a dapr`es ce qui prec`ede le degres de precision est 2n-1, on en deduit les equations suivantes:
Z +
a1 + a2 = ex dx = 1
0
Z +
a1 x1 + a2 x2 = xex dx = 1
0
Ce qui donne:
2+1 21
a1 = = 0, 85355 et a2 = = 0, 14644
2 2 2 2
d- Cas n=2, on:
+
Z
xex dx a1 x1 + a2 x2 = 0, 9238
0
Cas n=3, on a:
+
Z
xex dx a1 x1 + a2 x2 + a3 x3 = 0, 90644
0