INTEGRACION

Y DERIVACION
NUMERICAS

1. Aplíquese la regladel trapecio en un sólo intervalo para aproximar el va-

lor de la integral )r{2" - 1) d-x . Compárese con el valor exacto.

2. Aplíquese la regla del trapecio en un sólo interr¿alo para aproximar el va-

,z( I\
lor de la integral J. - *t . Comnárese con el valor exacto.
Ir" "' Jrlx

Compruébese cómo mejora la estimación al aplicar la regla del trapecio

en dos. cinco v diez subintervalos.

3. -- integrat
-- la
Esrímese el valor de -..'-o- -- l" .5i + x + 4 r/x aplicando la regla
lo del
3.tr. -
trapecio en el intervalo [0, 2]

4. Divídase el intervalo [0,21r) en doce partes iguales para estimar el valor

de la siguiente integral, utilizando la regla del trapecio

12n 1,
I
JO
2+cosx

161
3. Integración y derivación numéricas

5. Aplíquese la fórmula de integración numérica por interpolación en dos

puntos y compruébese que se obtiene la regla del trapecio.

6. Aplíquese la fórmula de integración numérica por interpolación en tres

puntos y compruébese que se obtiene la regla de Simpson.

Aplíquese la regla de Simpson para estimar el valor de la integral

¡5
f. x(-r- 4)= dx y compárese el resultado obtenido con el valor exacto.

8. Estímese el valor de la integral lu ¿^' dx mediante la regla de Simpson

d{

en tres subintervalos (esto es, con 7 nodos), tealizando las operaciones

con una calculadora científica.

9. Desarróllese un programa informático para estimar el valor de la integral

e6
I' e-r" dx mediante la regla de Simpson en n subintervalos, haciendo
J4

vaiar n desde t hasta 20.

10. Aplíquese la regla de Simpson en 1, 10, 20,30,..., 100 subintervalos p¿Ira

aproximar el valor de la integral

l"-Lar.
Jo 4x'+1296

Compiárese con el valor exacto.

11. Aplíquese la fórmula de integración numérica por interpolación en cuatro

puntos igualmente espaciados para obtener un regla de integración numé-
rica (regla 3/8).

12. Compruébense analíticamente los siguientes resultados, que aparecen en

la deducción de la fórmula del error para la regla de Simpson

f" A -roX,x - r'Xx- xr) dx = o

162
 Enunciados

f. U"- roX¡ - x,)' (x - x')l dx = #

en donde o 1 b y,como es habitual, xo = cl, xr:T' x, = b'

13. AplíqueselaregladeSimpsonenunsolointervaloparaestimarelvalor
dos fór-
dela integral sen x dx. Acote el error cometido mediante las
f
acotacio-
mulas del error para la regla de Simpson. compárense estas dos
calculando el valor de la
nes con el error exacto, que podemos conocer
integral sin aProximaciones.

la regla del
t4. Determínese, mediante la fórmula de acotación del error en
estimar por ese procedi-
trapecio, el número de intervalos necesario para
miento el valor de la integral

18
11, ¡1 ¿y

con un error de truncamiento menor que una centesrma'

15. Determínese, mediante la fórmula de acotación del error en

la regla de
Simpson, el número de intervalos necesario para estimar
por ese procedi-
miento el valor de la integral

l"-Ld,
ro 4x'+1296

con un error de truncamiento menor que una centesrma'

16. para estimar el valor de la integral

I* ,-" dx, se aplica la regla de Simpson
en el intervalo [-3, 3], dividido en cinco subintervalos iguales.
Determí-
nese una cota del error cometido.

163
3. Integracíón y derivación numéricas

17. LasreglasdeltrapecioydeSimpsonsehandeducidointegrandounpoli-
nomio"interpoladordegradosunoydosrespectivam€nte.Deduciruna
('t) dx, interpo-
nueva regla de integración numérica, para aproximar J f
exponencial de la for-
lando entre los extremos del intervalo una función
ma q(x) = Aeo'.

nodos'
18. Deducir la fórmula de cuadratura de Gauss con tres

19. Estímese el valor de la integrul Jt e-" dx mediante la fórmula de cua-

dratura de Gauss con 2 nodos'

20. Estímese el valor de la integral l-e -" dx

J-t
mediante la fórmula de cua-

dratura de Gauss con 3 nodos'

numéricamente la integral
21. Desarrolle un programa informático para calcular

l"-Ld,
Jo "Jx'+1296

mediantelafórmuladecuadraturadeGausscondosnodos,dividiendo
previamenteelintervaloenl'10,20'30'"''l00subintervalos'Compá-
10'
iense los resultados con los obtenidos en ejercicio

numéricamente la integral
22. Desarrolle un programa informático para calcular
5

f' -J-
rr l+x
¿"

mediantelatórmuladecuadraturadeGausscontresnodosencadasub-
intervalo,considerando2,4,6,8yl0subintervalos.Compáresecadare-
sultadoconelqueseobtienealconsiderarelmismonúmerototaldeno-
dos'peroaplicandolafórmuladecuadraturadeGausscondosnodosen
cadasubintervalo'Utilíceseparaelloelprogramadesarrolladoenelejer-
cicio 21'

164
 Sotruciones

divididas centradas con tamano

23. Aplíquese la fórmula de las diferencias
de la derivadaf'(0,25), en donde/
de paso /z = 0,01 para-ástimar el valor
es la función f lx) = !;r' Compárese con
el valor exacto'

24. Compruébesequelafórmuladelasdiferenciasdivididascentradaspro-
cuando se aplica a polino-
porciona resultados Á",o' para la derivada
mios de grado menor que tres'

)< Estímese el valor de la derivada en el

punto r = 0'9 de la función

f(x) =rím (1+-r) (t.;) [t.;)'

utilizandolafórmuladelasdiferenciascentradasdivididasconlncre-
mentoh=0,0|.Estúdiesecondetenimientocómoevaluarlafunciónl.pa- I
esta deri-
posible calcular el valor de
ra obtener la estimacián pedida. ¿Es
,tI

vada de forma exacta'| i$

tl

función/(x) = x e'co-s 1-en el pun- $

il
26. Estímese el valor de la derivada de la
to x Tt,interpolando con tres nodos
y tamaño de paso h = 0'0I ' t
=
de la funciónf r,x) = x ¿' cos x en
27. Estímese el varor de la derivada segunda
el punto x= Tr,interpolando i'"' nodos y tamaño de paso l¿ = 0'01'
"o"

Soluciones

para aproximar el va- I

1. Aplíquese la regla del trapecio en un sólo intervalo

el valor exacto'
lor de la integral lit.z* -1) r/"r ' Compárese con

Solución. Sea/(x) -2x- l' Porlaregla deltrapecio

dx=(4- I)
r(l)t f(4) =31!r7 =12'
f" ftrl
JrL-

165
3. lntegración y derivación numéricas

Es sencillo calcular en este caso el valor

exacto

f ft"l clx=Íx' - xl',=i : 12'

Noescasualqueelerrorresultesernulo.Laregladeltrapecioprocede
1 en vez del integrando
de integrar el polinomio interpolador de orden
dado.Peroelpolinomiointerpolado'deprimerordencoincideconelin-
grado l
tegrando .uundo este es un polinomio de '

para aproximar el va-

2. Aplíquese la regla deltrapecio en un sólo intervalo
.z( II .
Compárese con el valor exacto'
lor de la integral f ['" - 'x' +
;)O-
Compruébesecómomejoralaestimaciónalaplicarlaregladeltrapecio
en dos, cinco Y diez subintervalos'

Solución.Aplicamoslaregladeltrapecioenelintervalo[0,2]alinte-
grando .f (x) =2x- x' +I

,,(- , t) ,.. "f(VlQ\ =l*l=1.

l"lz"-x'+;ldr=t-
uu
\ z) z 2'2

Enestecasopodemoscalcularfacilmenteelvalorexactodelaintegral

,,( r\ I ,*,I l'='-7

[, -1" = j -t ??
J, - x' +t)ur= "1.=o
¡t" L"'

166
 Soluciones

,-,-t:;-ii-:"'a.

Figura 3'1' .f (x)=2x- xt +'; e interpolación lineal enl'2y 5

subintervalos.

4 1 : t,^ ^^- ^l valor nrra se

cr qulere
en relación con el "^l^- que
El error cometido, 1, "t grande
que la gráfica def una.parábola'
aproximar. En la frgura 3.1 podemos ver
*uy curvada en el iniervalo [0, 2]' alejándose bastante de la recta
"itá que interpola linealmente en ese intervalo'
y= |
Apliquemos ahora la regla del trapecio en dos subintervalos

,z( l) z(f(o)+f(2).,,,..)-1,1-',
I, -,="*;)ax= ,\-f * ¡rt)
):tnt='
¡t,

167
3. lntegración y derivación numéricas

Ahora dividimos el intervalo en cinco paftes iguales

¡z ( -,, )ld^=;l-
lZ"--r'*;
|.,__?l f\ot+ f(z)*i r f,jll= !=2,2s
l" L) ,\ - Á" ''') 25
y, finalmente, en diez Partes

l' ( r,-"' *
r
=
)'
z ( rot+ rg)*i
t)o' ro[ 2 -k' r rat)-
t1+il=
53:2,t2
x
,o
\
En la figura 3.1 pueden verse las aproximaciones de/correspondientes a
la interpolación lineal en 1,2 y 5 subintervalos.

3. Esrímese el valor de la integra, f t** dx aplicando la regla del

trapecio en el intervalo [0, 2]

Solución. Siendo .f (x)=# integrando, tenemos

*^el

Figura 3.2. .f (x):#"* y el polinomio interpolador de grado 1.

168
 Soluciones

t' _ _ !__ ¿*=2.f Q)+.f (2) = 1 *_L = 1= 0,78.

Jn 3x' +x+4 2 4 102 68
Utilizando métodos más eficaces que la regla del trapecio en un sólo in-
tervalo, puede comprobarse que el error cometido en la anterior estima-
ción es menor que 0.02. En la figura 3 podemos comprobar que en la pri-
mera mitad del intervalo se comete un error por defecto al reemplazar la
función/por el polinomio interpolador de grado 1, mientras que en la se-
gunda mitad del intervalo, el error es por exceso; obviamente, al integrar,
uno se compensa con el otro parcialmente.
4. Divídase el intervalo [0, 2n] en doce pafes iguales para estimar el valor
de la siguiente integral, utilizando la regla del trapecio

l'" |
Jo 2+cos"t ,r.r.

Solución. Escogemos los subinterr¿alos: [0, :], Ít,:1,...,[Y ,zrc].

Aplicando la regla del trapecio, dividiendo el recinto de integración en 12
partes iguales
/ \
f" '\¿v'v"'
f(*td*=Ll ftol+ftzrt*i rf ':ll
ro 61. 2 E''"')
en donde/es el integrando, .f ( x) : 1
. Por una parte
2 *.o, ,

/(o) =
+,,(z)= ;f,t,,?+)=?, t(T)= ),,(T)=?,
s"\ z
'\"( 6 / a-./¡
.f(n)=I , \
f etTt =!. fl 11o l =
2 .( tott\ 2 "( gft\ | "lsz) 2
3 \6/ 4+8 tl. u J:r't[?.,l=t'r[?,j=t'
l1*\ .L
.l ttt I
( 6 / a-^lt

169
3. Integración y derivaeión numéricas

^(kr\ , 62
luego ) fl l--!-
ii, ./ lor-lr'rs' Finalmente

l" f(rt ¿, =L( !*Z*9]

3,62i5ee
Jo 6[3 13 t5) =

Esta integral puede calcularse exactamente utilizando procedimientos

que superan los objetivos de este curso. El valor exacto es

z"fzr
f" f{r)o*=
y el error cometido arriba es menor que una millonésima.

5. Aplíquese la fórmula de integración numérica por interpolación en dos

puntos y compruébese que se obtiene la regla del trapecio.

Solución. Enesrccaso, tenemosn= I,a=xs (Jl =bporloque solo

hay que calcular aoY at

En consecuencia

¡b J f\u\+ f(b)
J" ftrldx= La^.f(x^\=aof(xol+a,.f(x,)=(b-a\,

6. Aplíquese la fórmula de integración numérica por interpolación en tres

puntos y compruébese que se obtiene la regla de Simpson.
 Soluciones

Solución. Tenemos xo = a < 'Y, =! < xt = b' es decir' los nodos defi-

nen en el intervalo [a, b] dos subintervalos

de amplitud + . Calculamos

numérica por interpolación

los coeficientes 4s, at Y azpara la integración
en tres Puntos.

l' t,^(x)rl¡ f' t"--1,)(.r--r,)d'r,',1b-a)t b-a

qu-
pn(.ro) (.ro-r, )(xn-x') l","lt't-01 o
eb

I P,(x) dx I
(¡-¡o) Q-x.)dt
,, to
u1
'--
P, ("rr )
=
(Ir ro) (r, - 12)

t.

I P, ( r) r/'t = )"
(, x) (x - x,) dr f (0-ül b-a
,' =@;), =- 6
U1
--'u
-
príxr) (,r" - xo ) (-r, - -x, )

En consecuencta

,t, . f{xo) +4f(x,)+.f('r,)

j"' t't"l d.Y =¿u,.f l'r rl = (b - o'T'

7.AplíqueselaregladeSimpsonparaestimarelvalordelaintegral
el resultado obtenido con el valor exacto.
f ,t" 4), dx y compárese
_

se tiene
Solución. Sea/(x) = x (x -4)2' Por la regla de Simpson
ri dx=(5-trtj-]f-u:()-rt- ftll + flll+ ft5l ,c ,,9t1?jj :T'
52

I, t't.tt 6
este caso
Por otra parte, es fácil calcular el valor exacto en

. [r I 8 l'- {,
ax=lof -1 rt +8-rt |t1 =:
f "t." - +l' rlx: Ji 1-r, -8x'+16¡)

171
g. lntegración y derivación numéricas

po-
Comprobamos que el error cometido es nulo' En un primer momento
que el lec-
demos pensar que se trata de una casualidad. Es recomendable
ejercicio sustituyendo el
to. n]odifique por su cuenta el enunciado de este
que el error
integrando poi cualquier polinomio de grado 3 y compruebe
sigue siendo cero.
interpo-
Detengámonos un momento en este asunto. Sea P2 el polinomio
nodos en los ex-
lador áe segundo orden de una funciónl escogiendo los
tremos y en el punto medio de un intervalo la, bl' De la misma manera
que la iegla del trapecio proporciona resultados exactos cuando
el inte-
jranoo es un polinomio de grado uno, la regla de Simpson tendrá un efror
polinomio
nulo al esrimar la integral 12 f f"¡rl-r, siempre que/(-r) sea un
degrado2.SólohayqueobservarquelaregladeSimpsonprocededein-
te{rar Pry que P2 =/, si/es un polinomio de segundo
grado' Sin embar-
gJ, este ejercicio, el integrando/es un polinomio de grado
3' En prin-
óipio no cabía esperar que el error fuera nulo' pues Pz + f
"n (de hecho Pz
que la expli-
tiene grado no mayor qúe 2). Más adelante comprobaremos
cación de este aparente enigma reside en que la diferencia, f (x) - Pz6),
respecto del
aunque no es iáénticamente cero, es un polinomio impar
centro del intervalo, por lo qu" fi{ f{r) - P, ('r)) dr = 0'

e6

8. Estímese el valor de la integral I e-" dx mediante la regla de Simpson

operaciones
en tres subintervalos (esto es, con 7 nodos), realizando las
con una calculadora científica.

Solución. Consideramos los nodos

xo= -6, xl = 4, xz= -2,f¡ = 0' x+=2' xs= 4' xa= 6'

La suma de los valores del integrando f("r) = ¿-" en los nodos

de índice

par es

i ft-rr-) - e'u +e' +et +e-.u =2(e' +e-'u)'

172
 Soluciones

mientras que cle la suma de los valores

de/en los nodos de índice impar
se obtiene

+ ., r - --16+lr( 16
)-'l (Yzt-t)=( '1-o =l+Ze''u,

En consecuencia

fh -,,,..,!rZto* +t,-'u)+2(l +2e'n)-e'o\=2,1155

Ju' "= 3
integral. En un curso avanzado de
No es sencillo carcular el valor de esta
cálculo se demuestru qu" j* e-' tlx = "'li '
Por lo tanto'

fr- , f* -h! , I .-,n - r,

rr, -.- dx ^ e-"'dx = < l0-'u
{t - J*
e '' dx =2 ), e-'
" l,: ¡o-'n

É
luego J_"e-'dx= , "li = 1.772453851 con un
error menor que l0-8' Con-

al aplicar ra regla de Simpson

cluimos observando que er error cometido
en tres subintervalo, *uy grande' casi una unidad'
"t
g. para estimar el valor de la integral
Desarróllese un programa informático
16
'-r tlx
-
meorante la regla de Simpson en n
subintervalos' haciendo
l*e
variar n desde t hasta 20'

incluido en el libro de
Solución. Modificamos ligeramente el código
teoría.

Function f (x)
f=Exp (-x^2)
End Functron
Function SimP(n)

173
 3. Integ ración y derivación numéricas

a=-6
h- 6
h- (b - a) I Q x n)
5I=U
SP=f(a)
Do
If k > 2 * n Then Exit Do
SI = SI + f (a + k * h)
k=kl-L
Sp = SP + f (a + k * h)
k=k+1
Loop
Simp= ((b-a) I G *n)) * (SP+2 * sr - ((f (a) +f (b) ) l2))
End Functron
Excel copiamos los distin-
De nuevo, en las celdas Al--A20 de la hoja de
I
tos valores de n- En la celda
gl- escribimos =Simp (A1) ' Finalmente'
en Edición-Rellenar-
seleccionamos las celdas Bl--820 y pulsamos
en lás celdas 81-B20 nos aparecen las
Hacia abajo. De esta manera,
aproximaciones de la integral'
,,.r. .. ..: ,, :, :_ :_ .

Sinpstmen¡ suninte*

I 8 t1 1,712149508
2 2,0009872'78 t2 t,7'72392133
a
J 2,71550897r 13 1,777464867
4 t,421843124 l4 t,'77245212
5 1,989171735 15 t,772454088
6 t,6'72428605 r6 t,772453823
7 1,813568907 t7 1,712453854
8 t,757149279 18 t,'7'72453851
9 t,7'71039826 l9 1,172453851
10 t,'711206823 20 t,'772453851

10.AplíqueselaregladeSimpsonenl'10'20'30'"''l00subintervalospara
aproximar el valor de la integral

174
 Soluciones

f"-Lrk.
Jo 4x'+1296

ComPárese con el valor exacto'

desanollado en el ejercicio 9'
modi-
Solución. Aplicamos el algoritmo
ficandoladefinicióndelafunción'estoes'cambiamoslasprimerasTlí-
neas Por las siguientes

Function f (x)
f = x / Sqr ((x ^ 2) + 1-296)

End Functíon
Function SimP (n)

a=0
b=7'7
Así obtenemos los siguientes resultados

I 49,12059174
2 49,07058593
a 49,01351809
-)
+ 49,00315837
5 49,00143',1Q2
10 49,00008347
20 49,00000514
30 49,00000101
40 49,00000032
50 49,00000013
60 49,00000006
70 49,00000003
80 49,00000002
90 49,00000001
r00 49,00000001

175
3. Integración y derivación numéricas

El valor exacto de esta inteeral es fácil de calcular

l"
"'1'l ¡f t"-llx,
r
+t2961'-" =-17n5-"ln% =g5- 36=49
Jo .rl¡;nX "'n - [ '' -- - L=o

-1t=7,1
11. Aplíquese la fórmula de integración numérica por interpolación en cuatro
puntos igualmente espaciados para obtener un regla de integración numé-
rica (regla 3/8).

Solución. Recordemos que al integrar el polinomio interpolador se ob-

tiene
rb¿lb
f .f tx\ rtx =Z o r.f (, r). en donde ou : !!#!!. En nuesrro caso.
Q o(x o)
poniendo xr= a +;; k = 0, 1,2, 3, resulta
"t"
qo(x) = (x * a - h) (x- a - 2h) (x - a - 3h)
q{x) = (x - a) (x - a - 2h) (x - a - 3h)
qz(x) = (x - a) (x - a- h) (x - a - 3h)
qz@) = (x - a) (x - a - h) (x - a - 2h)

J"
qo(x)dx
uo - -------i----l- -
8o\xo )
ta+Jh

= J" tr' -3x'a-6x'h+3xd +12xah+llxh2 - a' -6a'h*llah'-6h'¡dx _

-6hu
a
[t 4 ^1
I x -x' a-2x'h+--- x'a' + 6x'ah ¡l-)
3
- a' x - 6at ltx -llah2 x - 6h3 x
l4 2 "'¡'

-6h'
o
'r4fl.¡
4 J'
:-----------:- n
-6h' 8
_
I q,(x)dx
10
---------.----.-
'
ut
^
-
4,{x,)
-

176
 Soluciones

*to
("t - 3xt a - 5xt h +3xa' + lOxah + 6xh'
qr(a+ h)
l-' < 1
l]IA .r' - "' u-1 r' n + ) x' a' + 5x' ah+3x' h' - a
2lt'
0/ 1J

n' 9
h
.

=A-=-
2h' 8

n
I q,t¡ ldr'
.o
u.-
-- =
q2\.x2 )

- 'l -4ctth-3ah- )dx

fn"n (.' -3.tta- 4-rtlt+3xt¡t + Sxah+3'rfit -
- J,,

Qr\a +2h)

Ir x" -x'a- 4 ,ah+1 xth. -a.x-4a 'nx - unt x)

'h+3 xta'+4x
L¡ ax z
o l-l
-zn
9,4
_-Á'=2t
-2h' 8

lo q'1x¡rtx
.1
u1 -
-'fl -
qr(.r. )

*to
("t - 3xt a-3xt h +3xa' + 6xah + 2xlt' - a' -3a'h-2ah')dx
q.(a+3h)

Ir x' - x'a- xt lt+ 2 x'at +3xtah+x'h' - a' x -3at h, - Zont ,)

ti
3-

L4
6h'
o
--A n .J, I

= -
I

=-11.
-6h' 8 l

177
3. Integración y derivación numéricas

En consecuencia
) + 3.f(-r,) +3.f (x.) +./(¡,)
li rota"=f oo f G):3h19
es decir
pb f'( x ) + 3 f(x,) + 3.f (x.)+./(-t,)
I
tr
ft.r)r/r =1fi-u)L

12. compruébense analíticamente los siguientes resultados, que aparecen en

la deducción de la fórmula del error para la regla de Simpson
_rnX¡ - -r,X'r- xr) clx :o
I' G

l" l{-r - xo xx - -1, )'('t' -x' )l dx =+

en donde a < b y,como es habitual, xo = a' x, =T' x: = b'
solución. Es geométricamente evidente que la primera integral debe
pun-
valer cero, porque la gráfica del integrando es simétrica respecto del
to del plano (x', 0), de forma que el área de la región que queda por deba-
jo del eje,r se compensa con el área de la que está por arriba'
Para realizar la comprobación analíticamente, aplicamos en la primera
in-
(b
tegral el cambio de variable x = (b - a)u * a, de manera qu'e dx = - a)du

¡b lb a+b¡¡ t\ J
],, f,--{o)('r --r')(-r- r)dx = I:,* - a)(x -"i)6 - b)dx =

= f ttA - a )tt + o - ct) ((b - a)u + a - *¡ <fO - a)u + u - b) (b - a) du =

!Í) (b - a)u - (b - a)) du =

= (b -r) l, tt¡ - a)u) ((b - a) u -
= (b - a)'
, fl
u@ - l¡ f" - I) du : (b - a)o., .|ofl r'r -':I u') + \ u) du =
(
|

l'
= 1b - d¡o l!u' -'¡ r' + "' li=i
t, : (b -a)- (
i - i * ])= o

178
 Soluciones

variable para calcular la otra integral

Aplicamos el mismo cambio de

- ('r -'t, )' (-T - -r,)l c/'r = -f (-t - -ro) (x - -r')' (-x - x ) ch =
ij t,-, -ro)

: -il,,, - a) (-r - \r)' ( x - b) dx = (a - fi' I' u@ -'¡)' tu - 1) du =

.¡( , 5 ' - I ), =
= \a - h)' .|, u' -21l' +
I\
-ti 7' )u^
..[t, 5...-1"
l, +-u 1,,,-1"='-
-(t¡-bl'l ,a' --'1 =
1,,=u

,.,,i1-1*l-l)= ú-n1'
|
-(rr_,,' l5-t-:f 8,l l20,'
valor
un solo intervalo para estimar el
13. Aplíquese la regla de Simpson en
O"*t"r"O"fJr'"n-r¿lx'Acoteelerrorcometidomediantelasdosfór-
acotacio-
Simpson' Compárense estas dos
mulas del error para la regla de el valor de la
nes con el error exacto' que podemo'
tonot"i calculando
integral sin aProximaciones'

de Simpson
Solución. Aplicamos la regla

J r.n, a, -
offi =+ =z'os43ssro3

:b-t'e+ en donde a = 0'

del error es € <
La primera fórmula de acotación '
absoluto de la derivada tercera del inte-
b = It ! $&-Wa pota del valor
el intervalo [0, z]; escogemos Mt=
1' de ma-
,.undoü{,.ffii
nera que

e< "t' <0.51

- r92

179
3. Integración y derivación numéricas

La sesunda fórmula de acotación es e < t-,#' , siendo M+ una cota del

*)
valor absoluto de la derivada cuarta del integranOci¿ txl.- , -- cos x, en el
intervalo [0, zr]; de nuevo escogemos Ma = l, de forma que I
I

' jr
e '"
t'
<2880 <0,1
I

por-
Finalmente, determinamos el error con exactitud, 1o cual es posible
que podemos calcular el valor exacto de la integral

f t"n ,t rh: [-cos xl',=!, =2.

El valor exacto del effor es € =){ -2=0,094395103. Por supuesto, las

fórmulas de acotación del error son útiles cuando ignoramos el valor
exacto de la integral.

Determínese, mediante la fórmula de acotación del error en la regla

del
14.
procedi-
trapecio, el número de intervalos necesario para estimar por ese
miento el valor de la integral

¡'fi*r a'

con un error de truncamiento menor que una centeslma'

solución. Al estimar el valor de la integral mediante la regla del trape-

cio en m subintervalos, el error absoluto está acotado como sigue

_ - (8 -o)'M,
' =-- n*'
en donde l,f'(") I 3 M2paratodo -r e [0, 8]' siendo .f lx) = Vr.+
t Vamos

a hallar el menor valor de m fal que

180
 Soluciones

8'M. 1

lzr' - roo
se obtiene
Despejando m enladesigualdad'

nt>10
ry.
3

ParacalcularlacotaM2,derivamoselintegrando

.f (.r)=V.+I=1;¡*l)"'
.f''(r) =f tx + t)-"'

Í'"(x)=]r"*l)-"'= ;,;;1
de manera que
Podemos escoger M t =2o'

160
,rr >
^^,{
Bo!Í = - 3o:,szo143i
=
JF
intervalo l0' 8lcn 3l subinterralos'
por lo que bastadividir el la integral ( l3'29056581 )'
el valor i*""cle
En este ."'o, poo"*ot "uttutu' 13'28914568' lue-
l" 31 subinterualos sale
Al aplicar "n que una- centésima' De hecho'
'"gl;;i';;;tiobastante
go el error cometido es -"nn' para conse-
*t" con dividir en 12 subintervalos
puede .otptnuu"Jq* supuesto' esto últinro
que una centésima' Por
guir que el error sea menor
el valor exacto de la integral'
no lo sabría*o'
'ü"'*tiltamo'
error en la regla de
rnediante la fórmula de acotación del
15. Determínese'
Simpson,"tnnmerodeintervalotn"t"J*ioparaestimarporeseprocedi-
mi"nto el valor tle la integral

181
3. Integración y derivación numéricas

l"-Lrtx
Jo r/rr + 1296

con un error de truncamiento menor que una centésima'

Simp-
Solución. Al estimar el valor de la integral mediante la regla de
son en m subintervalos, el error absoluto está acotado
como sigue

t =ú-u)'M, -- -\77-O)'M1
zs8o;: 2880"f

cuarta del inte-

en donde M¿ es una cota del valor absoluto de la derivada
. Vamos a hallar el menor valor de m tal
que
grando (x).f ="..#*

775M1-<0,01
2880m'

Despejando m enla desigualdad, se obtiene

la derivadas:
El cálculo de la cota M4es laborioso. comenzamos calculando

\.-, -,-Jr,- '

fl-r):
.,
+1296
,: -t(x' +I296)-t''

.f '(x)= (x'1 +1296)-"' - l-t {x' 11296\-^''2' =

L

= (x2 +1296)-''2 ((-r' + 1296) - x') = 1296 (x' +7296)-'''

f " (x) :-3888-r (x' + 1296)-5''

.f"'(x) = ¡15552x')- 5038848) (x' +1296)-'''
.f (¡) = -(777 60 x' -
n,
7 5582720 -r) (x' + 1296)u''

182
 Soluciones

planteamos un problema
de máximos y
Para acotar la derivada cuarta'
*ini-or, hallar los extremos de la función
'75582720 x) (x2 + 1296)-et2
8(.r) = - (77',760 x' -
extremos
éstos.se alcanzarán en los
[0,77]' Sabemos que
en el intervalo
delintervalooenalgunodelospuntoscríticosdelafunción.Porlotanto'
necesitamos derivar de nuevo
+ 97955205120) 6 + 1296)-"t2
g@) = ¡st @) = (466560¡o - 90699 2640;
de igualar
resolviendo la ecuación resultante
Hallamos los puntos críticos
a cero la función g'(x)

gx6gg2640x)+97955205120)(x2+1296)-'1/2=0+
@66560x'-
+ 91955205120 = 0 +
= 466560x' - 90699z640¡'
2\466560)
30s 80ü f = 9'72 ú2a"'11
so6sg264ot 302 -r
i ^-^ ,o¡^[7
= ,'=972-324"11
=972!32447)
933 120

dentro del intervalo l0'771

Así obtenemos dos puntos críticos

,, = l8rFlT, "=n^FtT
t'(x)ldebe valores:
ser el mayor de los siguientes
El valor máximo de l,f

l f '' (o)l= o
(ll)l= 1,281431998 x 10-i
l.f " I

r' ('*ú1',T )l= s.utat5488x ro '

I

*f,: "''z = +' zooro8 r e3 x r o*

r' ( )\
I '

183
3. Integración y derivación numéricas

Por lo tanto, podemos escoger M+= 5 x 10-'. Finalmente

* r rW-=
\l 2880 x 0.01 rr r@
\| 28.80
= 4.655944402

por 1o que concluimos que basta dividir el intervalo 10,l7l en m = 5 sub-

intervalos.
En el ejercicio 10 comparamos el valor obtenido al aplicar la regla de
Simpson con el valor exacto. Obsérvese que con m = 5, efectivamente, el
enor es menor que una centésima. Podemos observar que esa cota se
consigue mejorar eon m = 4, pero la utilidad de cálculo desarrollado en
este ejercicio estriba en la posibilidad de efectuarlo cuando el valor exac-
to de la integral es desconocido.

16. Para estimar el valor de la integral l* ,-" clx, seaplica la regla de Simpson
en el intervalo [-3, 3], dividido .nlrn.o subintervalos iguales. Determí-
nese una cota del error cometido.

Solución. Por una parte,

" , fr " , ^ fF ,r'

f*- -,,
J-e dx= J_re dx+¿J\ e' dt
luego el error €1 cometido al integrar en el intervalo [-3, 3] en lugar de
hacerlo en toda la recta se puede acotar de la siguiente manera

,, =2 I^ u-" r1, <21,* r-'. a* =!r-' : 8,22732027 4 xl0-5

Por otra parte, el error €2 que se comete al aproximar la integral en el inter-

valo [-3, 3] mediante la regla de Simpson en cinco subintervalos verifica

(3-(4))5 M,
-1
€2 =27 *"
2880x5u 6250

184
 Soluciones

cuarta de
valor absoluto de la derivada
en donde M¿, es una cota del la
3l' Al derivar cuatro^veces
la función .f(x) = ¿" "l "l intervalo t-3' se obtiene f
o' (') :
u lu
función, sacando t'actor común "*pon"ncial' en el intervalo [-3' 3]
de esta función
= 4 (xo - l\f + Z+) u:i.gtmáximo
un punto crítico' Pa-
de los extremos' o bien en
se alcanza, bien en alguno a cero y re-
derivamos de nuevo' igualamos
ra hallar los puntos
"'J'i*''
solvemos la ecuación resultante

(¡) = -8.x (4xo -20xt + 15) e-.' = 0

f ',

la ecuación
cuatro son las soluciones de
Una solución es x = 0, las otras
ui"r"¿tu¿" 4x4 -2ox2 + 15 = 0' operando
2A! 5 r-Ji0
2x4

de la derivada cuarta son

*:|^@+ffi
De manera que los puntos críticos '

x=:úortffi *=-l.1llllffi' x=0' En

x=-l.FoJZF, ' 3] es el mayor
(x) en el intervalo [-3'
consecuencia, el valor máximo/o)
de los siguientes valores

='''**'oo*
\,.,[1.rr..m)l
=
l
r-
[ ; +".'m) \

o"or""
(; l'-'[ ;ño-@)\=''
=
1'' a--um-¡) I

t
l./-'to) | =
l./-'t:l |
= l f" f-:l |
= 0,108106e884

Mq=8' de manera que

Podemos entonces escoger
M'4=0.03456
o,.27
-'- 625A

185
I
3. lntegración y derivación numéricas

Finalmente, el error total verifica

e(er*e2<0,035
Podemos hacer una comprobación del resultado obtenido recordando que
el valor exacto de la intesral es

Jl e-" dx = ',fi = 1,772453857

mientras que, de la aplicación de la regla de Simpson en 5 subintervalos

sobre el intervalo [-3, 3] se obtiene

r+3

J_,
e-' dx =1,773647733

luego el enoÍ exacro es 0,001 193882.

17. Las reglas del trapecio

y de Simpson se han deducido integrando un poli-
nomio interpolador de grado uno y dos respectivamente. Deducir una
nueva regla de integración numérica, para aproximar l" .f {xl d-r:, interpo-

lando entre los extremos del intervalo una función exponencial de la for-
ma q(x) = Aeu'.

Solución. Se trata de integrar la función q(x) = Aeu' en vez de /(r), es-

f
cogiendo las constantes A y B de manera que q y coincidan en los extre-
mos del intervalo.

FD pb ¡b ApBb - AeBo
Ae'' dx=lleu^1)='":L
I' t<xlax= l¡otxldx= J" u-
La función q debe verificar q(a) = f (a), q(b) = f (b), es decir, debe coinci-
dir con/en los extremos del intervalo Ia, bl. Luego

Ora. = f (a)
Aeto =f (b)

186
 Soluciones

por la primera' resulta

Dividiendo la segunda ecuación
.f (b\
'^Srt,-", -.¡ta\
-

tomando logaritmos

B(b- u\= Int.f tb)) - 1¡1 f1a))

'"[#)=
de donde resulta
ln(f(b))-lnt.ft¿¡))
D = --------
b-a

luego
eb
rl d, =
Aer''' - or','= -ñf Ú\- f tu\
*;1'";
l" f t

-T-
regla de Lyday
Finalmente, obtenemos la llamada
.b tb-u)(f(D- f\u\)
\-'l-".' In(.f(b))_lnt.fta))
IJa.'f(.\r¿¿.t=---;-

Gauss con tres nodos'

18. Deducir la fórmula de cuadratura de

Solución.Razonamostalycomoseindicaenlasección6detlibrode
teoría.
a) h(x) = x3 * bzxz* bé * bs'
b) Planteamos el sistema
J',
t"' * h,.rt + h,'r + bn)dr = o
r{x' + 1t..1' + h,x + h,)¿/x = 0
J',
+ ho)tlx = o
f, "t {"t + h"-r'
+ b,x
!

187
3. lntegración y derivación numéricas

que operando resulta

g
!b, +2bo =
i+!u,=o
?u,*i bo = o

Lasoluciónes bo = 0,br=-1,br=0.
c) Las raíces del polinomio Q("r) = x' - -1-r son

IT
- - t-
xt=0, ^'-15

ü Los pesos de la fórmula de cuadratura de Gauss con tres nodos son

)-
a, _L-:-
oo: qo(rl dx
I'-, Jl, (" - x, (x - xr) dx _ J1,(¡ -o)(¡ -.1[t,>
q

Jr, rJx, -,; - (-f -orr-^E -^f t 9g

)
^(x) = -
)

I'-,q,(x) dx Il,("-ro)(x- Il,('*rE¡1x-f¡dr -8rs au

xr)dx
ot=-;Jrr\ = (x; xr)G, - (o+1f,X0-\g) :::_1 -q
ü 5-

q,(x) dx
J:, Jl,(x-.ro)(x--r,)d,r Ji,(x+\E )(, -o) ¿' _ :\
2-
3 _"
-'¿
a^=%:¿
qrkr) - -
(xr- xo)(xr-xr) r.€ *rfiltf -ol
- ó -O
a'

Finalmente

f- r <u>¿u = l rr-f I *
| rrol *
f rt
^f,
I

188
 Soluciones

9. Estímese el valor de la integr al fre-r' dx mediante

la fórmula de cua-

dratura de Gauss con 2 nodos.

Solución' sea g el integrando 8@) = ¿-"' Mediante el cambio de varia-

ble x ='a?"*'-?=?¡u+| se tiene

(t r)'
-t -ü+ |

: e'"'du
f-,s G)
ax
:l'-,
g 1)u + |) du =
; I'-,

Si
,i:,.+)'
.f(r):¿ l"l , entonces
:

I l"r+ll
ll ^. .,
(.1l+ ¡1.'.)=¿-i'n- "
t tl
+¿ r'*
"' = l'029349143
J-,.t\unu=.f
Finalmente
?¡l
J_,e-"
dx =;J ,.f
tu\du = 1.544024614

20. Estímese el valor de la integr al fre-" dx mediante la fórmula de cua-

dratura de Gauss con 3 nodos.

solucién. Definamos/y g como en el ejercicio 19. Entonces

I'-,rtu¡du=|rr_q,.ir rot*J rtf I =

= j *ir-i +;r-i-*ro = 1,085839-57
"-i.i'*
Finalmente

(u) du = 1,6287 59355

f,r- " dx = : I'-,f

189
 Integración derivacíón
3. Y
'Y-éñ"u
que
Se puede demostrar que el error de esta aproxima:tó".tt menor
e- '' dx = l'628 con cua-
2 x l|u,por lo que se puede dar el resultado J-,

tro cifras signifi cativas'

2l.Desarrolleunprogramainformáticoparacalcularnuméricamentelaintegral

l"-Lrt"
Jo 4x'+1296
de Gauss con dos nodos' dividiendo
mediante la fórmula de cuadratura
previamente el interval'o en 1' 10'
20'30' "'' 100 subintervalos' Compá-
en ejercicio l0'
rense los resultados con los obtenidos
A1t-F11 y es-
Solución. En una página en blanco de Excel pulsamos
que se abrirá
cribimos las siguientet"in't*"tiones en la ventana
Function I $) As Double
g=x I Sqr((x ^ 2) + 1296)
End Functron
Function gaud(n) As Double
gaud = 0

A- a

For k= 1 To n
^-á
d=c+ ((b-a) ln) (-0'5'713502693) )+,,,s ,.\\
((d + cl /¿tt
f L=s(( ((d - c) l2) * + ((d + c) I 2\)
f 2= 9(( ((d- c) I 2) * (0 '5'713502693))
(f 1 + f 2)
gaud.=qaud + ((d - c) I 2) *
Next
End Function
R de la hoja de cálculo los va-
A continuaciÓn, escribimos en la columna
(1, 2, 10, 20' "'' 100) y en la casilla 81
"',
lores del número oe intervalos
Finalmente' seleccionamos en la
columna B
ponemos = gaud(A1 ) '
mnmscasillur.o.no-huyamosrellenadoenlacolumnaAypulsamos

190
 Soluciones

Así obtenemos los siguientes re-

sdición-Rellenar-Hacia abaj o '
sultados, que se comparan con los del
ejercicio 10'

Simpson,(2n + L nodos

48,9'1548302 49,12059174
I
48,95223413 49,07058593
2
48,99082570 49,01351809
J
48,99',7467',78 49,00375837
4
48,99903605 49,00143702
5
48,99994428 49,00008347
10
48,99999657 49,00000514
20
48,99999933 49,00000101
30
48,99999979 49,00000032
40
48,99999991 49,00000013
50
48,99999996 49,00000006
60
48,99999998 49,00000003
70
48,99999999 49,00000002
80
48,99999999 49,00000001
90
48,99999999 49,00000001
100

SeobservaenestecasoquelafórmuladecuadraturadeGausscondos
nodosencadasubinterval.oproporcionaunosresultadosligeramenrcme-
jores que la regla de SimPson'

22.Desarrolleunprogramainformáticoparacalcularnuméricamentelaintegral
.¡ r
J,;,r"
mediantelafórmuladeouadraturadeGausscontresnodosencadasub-
Compárese cada re-
intervalo, considerand o 2, 4,6, 8 y 10 subintervalos.
sultadoconelqueseobtienealconsiderarelmismonúmerototaldeno-
de Gauss con dos nodos en
dos, pero aplicando la fórmula de cuadratura

I
191
3. lntegración y derivación numéricas

cada subintervalo. Utilícese para ello el programa desarrollado en el ejer-

cicio 2t.

Solución. Partiendo del programa desarrollado en el ejercicio 2I, a¡la-

dimos las líneas necesaria para que aplique también la fórmula de cua-
dratura de Gauss con tres nodos en cada subintervalo' Quedará
Function g(x) As Double
g-(x^5) lQ+x^3)
End Function
Function gaud(n) As Double
gaud, = 0
--1
¡.- ?

d=a
For k=1 To n
^- A
d=c* ((b - a) / n)
f I= s((((d- c) / 2) x (-0.s773502693) ) + ((d+ c) / 2))
f 2=s( (((d-c) /2) * (0.s173502693))+ ((d+c) /2\)
gaud= gaud* ((d - c) / 2) * (f L+ f 2)
Next
End Function
Function gaut (n) As Double
gaut = 0
1
^-
d- L
1r- ?

d=a
For k= 1 To n
^- A
d=c* ((b-a) /n)
h1-= s(( ((d - c) I 2) * (-0 .7745966692)) + ((d+ c) / 2))
h2=9((d+c) /2)
h3 = s(( ((d - c) / 2) * (0 .1145966692)) + ((d+ c) / 2))
9aut.=gaut+ ((d-c) /2) * ((5 /9) *h1* (8/9) *h2+
(5/9) *h3)
Next
End Function

192
 Soluciones

Ahora en las celdas A1-A5 de la hoja de cálculo Excel escibimos 6,

12,
nodos empleados en
18,24y 30, es decir, en cada caso el número total de
la fórmula de cuadratura de Gauss con tres nodos en cada intervalo,
divi-
diendo el intervalo [1, 3] en 2,4,6,8 y 10 subintervalos respectivamente'
En la celda Bl- introducimos = gaud(A1/2) , mientras que en la celda
C1 ponemos - gaut(A1l3). Finalmente, seleccionamos desde Bl-
j-a abaj
hasta C5 y pulsamos Edición-Rellenar-Hac o'

Gauss con 2 nodos Gauss con 3 nodos

en cada subintervalo en cada subintervalo
n/2 subintervalos z/3 subintervalos

6 7,18684630 1,78699',731
l2 1,78697461 7,18698153
18 7,78697974 1,78698095
.A
L.+ 7,78698054 7,78698090
30 7,',78698075 7.78698089

Es interesante compafar estos resultados con el valor exacto de la inte-

gral. Efectuando la división e integrando después se obtiene

I r' , ,( , x' ,.. -1r'-ltn,l+x,)

'1

I,u"=Jl" r*r..)o'=t^ 3¡r'\¡'¡'r '

lueso
¡3f "
5
)Á1
,d.r= '-! - -tn2 - :-ln7 = 7, 78698089
1

|Jt l+ir' -1 -1
I

con nueve cifras decimales exactas'

23. Aplíquese la fórmula de las diferencias divididas centradas con tamaño

depasoh=0,0Iparaestimarelvalordeladerivadaf,(0,25),endonde/
es la función f (8 = {x. Compárese con el valor exacto'

193
 g. lntegración y derivacíón numéricas

Solución. Calculamos con cuatro cifras decimales

0,5099 - 0,4899
r''rt'l )lt = f t0,261- f \0,24) _Jofr-"{07 = 1,0000
0'02 0,02 0,02

de/es .f '(x) = luego el valor exacto de

Por otra parte, la derivada
*
/'(0,25) es

.f'(0,2s)=#=t
Hemosobtenidounerrornuloporqueloserroresdetruncamientoyde
numerador'
redondeo se han cancelado en la resta del

divididas centradas pro-

i
24. Compruébese que la fórmula de las diferencias
porcionaresultadosexactosparaladerivadacuandoseaplicaapolino-
mios de grado menor que tres'

Consideremos la función,f (x) = of + bx + c' en donde a' b y

Solución.
Csonconstantes,querepresentaunpolinomioarbitrariodegradomenor
que tres.
Sabemos que su derivada es/'
(x) = 2ax + b' Apliquemos ahora la fórmu-
la de las diferencias divididas centradas'

f (x+h)-.f (x-h) +b(x+h)-a(x-h)2 -b(x-h) =2ax+b

r.-=
-a(x+h)2 2h
2h

Estímese el valor de la derivada en el

punto x = 0'9 de la funciÓn
25.

(l +x) |
(,+-") l,t+;,)
ftx) = lím
n4e
I
./
1...1 l'
\ - ,/
\

utilizandolafórmuladelasdiferenciascentradasdivididasconincre- pa-
cómo evaluar la función.f
mento h = 0,01. Estúdiese con detenimiento

194
 Soluciones

deri-
posible calcular el valor de esta
ra obtener la estimación pedida' ¿Es
vada de forma exacta''!

Comencemos observando que

la función / no está generada
Solución. básicas' pues
que caracterizan a las funciones
de acuerdo con las norrnas de polinomios'
aparece un producto infinito
en la expresión que la define de derivadas para
direclamente la tabla
Por lo tanto, no podemos aplicar tratada
hallar/'(x)' Sin "tuu'gá' ia función/puede ser razonablemente
mediante métodos numéricos' 1 r:-.:r tA¡c ¡¡
Aplicandolasfórmuladelasdiferenciascentradasdivididasconrncre-
mentoh=0'01 tenemos

f (0, 9 * 0, 0 1): /(0f:991) = 5( f(0.91) - f(0,89))

.f '(0,9) =
2h

habrá que computal el pro-

para estimar los valores de/(0,91) y/(0,89),
m de factores en la definión de /'
esto es
ducto de un número *nti"ni"
f
"\ I x)
f(x)=0+r)[r.;,J It*T).
decidir qué número m es suficiente'
Por supuesto, es un asunto delicado
gu.tu;u grande como para que
escoger m lo suficientemente
I

'.-f)l'.=-- I l'.*l<
r'oor
f
\''r'J[^'z'.'J'\ 2") :

se puede conseguir'
para todo , , *,y x e [0'8 ' I ]' por ejemplo' Esto
es la serie de Taylor de
porque I + / < e' paratoio 01recuérdese cómo
't
función exponencial), lo que implica
que
ia
( x )f.. t )1,*¿)<e''i-' :e"''-2^t.rt)'^t
['.*.i['.*,J [.'*tJ-'
Bastará escoger m = l2,pues e't"' < l'001'

195
3. lntegración y derivación numéricas

Sobre una página en blanco Microsoft Excel prisamos

Alt-Fll' para
Insertar de este editor se-
acceder at e¿ito. de visual Basic; en el menú
en la ventana en
leccionamos Módulo y escribimos el siguiente código
blanco que aParecerá:
Function f (x, n)
Dim k As Integer
F-1
Fork=0Ton
f = f * (1 * (x / (2 ^ k) ))
Next k
End Functron
Functionder (a, h, n)
der=(f(a+h,n) - t (a - h, n) ) / (2 * h)
End Functron
A continuación, escribimos en las celdas A1 ' A2 y A3
respectivamente
elpunto(número)aenelquequeremosestimarladerivada,elincremen-
fo^h,yelnúmerodefactoresnempleadosparaevaluarlafunciónlFinal.
mentl, en la celda 44 escribimos = der (A1 ; A2 ; A3 )
y obtendremos en
esaceldalaestimacióndef'(a),calculadaconincrementolzyevaluando
f con n factores. en la A2 ponemos 0'01 y
En nuestro caso, ponemos 0,9 en la celda ,\1 '
5,48. Es interesante
en la A3 escribimos 12. Así obtenemos que.f'(l) =
comprobarqueestaaproximaciónpermanecebastanteestabletantosiau-
mentamos t?, como sr disminuimos h'
Nota.Aunquenoestrivial,laderivaclade/puedecalcularseconexacti-
md de la siguiente manera

-a-
.f(.r)=f[tt+2 ".rt lnt.f(;)l=) tntl+2"x) =)
il tu ' ¡=ir

=i
f'(r) 2-" :i i
- ft.rl ftt+2"t i, r-,f\ñ r-r-,.r)* )l=
:á,-,,-[i {2'-')'')=¿ #
196
 Soluciones

de donde

(-{,g)* )
.f'(o,g)- t,o,nrfá T;; )

Ya hemos visto que para aproximar/(0,9) es suficiente tomar 12 factores

en la fórmula que define a la funciónf sin embargo, para alcanzar un re-
sultado tan preciso como el 5.48 calculado antes, necesitamos computar
más de 60 sumandos en la serie anterior. En efecto, si tomamos n suman-
dos, la parte despreciada de la serie será

iL (-o.e)*
i =
l_Z_i_, o?.r, ".-,
,_o.s),=
(1.?):'
l+0.9

Para conseguir que esa cantidad sea menor que una milésima, ¡z debe va-
ler más de 60

It-o.g'"-'l ,, togto.oolg)
o'ool 4+ tt+ =59.5.
lñ-l< roe(o.e)

En conclusión, incluso cuando es factible obtener la expresión exacta de

la derivada, puede resultar conveniente aproximarla mediante intelpola-
ción (esto es lo que realmente hacemos al aplicar la fórmula de las dife-
rencias divididas centradas).

26. Estímese el valor de la derivada de la función/(-r) = x e" cos f en el pun-

fo x = 78, interpolando con tres nodos y tamaño de paso h = 0,0I'

solución. Aplicamos la fórmula correspondiente para obtener

Í'(n)
. =-3f \n)+4f (n+0,0\)-.f(n+0,02) :
4,62447 568
_218,09589 -294,6372515 +7 = -95.844291
a,02

197
g. Integración y derivación numéricas

exacto:
Comparemos el resultado obtenido con el valor

.f '(x) =¿.' cos -x + f¿' cos x - xe' sin x

Por lo tanto, -f '(tt) = -95'83932262'

función,f (x) = 'r ¿' cos f en

27. Estímese el valor de la derivada segunda de la
de paso h = 0'01'
el punto x = ft, interpolando con tres nodos y tamaño

Solución. Aplicamos la fórmula de la derivada segunda

- 0,0t¡ - (n)+.f(zr + 0'01)

f"(n) =
ftzr 2.f
-
1,7 4257 488 + | 45,.3-9-7 26 - 7 3' 6593 1288
- -7 = 46,27 7 6
0,0001

exacto:
Comparemos el resultado obtenido con el valor

.f " (x) = 2e' cos x -2e* sin ¡ - 2-r¿' sin ¡

Por lo fanto' 'f "(7t) = 46'2813852'

198