1) Trace o gráfico e isole as três primeiras raízes positivas da função: f(x)=5 sen(x²)-exp(x/10)

Para f(x)=0, temos 5 sen(x²)=exp(x/10)

Fazemos g(x)= 5 sen(x²) e h(x)=exp(x/10)

x 0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1 1,1 1,2 1,3 1,4
g(x) 0 0,05 0,20 0,45 0,80 1,24 1,76 2,35 2,99 3,62 4,21 4,68 4,96 4,96 4,63
h(x) 1 1,01 1,02 1,03 1,04 1,05 1,06 1,07 1,08 1,09 1,11 1,12 1,13 1,14 1,15
Sinal f(x) - - - - - + + + + + + + + + +

x 1,5 1,6 1,7 1,8 1,9 2 2,1 2,2 2,3 2,4 2,5 2,6
g(x) 3,89 2,75 1,24 -0,49 -2,26 -3,78 -4,77 -4,96 -4,19 -2,50 -0,17 2,29
h(x) 1,16 1,17 1,19 1,20 1,21 1,22 1,23 1,25 1,26 1,27 1,28 1,30
Sinal f(x) + + + - - - - - - - - +

Da tabela acima, obtemos as raízes: x1∈[0.4, 0.5], x2∈[1.7, 1.8], x3∈[2.5, 2.6].

Abaixo usamos o método da bissecção para aproximar x1, com erro<10-5.

Iteração a b c f(a) f(b) f(c) b-a

1 0,4 0,5 0,45 -0,24421974 0,185749 -0,04043 0,1
2 0,45 0,5 0,475 -0,04043349 0,185749 0,069932 0,05
3 0,45 0,475 0,4625 -0,04043349 0,069932 0,014057 0,025
4 0,45 0,4625 0,45625 -0,04043349 0,014057 -0,01336 0,0125
5 0,45625 0,4625 0,459375 -0,0133621 0,014057 0,000304 0,00625
6 0,45625 0,459375 0,457813 -0,0133621 0,000304 -0,00654 0,003125
7 0,4578125 0,459375 0,458594 -0,00653976 0,000304 -0,00312 0,001562
8 0,45859375 0,459375 0,458984 -0,00312045 0,000304 -0,00141 0,000781
9 0,45898438 0,459375 0,45918 -0,00140875 0,000304 -0,00055 0,000391
10 0,45917969 0,459375 0,459277 -0,0005524 0,000304 -0,00012 0,000195
11 0,45927734 0,459375 0,459326 -0,00012409 0,000304 9,01E-05 9,77E-05
12 0,45927734 0,45932617 0,459302 -0,00012409 9,01E-05 -1,7E-05 4,88E-05
13 0,45930176 0,45932617 0,459314 -1,7002E-05 9,01E-05 3,65E-05 2,44E-05
14 0,45930176 0,45931396 0,459308 -1,7002E-05 3,65E-05 9,77E-06 1,22E-05
15 0,45930176 0,45930786 0,459305 -1,7002E-05 9,77E-06 -3,6E-06 6,1E-06
 Na tabela abaixo usamos o método da bissecção para aproximar x2, com erro<10-5.

Iteração a b c f(a) f(b) f(c) b-a

1 1,7 1,8 1,75 0,059429 -1,68846 -0,7962 0,1
2 1,7 1,75 1,725 0,059429 -0,7962 -0,36224 0,05
3 1,7 1,725 1,7125 0,059429 -0,36224 -0,14969 0,025
4 1,7 1,7125 1,70625 0,059429 -0,14969 -0,04468 0,0125
5 1,7 1,70625 1,703125 0,059429 -0,04468 0,007489 0,00625
6 1,703125 1,70625 1,704688 0,007489 -0,04468 -0,01857 0,003125
7 1,703125 1,7046875 1,703906 0,007489 -0,01857 -0,00553 0,001562
8 1,703125 1,70390625 1,703516 0,007489 -0,00553 0,000981 0,000781
9 1,70351563 1,70390625 1,703711 0,000981 -0,00553 -0,00228 0,000391
10 1,70351563 1,70371094 1,703613 0,000981 -0,00228 -0,00065 0,000195
11 1,70351563 1,70361328 1,703564 0,000981 -0,00065 0,000167 9,77E-05
12 1,70356445 1,70361328 1,703589 0,000167 -0,00065 -0,00024 4,88E-05
13 1,70356445 1,70358887 1,703577 0,000167 -0,00024 -3,7E-05 2,44E-05
14 1,70356445 1,70357666 1,703571 0,000167 -3,7E-05 6,51E-05 1,22E-05
15 1,70357056 1,70357666 1,703574 6,51E-05 -3,7E-05 1,42E-05 6,1E-06
Na tabela abaixo usamos o método da bissecção para aproximar x3, com erro<10-5.

Iteração a b c f(a) f(b) f(c) b-a

1 2,5 2,6 2,55 -1,4499215 0,997827 -0,20266 0,1
2 2,55 2,6 2,575 -0,20265771 0,997827 0,408766 0,05
3 2,55 2,575 2,5625 -0,20265771 0,408766 0,105173 0,025
4 2,55 2,5625 2,55625 -0,20265771 0,105173 -0,0483 0,0125
5 2,55625 2,5625 2,559375 -0,04829677 0,105173 0,02856 0,00625
6 2,55625 2,559375 2,557813 -0,04829677 0,02856 -0,00984 0,003125
7 2,5578125 2,559375 2,558594 -0,00983926 0,02856 0,009368 0,001563
8 2,5578125 2,55859375 2,558203 -0,00983926 0,009368 -0,00023 0,000781
9 2,55820313 2,55859375 2,558398 -0,0002339 0,009368 0,004567 0,000391
10 2,55820313 2,55839844 2,558301 -0,0002339 0,004567 0,002167 0,000195
11 2,55820313 2,55830078 2,558252 -0,0002339 0,002167 0,000967 9,77E-05
12 2,55820313 2,55825195 2,558228 -0,0002339 0,000967 0,000366 4,88E-05
13 2,55820313 2,55822754 2,558215 -0,0002339 0,000366 6,62E-05 2,44E-05
14 2,55820313 2,55821533 2,558209 -0,0002339 6,62E-05 -8,4E-05 1,22E-05
15 2,55820923 2,55821533 2,558212 -8,3842E-05 6,62E-05 -8,8E-06 6,1E-06

Obtivemos: x1=0.459305, x2=1.703574, x3=2.558212

2) Verifique (analiticamente) que a única solução real da equação: x ex = 10 é ponto fixo da seguinte
função: 𝑔(𝑥) = ln(10/𝑥)

x ex=10 pode ficar: ex=10/x

Desenhando essas duas funções podemos observar que possui somente uma solução real, que é
positiva.
 Comparemos agora com a função g(x)=ln(10/x), que tem ponto fixo em g(x)=x, portanto

x=ln(10/x).

Se aplicamos exponencial dos dois lados da igualdade, obtemos ex=10/x.

Isso mostra que a raiz real da equação inicial é ponto fixo da função g(x)

3) Encontre as raízes do polinômio 𝑓𝑓(𝑥𝑥) = 𝑥4 − 4𝑥2 + 4 através do método de Newton. O que

você observa em relação ao erro obtido? Compare com a situação do Problema 3.2.4

Calculando a derivada: f’(x)=4x³-8x.

f(x ) xn 4 −4xn ²+4
Método de Newton -> xn+1 = xn − f´(xn ) = xn − 4xn ³−8xn
n

Adotando-se x0=1, têm-se:

n xn xn+1 f(xn+1) xn+1- xn

0 1 1,25 1 0,25
1 1,25 1,3375 0,191406 0,0875
2 1,3375 1,376957 0,044561 0,039457
3 1,376957 1,395837 0,010814 0,018880
4 1,395837 1,405086 0,002667 0,009249
5 1,405086 1,409665 0,000662 0,004579
6 1,409665 1,411943 0,000165 0,002278
7 1,411943 1,413079 0,000041 0,001136
8 1,413079 1,413647 0,000010 0,000567
9 1,413647 1,413930 0,000003 0,000284
10 1,413930 1,414072 0,000001 0,000142
11 1,414072 1,414143 0,000000 0,000071
12 1,414143 1,414178 0,000000 0,000035
13 1,414178 1,414196 0,000000 0,000018
14 1,414196 1,414205 0,000000 0,000009
15 1,414205 1,414209 0,000000 0,000004
 16 1,414209 1,414211 0,000000 0,000002
17 1,414211 1,414212 0,000000 0,000001
18 1,414212 1,414213 0,000000 0,000001

Para encontrar a outra raíz, adotamos x0=-1.

n x xn+1 f(xn+1) xn+1- xn

0 -1 -1,25 1 -0,25
1 -1,25 -1,3375 0,191406 -0,0875
2 -1,3375 -1,376957 0,044561 -0,039457
3 -1,376957 -1,395837 0,010814 -0,018880
4 -1,395837 -1,405086 0,002667 -0,009249
5 -1,405086 -1,409665 0,000662 -0,004579
6 -1,409665 -1,411943 0,000165 -0,002278
7 -1,411943 -1,413079 0,000041 -0,001136
8 -1,413079 -1,413647 0,000010 -0,000567
9 -1,413647 -1,413930 0,000003 -0,000284
10 -1,413930 -1,414072 0,000001 -0,000142
11 -1,414072 -1,414143 0,000000 -0,000071
12 -1,414143 -1,414178 0,000000 -0,000035
13 -1,414178 -1,414196 0,000000 -0,000018
14 -1,414196 -1,414205 0,000000 -0,000009
15 -1,414205 -1,414209 0,000000 -0,000004
16 -1,414209 -1,414211 0,000000 -0,000002
17 -1,414211 -1,414212 0,000000 -0,000001
18 -1,414212 -1,414213 0,000000 -0,000001

As raízes obtidas são x1=1,414213 e x2= x2=-1,414213. Em comparações com as raízes √2 e -√2,
temos erro inferior a 10-6.
Em comparação ao método usado no problema 3.2.4, o método de Newton é adequado para
encontrar as raízes pois não depende de valores com sinais diferentes para aproximar a solução e não tem
problema para encontrar raízes duplas.

3.2.4 O polinômio f(x)=x4-4x²+4 possui raízes duplas em √2 e −√2, O método da bisseção pode ser
aplicado a f? Explique.

f(x)=x4-4x²+4.
 Nota-se que f(x)>0 para x∈ℝ, então o método da bissecção não pode ser aplicado diretamente a f(x),
pois o método consiste na busca de uma raíz entre um ponto em que f(x)>0 e outro em que f(x)<0.

4) Aplique o método das secantes para resolver a equação 𝑒−𝑥2 = 2x.

f(x)=𝑒−𝑥2-2x

Método das secantes

x𝑛 − x𝑛−1
x𝑛+1 = x𝑛 − f(x𝑛 )
f(x𝑛 ) − f(x𝑛−1 )

Usando como pontos iniciais x0=0 e x1=1, temos:

n xn-1 xn f(xn-1) f(xn) xn+1 f(xn+1)

1 0 1 1 -1,63212 0,379921807 0,105749
2 1 0,37992181 -1,63212056 0,105749 0,41765353 0,004625
3 0,37992181 0,41765353 0,10574928 0,004625 0,419379155 -3,9E-05
4 0,41765353 0,41937916 0,00462483 -3,9E-05 0,419364819 1,34E-08
5 0,41937916 0,41936482 -3,8744E-05 1,34E-08 0,419364824 3,85E-14
6 0,41936482 0,41936482 1,3365E-08 3,85E-14 0,419364824 0

A raiz obtida é x1=0,419364824.