LLOCLLE BELITO JOSE LUIS 20191257F
2. Cuadratura Gaussiana
( a ) Encuentre el polinomio de legendre Ln (x) para n=4 por la formula de recurrencia .
Calcule sus raíces de forma algebraica y estime el valor de r 1 , r 2 , r 3 ,r 4 con 8 cifras decimales .
Fórmula de recurrencia :
2n+1 n
Ln+1 ( x )= x ln ( x )− L ( x)
n+ 1 n+1 n−1
1
L0 ( x ) =a : ∫ a2 dx=2 → a=1
−1
L0 ( x ) =1
1
2
L1 ( x )=bx+ c : ∫ ( bx +c )2 dx=
−1 3
1 1
∫ L1 ( x ) L0 ( x ) dx=∫ ( bx +c ) dx=0 → b=1 , c=0 , L1 ( x )=x
−1 −1
3 1 3 x 2−1
L2 ( x ) = x L1 ( x ) − L0 ( x ) =
2 2 2
5 2 5 x 3−3 x
L3 ( x ) = x L 2 ( x ) − L1 ( x ) =
3 3 2
7 3 35 x 4 −30 x 2+3
L4 ( x ) = x L 3 ( x ) − L2 ( x ) =
4 4 8
35 x 4 −30 x 2+3
L4 x =
( ) =0
8
35 x 4−30 x 2 +3=0
15+2 √ 30
r 1=
√ 35
=0.86113631
15+2 √ 30
r 2=−
√ 35
=−0.86113631
15−2 √ 30
r 3=
√ 35
=0.33998104
15−2 √30
r 4 =−
√ 35
=−0.33998104
(b)Use las formulas de cuadratura Gauss−Legendre :
2
W i=
(1−r 2i )¿ ¿
para calcular loscoeficientes w1 , w2 , w3 , w 4 en :
1
∫ f ( x ) dx=w 1 . f ( r 1) + w2 . f ( r 2 ) +w 3 . f ( r 3 ) +w 4 . f ( r 4 )
−1
2
W 1=
(1−r 21 )¿ ¿
2
W 2=
(1−r 22 )¿ ¿
2
W 3=
(1−r 23 )¿ ¿
2
W 4=
(1−r 24 )¿ ¿
(c )Calcule de forma exacta usando lasreglas de integración el valor de la siguiente integral:
5
x −1
∫ dx
1 x ( 1+ x 2 )
5
5
)|
2
ln ( x +1 )
∫ x −1 2 dx= −ln| x|+ 2 + arctan(x )
(
1 x ( 1+ x ) 1
5
x −1 ln|26| 2 ln |2|+ π
∫ dx= −ln|5|+arctan ( 5 ) −
1
2
x ( 1+ x ) 2 4
5
x −1
∫ dx=0.2610393698442355
1 x ( 1+ x 2 )
(d ) Implemente un algoritmo de Matlab que permite calcular el valor de la integral en(c)usando la
formulade traslación en la formula de(b).
(e )Calcule el error absoluto para los valores real y aproximado obtenidos en( c) y (d )respectivamente .
ERROR=|0.26103936−0.26256582|=0.00152646