MS211C/F - C Alculo Num Erico PROVA 1 (12/05/2022) Nome: RA
MS211C/F - C Alculo Num Erico PROVA 1 (12/05/2022) Nome: RA
MS211C/F - C Alculo Num Erico PROVA 1 (12/05/2022) Nome: RA
13 2
1. Seja p : R → R, dada por p(x) = −1 + 3
x − 2x4 .
(a) Partindo de x0 = 1, mostre que a sequência exata de iterandos gerada pelo método de Newton
não converge a um zero de p.
(b) Um estudante, duvidando do afirmado no item acima, resolveu fazer um experimento numérico
no Octave. As primeiras 21 iterações de Newton resultaram nos dados exibidos na tabela a
seguir.
k xk f (xk ) k xk f (xk ) k xk f (xk )
1 −1.00000 1.3333 · 100 8 1.00062 1.3337 · 100 15 1.67714 −4.6349 · 100
2 1.00000 1.3333 · 100 9 −1.02894 1.3460 · 100 16 1.47740 −1.0700 · 100
3 −1.00000 1.3333 · 100 10 5.61611 −1.8539 · 103 17 1.39505 −1.4173 · 10−1
4 1.00000 1.3333 · 100 11 4.26129 −5.8178 · 102 18 1.38033 −4.0848 · 10−3
5 −1.00000 1.3333 · 100 12 3.26184 −1.8130 · 102 19 1.37988 −3.7517 · 10−6
6 1.00000 1.3333 · 100 13 2.53481 −5.5725 · 101 20 1.37988 −3.1761 · 10−12
7 −1.00001 1.3333 · 100 14 2.02039 −1.6637 · 101 21 1.37988 1.7764 · 10−15
2 +2x
2. Seja f (x) = x − e−x .
(a) Prove que f tem pelo menos um zero. Exiba um intervalo finito que contenha um único zero
de f e tenha comprimento menor ou igual a 0.5.
(b) Aplique um dos métodos visto nesta disciplina para estimar o valor do zero de f , com pelo
menos 4 casas decimais corretas.
2 4 1 5
3. Considere o sistema Ax = b, com A = 5 −2 3 e b = −13 .
1 6 4 0
Abaixo está exibida a sequência de operações elementares utilizadas no escalonamento de A, com
pivoteamento parcial.
`2 2 4 1 5 −2.0 3.0 5 −2.0 3.0
`1 5 −2 3 ∼ `2 − 0.4`1 2 4.0 1.0 ∼ `3 0 4.8 −0.2
1 6 4 `3 − 0.2`1 1 6.0 4.0 `2 0 6.4 3.4
5 −2.0 3.0 5 −2.0 3.0
∼ 0 6.4 3.4 ∼ 0 6.4 3.4
`3 − 0.75`2 0 4.8 −0.2 0 0 −2.75
(a) Caso a matriz A fosse a seguir, seria possı́vel assegurar, de alguma forma, a convergência
da sequência gerada pelo método de Jacobi? Se sim, você pode dizer alguma coisa sobre a
velocidade com que o erro reduziria a cada passo?
10 3 1 5
2 5 1 1
A=
1
.
4 0 9
3 1 8 1
(b) Escreva um algoritmo para o método de Gauss-Seidel. Deixe claro quais outras informações
iniciais, além de A e b, seu algoritmo necessita. Qual o custo computacional, em termos do
número de operações aritméticas de ponto flutuante, de cada iteração do algoritmo?