Introdução ao Cálculo Numérico

Prévia do material em texto

�PAGE \* MERGEFORMAT�42�
CÁLCULO NUMÉRICO
Bibliografia
- RUGGIERO, M. A. G.; LOPES, V. L. R. Cálculo numérico: aspectos teóricos e computacionais. São Paulo: McGraw-Hill, 1997. 
- SPERANDIO, Décio; MENDES, João Teixeira e Silva; LUIZ H. M.. Cálculo Numérico Características matemáticas e computacionais. São Paulo: Prantice Hall, 2003
- CLAUDIO, Dalcídio Moraes; MARINS, Jussara Maria. Cálculo numérico computacional - teoria e prática. São Paulo : Atlas, 1989.
Introdução
Cálculo Numérico trata de métodos numéricos para resolução de problemas em diversas áreas. O Cálculo Numérico corresponde a um conjunto de ferramentas ou métodos usados para se obter a solução de problemas matemáticos de forma aproximada. Esses métodos se aplicam principalmente a problemas que não apresentam uma solução exata, portanto precisam ser resolvidos numericamente.
Para chegarmos a um resultado numérico, é necessário percorrer uma seqüência de passos preestabelecidos:
Levantamento de dados;
Construção de modelo matemático;
Escolha de método numérico;
Implementação Computacional;
Análise dos resultados;
Se necessário, reformular o modelo.
Em cada passo pode existir uma parcela de erro que se acumula ao final do processo. Os erros na fase de levantamento de dados ocorrem, por exemplo, por causa de aparelhos imprecisos ou leitura mal feita de resultados experimentais. Os erros na fase de modelagem são causados, geralmente por simplificações no modelo usado.
Assim, pode ocorrer que os resultados finais estejam distantes do que se esperaria obter ainda que todas as fases da resolução tenham sido realizadas corretamente. Os resultados dependem também:
Da precisão da máquina.
Da forma de representação do número.
Das operações numéricas efetuadas.
Por exemplo, na interação usuário com o computador os dados de entrada são enviados para o computador no sistema decimal e é convertido pelo computador em sistema binário e as operações são realizadas em binário, onde os resultados são convertidos novamente em decimal. Todo esse processo de conversão é fonte de erros.
Exemplos:
Calcule a área e o perímetro de uma circunferência de raio 100m, considerando:
( = 3,14
( = 3,1416
( = 3,141592654
É possível obter a área exata da circunferência? Qual o melhor resultado?
Calcule 1/3 *3 com uma máquina de bolso e com uma máquina científica.
1. ERROS
Os erros de truncamento e de arredondamento já conhecidos não podem ser eliminados ou mesmo minimizados, pois dependem do Sistema Operacional da máquina usada.
Erro absoluto e Erro Relativo
Seja x o valor exato de um número e 
 seu valor aproximado, calculado por algum método numérico.
O erro absoluto nessa aproximação é dado por Eax = x -
 .
Em geral não conhecemos o valor de x, apenas seu valor aproximado 
, e podemos apenas estimar o erro absoluto.
Exemplos:
Sabemos que ( ( (3,14; 3,15), logo o erro estimado será EA( < 0,01.
Considerando 
 = 2112,9 com EAx < 0,1 e 
=5,3 com EAy < 0,1, qual o intervalo de aproximação de x e de y? Observando que os limitantes superiores do erro são os mesmos podemos dizer que a precisão na representação é a mesma? 
Como os números são de ordens de grandezas diferentes precisamos de uma outra medida para compara-los; precisamos do Erro Relativo. ERx = EAx / x
No exemplo 2 temos, ERx = 0,1/2112,9 ( 4,7*10-5 e ERy = 0,1/ 5,3 ( 0,2*10-1. Logo, temos que x é representado com maior precisão.
QUANTO MENOR O ERRO RELATIVO MAIOR A PRECISÃO.
Erros nas operações
Temos que:	 EAx = x - 
( x = 
+EAx 
EAy = y - 
( y = 
+EAy 
Na adição: x+y = 
+
+EAx +EAy, ou seja, EAx+y = EAx + EAy.
Na multiplicação: x*y =(
+EAx)*(
+EAy) = 
�� EMBED Equation.3 +
.EAy +
EAx + EAxEAy
Como EAxEAy é pequeno, será desprezado.	Logo, EAx*y = 
.EAy +
EAx.
 
Observamos que o erro na multiplicação é maior que na adição.
Exemplos:
1) Considere uma máquina que permite até 4 dígitos, x = 0,3491*104 e y = 0,2789*100
Calcule (x+y) –x.
2) Considere máquina que permite até 8 dígitos, x = 3,00000053 e y = 2,00000021. Calcule x+y e x*y nessa máquina.
3) Com 3 dígitos significativos calcule 15,9*(4,99+0,02) e depois (15,9*4,99)+(15,9*0,02)
Conversão dos Números nos Sistemas Decimal e Binário
Números Inteiros
Considere os números 
 e 
. Cada um destes números está representado em uma base diferente, um na base 10 e o outro na base 2, respectivamente.
Estes números podem ser decompostos, de acordo com suas bases, da seguinte forma:
e
De modo geral, pode-se escrever um número na base 
, 
 com 
, na forma:
Desta forma, ficará mais simples analisar uma forma de converter um número do sistema binário para o sistema decimal. Vejamos como converter o número
:
Exemplo (3) Converter os números abaixo do sistema binário para o sistema decimal.
(a) 
(b) 
Vamos agora criar um algoritmo para esta conversão. Considere, por exemplo, o número 
, isto é, 
. Agora vamos colocar o número 2 em evidência:
Do mesmo modo como fizemos no exemplo acima, pode-se criar um algoritmo de forma genérica para esta conversão. De modo geral, a representação de um número 
 na base 10, que denotaremos por 
 é obtida da seguinte forma:
Assim, pelo algoritmo, a forma de conversão do número 
 para a base 10 será:
Agora, veremos um processo prático para converter números inteiros no sistema decimal para o sistema binário. Para isso, iremos dividir o número dado na base 10 sucessivamente por 2 até chegarmos em um quociente 0, e depois disso, para escrever o número na base 2, utilizaremos o resto destas divisões da última divisão para a primeira, ou seja, os restos serão tomados de traz para frente. 
Exemplo (4) Converter o número 
 da base 10 para a base 2, utilizando o processo prático.
Exemplo (5) Converter os números abaixo da base 10 para a base 2, utilizando o processo prático.
(a) 
(b) 
Agora, vamos criar um algoritmo. Considerando o processo prático descrito anteriormente e que 
. Desta forma, vamos admitir que:
Assim, temos que o resto desta divisão (237 por 2) é 1, e este valor será tomado como 
, ou seja, 
. Agora, iremos repetir este processo para o quociente 118, e este valor será o 
, ou seja, 
. Daí,
Desta forma, vemos que o resto da divisão é zero, assim temos que 
. Continuando o processo temos:
N2 = 	N3=	N4=	N5=	N6=	N7=
Portanto, 
De modo geral, ao considerarmos um número N qualquer na base 10, 
, e sua representação binária dada por 
, podemos utilizar o seguinte algoritmo, onde a cada k, se obtém o dígito binário 
. Vejamos:
Passo 0: Para 
, temos 
.
Passo 1: Obtenha 
 e 
 tais que 
.
 Faça 
Passo 2: Se 
, pare. 
 Caso contrário, faça 
.
 Faça 
 e volte para o Passo 1.
Números Fracionários
Números fracionários são números da forma:
, 
 ou 
Como se pode notar acima estes números podem ser finitos (
) ou infinitos (
 e 
). Primeiramente trataremos da conversão de números fracionários da base decimal para a base binária. 
Considere 
. Sua representação na base binária será dada por 
. Vamos tomar este número como exemplo para descobrir uma forma de conversão para o sistema binário.
Assim, 
Multiplicando toda a igualdade por 2, temos:
Daí, temos que 
, uma vez que a parte inteira do número 
 é zero. Agora, vamos aplicar este mesmo processo até que a parte fracionária seja igual a zero. Quando isso ocorrer, a conversão se encerra. Vejamos: 
Portanto, 
 . Agora, nem sempre este processo é finito. Veja o exemplo a seguir:
Exemplo (6) Converter o número 
 da base 10 para a base 2, utilizando o processo anterior.
Assim, 
Analisando o exemplo acima, podemos notar que a representação não finita de certos números pode causar erros aparentementeinexplicáveis nas máquinas, como foi o caso do exemplo (2). Este tipo de erro ocorre, pois as máquinas não têm a capacidade de guardar infinitos dígitos na mantissa.
Agora vamos fazer a conversão do sistema binário para o decimal. Para isso, considere 
. Sua representação no sistema decimal será obtida de forma muito semelhante ao que foi feito no processo de conversão para números inteiros. De modo geral, seguiremos o seguinte algoritmo: defina 
 e multiplique o número 
 por 10, mas na base binária, é dado por 
, para assim obtermos o dígito 
, que é a parte inteira deste produto, mas convertido para a base decimal. Vejamos este processo através de um exemplo.
Exemplo (7) Converter o número 
 da base binária para a base decimal.
Exemplo (8) Converter 
 da base 2 para a base 10. (Utilize 6 casas decimais).
Aritmética de Ponto Flutuante 
Uma máquina representa um número real em um sistema chamado aritmética de ponto flutuante, isto é, 		
Onde,
 base em que a máquina opera;
 número de dígitos da mantissa;
 com 
 e 
;
expoente no intervalo 
, com 
 
 e 
.
Em qualquer máquina, poucos números reais são representados exatamente, e, portanto, a representação de um número real será feita por truncamento ou arredondamento. Vejamos um exemplo:
Exemplo (9) Considere uma máquina que opera no sistema 
, 
 e 
. De que forma os números ficam representados nesta máquina (em módulo) ?
Considere o conjunto 
. Dado 
 pode ocorrer um dos seguintes casos:
I – 
: 
II – 
: 
III – 
: 
Observações:
Algumas linguagens de programação permitem que as variáveis sejam declaradas em precisão dupla, ou seja, uma variável representada no sistema de aritmética de ponto flutuante da máquina com aproximadamente o dobro de dígitos disponíveis na mantissa.
O zero em ponto flutuante é, em geral, representado com o menor expoente possível a máquina. Isto ocorre, pois a representação do zero por uma mantissa nula e um expoente qualquer para a base 
 pode acarretar perda de dígitos significativos no resultado da adição deste zero a u outro número.
Exemplo (10) Representar os números num sistema de aritmética de ponto flutuante onde 
, 
 e 
.
Solução:
	
	Representação obtida por arredondamento
	Representação obtida por truncamento
	1,25
	
	
	10,065
	
	
	-826,95
	
	
	2,71828...
	
	
	0,000008
	
	
	817235,92
	
	
Erros De Arredondamento e Truncamento em um Sistema de Aritmética de Ponto Flutuante
Seja um sistema que opera em aritmética de ponto flutuante de 
 dígitos, na base 10, e seja 
 escrito na forma:
onde, 
 e 
.
Exemplo (13) Se 
 e 
, represente 
 como acima.
A grande questão do exemplo anterior é como considerar a parcela de 
 na representação do número uma vez que esta parcela não pode ser incorporada totalmente na mantissa e definir o erro absoluto (ou relativo) máximo cometido. Para isso, temos dois critérios: o truncamento e o arredondamento.
No truncamento, 
 é desprezado e 
. Neste caso, temos:
No arredondamento, 
 é modificado para levar 
 em consideração. A forma de arredondamento mais utilizada é a seguinte:
Portanto, se 
, 
 é desprezado, caso contrário, soma-se 1 ao último dígito de 
.
Assim, se 
, temos:
E, se 
, teremos:
	e	 
.
A prova desta última desigualdade pode ser encontrada nos livros textos da bibliografia básica deste curso.
Observação: apesar de erros menor ocorrerem quando se utiliza arredondamento, este acarreta em um tempo de execução maior, e por isso, o truncamento é mais utilizado.
Análise de Erros nas Operações de Aritmética de Ponto Flutuante
Dada uma sequência de operações, como por exemplo, 
, é importante saber como o erro se propaga ao longo das operações seguintes.
O erro total de em uma operação é composto pelo erro das parcelas e pelo erro no resultado da operação. 
Para os exemplos 14, 15, 16 e 17, vamos considerar um sistema de aritmética de ponto flutuante com 
, 
 e com acumulador de precisão dupla.
Exemplo (14) Dados 
 e 
, obter os valores aproximados e exatos de:
a) 
		b) 
No exemplo acima se pode notar que mesmo que as parcelas de uma operação estejam representadas exatamente no sistema, o resultado final armazenado pode não ser exato.
Fórmulas de Erros nos Fatores
Na maioria dos sistemas, o resultado exato da operação denotado por OP é normalizado e depois arredondado ou truncado para 
 dígitos, obtendo assim o resultado aproximado 
 que é armazenado na memória da máquina.
Já vimos que o erro no resultado de uma operação, supondo que as parcelas estão exatamente representadas será:
 no truncamento
e
 no arredondamento.
Vamos agora ver as fórmulas para os erros absolutos e relativos nas operações de aritmética de ponto flutuante com erros nas parcelas. Estas fórmulas de erros não serão provadas, mas fica a critério do aluno olhar a demonstração nos livros da bibliografia.
Para estas fórmulas, vamos supor que erro final é arredondado.
Sejam 
 e 
, tais que 
 e 
.
ADIÇÃO: 
Erro absoluto: 
 e Erro relativo: 
.
SUBTRAÇÃO: 
Erro absoluto: 
 e Erro relativo: 
.
MULTIPLICAÇÃO: 
Erro absoluto: 
 e Erro relativo: 
.
DIVISÃO: 
Erro absoluto: 
 e Erro relativo: 
Note que o erro todas as fórmulas acima foram escritas sem considerar o erro de arredondamento (truncamento) no resultado final. E não se esqueça que a análise completa da propagação de erros se faz considerando os erros nas parcelas e no resultado de cada operação efetuada.
Exemplo (15) Suponha que 
, 
, 
, e 
 estejam representados exatamente. Qual é o erro no cálculo de 
? 
Exemplo (16) Sejam 
 e 
. Determine o erro relativo da operação 
, que é denominado cancelamento subtrativo.
Exemplo (17) Considerando o valor de 
 do exemplo acima, determine o valor de 
, se 
 e 
. Determine também o erro relativo desta operação.
EXERCÍCIOS
(1) Converter para decimal os seguintes números binários:
a) 10011	 b) 11100010		c) 1000001 	d) 1,1 	e) 1100,01	 f) 1000,001
(2) Converter para binário os seguintes números decimais:
a) 23 		b) 2615 	c) 2,5 		d) 0,1 		e) 3,8 		f) 10,05
(3) Seja um sistema de aritmética de ponto flutuante de quatro dígitos, base decimal e com acumulador de precisão dupla. Dados os números:
	
	e	
efetue as seguintes operações e obtenha o erro relativo no resultado, supondo que x, y e z estão exatamente representados:
a) 
	b) 
	c) 
	d) 
	e) 
(4) Suponha que x é representado no computador por 
, onde 
 é obtido por arredondamento, obtenha os limites superiores para os erros relativos de 
		e	
(5) Considere uma máquina cujo sistema de representação de números é definido por: base 10, (
), quatro dígitos na mantissa (
) e expoente no intervalo 
. Pede-se:
a) Qual o menor e o maior, em módulo, número representado nessa máquina?
b) Como será representado nessa máquina o número 73758 nessa máquina, se for usado o arredondamento? E se for usado truncamento?
c) Se a = 42450 e b = 3 qual o resultado de 
 se for usado o arredondamento? E se for usado truncamento? d) Repetir o item d) para a operação: 
.
2. ZEROS DE FUNÇÕES REAIS
Vamos ver métodos numéricos para resolução de equações não lineares. Um número real ( é um zero da função f(x) ou uma raiz da equação f(x) =0 se f(() =0.Graficamente, zeros reais são representados pelas abcissas dos pontos onde a curva intercepta o eixo x. 
Como obter raízes reais de uma equação qualquer?
Fora equações do segundo grau, é muito difícil encontrar as raízes de equações de grau mais alto, e às vezes nos contentamos com aproximações para esses zeros. Mas, com os métodos numéricos seremos capazes de a menos de limitação de máquina, encontrar os zeros de uma função com qualquer precisão determinada. A idéia é partirde uma aproximação inicial para a raiz e em seguida refinar essa aproximação através de um processo iterativo.
Observação: Iteração de modo geral significa repetição de um processo. A maioria dos métodos numéricos é iterativa. Elementos de métodos iterativos:
Tentativa inicial: primeira aproximação para a solução.
Equação de recorrência: equação por meio da qual, partindo-se da tentativa inicial, são realizadas as iterações ou aproximações sucessivas, para a solução desejada.
Teste de parada: a partir do qual o processo iterativo é finalizado.
Os métodos iterativos, assim, constam de duas fases:
Localização ou isolamento das raízes, onde encontramos um intervalo que contém a raiz.
Refinamento onde após escolhidas as aproximações iniciais no intervalo, melhoramos sucessivamente até obter uma aproximação para a raiz dentro de uma precisão ( prefixada.
FASE 1: ISOLAMENTO DE RAÍZES: Analise teórica ou gráfica de f(x)
Teorema: Seja f(x) uma função contínua num intervalo [a,b]. Se f(a). f(b) <0 então existe pelo menos um ponto x = ( entre a e b que é zero de f(x).
Observação: Se, além disso, existir f’(x) e f´(x) preservar sinal em (a,b), este intervalo contém um único zero. 
	Uma forma de isolar as raízes de f(x) usando resultados anteriores é tabelar f(x) para vários valores de x e analisar as mudanças de sinal de f(x) e o sinal da derivada nos intervalos em que f(x) mudou de sinal.
Exemplos:
1) f(x) = x³-9x+3
	X
	-(
	-100
	-10
	-5
	-3
	-1
	0
	1
	2
	3
	4
	5
	f(x)
	-(
	-
	-
	-
	+
	+
	+
	-
	-
	+
	+
	+
Sabendo que f(x) é contínua para qualquer x real e observando as variações de sinal, podemos concluir que cada um dos intervalos I1 = [-5,-3], I2 = [0,1] e I3 = [2,3] contém pelo menos um zero. Como f(x) é um polinômio de grau 3 podemos afirmar que cada intervalo contém um único zero de f(x). Assim, localizamos todas as raízes de f(x).
�
2) f(x) = x+ex
	X
	-2
	-1
	0
	1
	2
	3
	f(x)
	-
	-
	+
	+
	+
	+
3) f(x) =
	X
	-1
	0
	1
	2
	3
	f(x)
	-
	-
	+
	+
	+
�
Obs: Se f(a).f(b) > 0 então podemos ter várias situações no intervalo [a,b], e devemos fazer uma criteriosa analise gráfica de f(x) para obter aproximações para a raiz.
4) Faça você: f(x) = x³-8x+6
	x
	-5
	-4
	-3
	-2
	-1
	0
	1
	2
	3
	4
	5
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
Logo, as raízes estão nos intervalos: ______________________.
FASE 2: REFINAMENTO DAS RAÍZES: Os métodos iterativos para obter zeros da função fornece apenas uma aproximação da raiz exata. Vamos estudar cinco métodos: Bisseção, Falsa posição, ponto fixo e Newton e secante. 
Critério de parada
Quando parar as iterações? Quando a aproximação xk está suficientemente próxima da raiz exata?
Existem dois testes a serem considerados:
|x-(|< ( (que equivale a |b-a|<( após as reduções do intervalo.
|f(x)| <(( que pode ser |f(a)|<( e |f(b)| <(.
2.1. Método da Bisseção
 Seja a função f(x) contínua no intervalo [a,b] e tal que f(a).f(b) <0. Vamos supor que f(x) possui uma única raiz no intervalo (a,b). O Objetivo deste método é reduzir a amplitude do intervalo que contem a raiz até atingir a precisão requerida : b-a<(, usando para isto a sucessiva divisão de [a,b] ao meio.
Após cada divisão é escolhido um novo intervalo (a,b) que possui a raiz ( procurada. Para isto devemos avaliar os valores de f em a, b e no ponto médio m.
Logo para este método temos como entradas: extremos do intervalo (a,b), precisão ( solicitada para a aproximação. Como saídas a aproximação para o zero ( da função, ou seja, a solução.
Exemplos: 
1) Seja f(x) = x.log(x) –1, determine uma aproximação para o zero desta função no intervalo [2,3] com precisão ( = 10-2.
Solução:
	x
	1
	2
	3
	4
	f(x)
	-1
	-0,3979
	0,4314
	1,40824
Logo, existe um zero de f em (2,3).
1a iteração: x1 = (a+b)/2 = (2+3)/2 = 2,5
f(2) = -0,3979
f(2,5) = 2,5.log2,5 –1 = -0,00515
f(3) = 0,4314.			Zero procurado está em (2,5;3).
Critério de parada: |f(x1)| =|-0,00515 < (. Logo a solução será ( = 2,5.
Outro modo rápido de resolver esse exercício é usando uma tabela.
	a
	b
	f(a)
	f(b)
	xk
	f(xk)
	2
	3
	-0,3979
	0,4314
	2,5
	-0,00515< (
2) Sabendo que f(x)=x+ex, determine uma aproximação para o zero desta função no intervalo [-1, 0] com precisão ( = 10-2 
	a
	b
	f(a)
	f(b)
	xk
	f(xk)
	-1,000000
	0,000000
	-0,632121
	1,000000
	-0,500000
	0,106531
	-1,000000
	-0,500000
	-0,632121
	0,106531
	-0,750000
	-0,277633
	-0,750000
	-0,500000
	-0,277633
	0,106531
	-0,625000
	-0,089739
	-0,625000
	-0,500000
	-0,089739
	0,106531
	-0,562500
	0,007283< (
Logo, a solução será ( = -0,5625.
 
3)Calcular uma aproximação da raiz de f(x) = ex –sen(x)-2 pelo método da Bisecção com ( = 10-2.
	x
	-4
	-3
	-2
	-1
	0
	1
	2
	3
	4
	f(x)
	-2,738487
	-1,809093
	-0,955367
	-0,790650
	-1
	-0,123189
	4,479759
	17,94417
	53,35
Logo, a raiz está entre 1 e 2. 
a=1( f(a) = f(1) = e1-sen(1)-2 = -0,123189
b=2( f(b) = f(2) = e2-sen(2)-2 = 4,479759
1a iteração: x1 = (a+b)/2 = (1+2)/2 = 1,5
f(1) = -0,123189
f(1,5) = 1,484194		a= 1 e b = 1,5
f(2) = 4,479759
Critério de parada: |f(x1)|=1,484194 >(
2a iteração: x2 = (a+b)/2 = (1+1,5)/2 =1,25
f(1) = -0,123189
f(1,25) = 0,541358
f(1,5) = 1,484194		a= 1 e b = 1,25
Critério de parada: |f(x2)|= 0,541358> (
3a iteração: X3 = (a+b)/2 = (1+1,25)/2 =1,125
f(1) = -0,123189
f(1,125) =0,177949		a= 1 e b = 1,125
f(1,25) = 0,541358
Critério de parada: |f(x3)|= 0,177949 > (
4a iteração: X4 = (a+b)/2 = (1+1,125)/2 =1,0625
f(1) = -0,123189
f(1,0625) = 0,020021
f(1,125) =0,177949		a= 1 e b = 1,0625
Critério de parada: |f(x4)|= 0,020021 > (
5a iteração: X5 = (a+b)/2 = (1+1,0625)/2 =1,03125
f(1) = -0,123189
f(1,03125) =- 0,053372		a = 1,03125 e b = 1,0625
f(1,0625) = 0,020021
Critério de parada: |f(x5)|= 0,053372 > (
6a iteração: X6 = (a+b)/2 = (1,03125+1,0625)/2 =1,046875
f(1,03125) =-0,053372
f(1,046875) = -0,017129
f(1,0625) = 0,020021		a= 1,046875 e b = 1,0625
Critério de parada: |f(x6)|= 0,017129 > (
7a iteração: X7 = (a+b)/2 = (1,046875+1,0625)/2 =1,0546875
f(1,046875) = -0,017129
f(1,0546875) = 0,001332
f(1,0625) = 0,020021		a= 1,046875 e b = 1,0546875
Critério de parada: |f(x7)| = 0,001332 <(
Logo, a raiz será ( = 1,0546875.
4) Dada a função f(x) = x³+3x-1 pesquise a existência de raízes reais em intervalos. E determine uma aproximação para cada raiz real com precisão e = 0,1.
5) Determine uma aproximação para a raiz da função f(x) = x² - sen(x), com 4 iterações. Calcule o erro cometido.
Vantagens: Envolve operações aritméticas simples. É possível estimar o número de iterações. Desvantagens: convergência lenta.
2.2. Método da Falsa Posição
Seja f(x) contínua no intervalo [a, b] e tal que f(a).f(b) <0. Queremos encontrar a raiz de uma função f(x), ou seja, encontrar o xi na qual f(xi) mais se aproxime de zero. Precisamos, primeiro, localizar a raiz, ou seja, definir o intervalo onde ela se encontra: [a,b]. Podemos ter a raiz mais próxima de a do que de b, ou vice versa, logo, uma vez definido o intervalo, traçamos uma reta que passe por x1= a e x2 = b. 
Esta reta irá cortar o eixo X em um ponto, x3, que é uma aproximação da raiz. Assim, em vez de tomar a média aritmética entre a e b, o método da falsa posição toma a média ponderada entre a e b com pesos |f(b)| e |f(a)|, respectivamente:
 ou 
 
	Repete-se o processo até que | f(xk+1) | < e. 
	O Método da Falsa Posição é um excelente método quando o intervalo escolhido para a raiz é pouco preciso. Em último caso podemos simplesmente escolher dois pontos x que tenham f(x) de sinais opostos e aplicar o método. O método pode ser aplicado a qualquer função e, assim como o Métododa Bisseção, a raiz irá sempre convergir (desde que considerado um intervalo contínuo). Além disso, o método utiliza apenas matemática elementar, não exigindo conhecimentos de integração ou derivação. Em contrapartida, sua convergência é lenta.
Exemplos: Usando o método da falsa posição:
1) Ache a raiz de f(x) = sen (x)-x/2 , no intervalo [π/2, π]com e = 10-2 .
2) Ache as raízes de f(x) = x² - 3, com e = 10-2 . 
3) Determine o zero de f(x) = x.log(x) -1 em [2,3] com 5 iterações.
4) No intervalo [0,1] e e = 5x10-4, determine o zero de f(x) = x³-9x+3.
	
	Exercícios Resolvidos usando 
1) sen(x)-x/2 =0 [pi/2; pi] e = 0,01
	
	
	k
	a
	f(a)
	b
	f(b)
	x
	f(x)
	
	1
	1,570796
	0,214602
	3,141593
	-1,570796
	1,759603
	0,102427
	
	2
	1,759603
	0,102427
	3,141593
	-1,570796
	1,844202
	0,040756
	
	3
	1,844202
	0,040756
	3,141593
	-1,570796
	1,877013
	0,014975
	
	4
	1,877013
	0,014974
	3,141593
	-1,570796
	1,888954
	0,005336< (
	
	
	
	
	
	
	Logo x = 1,888954 é solução
	2)f(x) = x² - 3, com e = 10-2 .
	
	
	
	
	
	x
	0
	1
	2
	3
	
	
	
	f(x)
	-3
	-2
	1
	6
	
	raiz em [1,2]
	
	
	
	
	
	
	
	
	k
	a
	f(a)
	b
	f(b)
	x
	f(x)
	
	1
	1,000000
	-2,000000
	2,000000
	1,000000
	1,666667
	-0,222222
	
	2
	1,666667
	-0,222221
	2,000000
	1,000000
	1,727273
	-0,016529
	
	3
	1,727273
	-0,016528
	2,000000
	1,000000
	1,731707
	-0,001190< (
	
	
	
	
	
	
	Logo x = 1,731707 é solução
	3) f(x) = x.log(x) -1 em [2,3] com 5 iterações.
	
	
	
	k
	a
	f(a)
	b
	f(b)
	x
	f(x)
	
	1
	2,000000
	-0,397940
	3,000000
	0,431364
	2,479848
	-0,021886
	
	2
	2,479848
	-0,021886
	3,000000
	0,431364
	2,504964
	-0,001016
	
	3
	2,504964
	-0,001017
	3,000000
	0,431364
	2,506128
	-0,000047
	
	4
	2,506128
	-0,000047
	3,000000
	0,431364
	2,506182
	-0,000002
	
	5
	2,506182
	-0,000002
	3,000000
	0,431364
	2,506184
	-0,0000001
	< 10-6
	
	erro ( sempre positivo
	
	Logo x = 2,506184 é solução
	
4) [0,1] e = 5x10-4, f(x) = x³-9x+3.
	
	
	
	
	k
	a
	f(a)
	b
	f(b)
	x
	f(x)
	
	1
	0,000000
	3,000000
	1,000000
	-5,000000
	0,375000
	-0,322266
	
	2
	0,000000
	3,000000
	0,375000
	-0,322266
	0,338624
	-0,008790
	
	3
	0,000000
	3,000000
	0,338624
	-0,008787
	0,337635
	-0,000226< (
	
	
	
	
	
	
	Logo x = 0,337635 é solução
2.3. Método do Ponto Fixo (MPF)
Seja f(x) função contínua em [a,b], intervalo que contém uma raiz da equação f(x) = 0.
O MPF consiste em transformar f(x) =0 em uma equação equivalente x = ((x) e a partir de uma aproximação inicial x0 gerar a seqüência {xk} de aproximações para ( (raiz) pela relação xk+1 = ((xk), pois a função ((x) é tal que f(() = 0 ( ((() = (. Logo, transformamos o problema de determinar um zero da função em determinar um ponto fixo da função. A função ((x) é chamada de função de iteração para a equação f(x)=0.
Exemplo: Para a equação x²+x-2 = 0 temos várias funções de iteração:
x = 2- x²( (1 (x) = 2 – x²
x = 
( (2 (x) =
x =
( (3 (x) = 
x = 
( (4(x) = 
A forma geral das funções de iteração ((x) é ((x) = x+A(x).f(x), onde A(x) é uma função tal que A(() (0. Logo, vemos que dada equação f(x) =0 existem infinitas funções de iteração ((x). Podemos ver, graficamente, que para certas ((x) o processo pode gerar uma seqüência que diverge de (
Exemplo: Seja f(x) = x² +x-6 e tomemos (1(x) = 6-x² e (2(x) = 
 e escolhemos x0 =1,5 e vamos tentar encontrar uma aproximação para o zero de f(x), ( =2.
O processo iterativo para (1(x) = 6-x² é xk+1 = (1(xk)= 6 – xk2, k = 0,1,2,3,...
x0 =1,5
x1= (1(x0) = 6-x02 = 6-1,52 = 3,75
x2= (1(x1) = 6-x12 = 6-3,752 = -8,0625
x3= (1(x2) = 6-x22 = 6-(-8,0625)2 = -59,0039
E a seqüência não está convergindo.
O processo iterativo para (1(x) =
 é xk+1 = (2(xk)= 
; k = 0,1,2,3,...
x0 =1,5
x1= (2(x0) =
= 
=2,1213
x2= (2(x1) = 
= 
 = 1,9694
x3= (2(x2) = 
=
 = 2,0076
x4= (2(x3) = 
=
 = 1,9981
x5= (2(x4) = 
=
 = 2,0005
E a solução está convergindo para a solução exata ( = 2.
Assim, vemos que nem toda função de iteração gera uma seqüência convergente para a solução procurada. O teorema a seguir fornece condições suficientes para que o processo seja convergente.
Teorema: Seja ( uma raiz da equação f(x) =0, isolada num intervalo I centrado em (. Seja ((x) uma função de iteração para a equação f(x) =0. Se
((x) e (’(x) são contínuas em I,
|(’(x)| < 1, para todo x em I,
x0 ( I.
Então a seqüência xk gerada pelo processo xk+1 = ((xk) converge para (.
Exercício: Verifique a convergência das funções de iteração (1(x) = 6-x² e (2(x) = 
 para a função f(x) = x²+x-6.
Temos que quanto menor for |(’(x)| mais rápida será a convergência. Então, basta tomarmos (’(() = 0. Assim, podemos construir uma função ((x) para a qual a convergência da seqüência xk está garantida, desde que x(I = (a,b). 
((x) = x+A(x).f(x)
(’(x) = 1+A’(x).f(x)+A(x).f ’(x)
(’(() = 1+A’(().f(()+A(().f ’((), como f(() = 0
(’(() = 1+A(().f ’((), fazendo (’(() = 0, temos 
1+A(().f ’(() = 0 ( A(() = 
. 
Para satisfazer isso, basta escolhermos A(x) = 
.
Logo, ((x) = x+A(x).f(x) =
, que será a função de iteração do método de Newton.
2.4. Método de Newton- Raphson
No método de Newton, na tentativa de garantir e acelerar a convergência do MPF, escolhemos a função de iteração ((x) tal que (`(() = 0, ou seja, escolhemos a função de iteração ((x) =
, com o processo iterativo xk+1 = 
, k = 0,1,2,3...
Critério de parada
Seja xk+1 = ((xk) a partir de uma aproximação inicial x0 dada. Se desejamos uma precisão (>0 na aproximação obtida o ideal é que |xk+1 -(|<(. Mas, como ( em geral é desconhecido não podemos fazer este cálculo. Assim, adotamos como critério de parada o cálculo da distância relativa entre dois termos consecutivos da sequência de aproximações: |f(xk+1)| < (.
Exemplos: Encontre um zero de f(x) = x²+x-6 pelo método de Newton com precisão ( = 0,01 a partir da aproximação inicial x0 = 1,5.
2.5. Método da Secante
Ao estudarmos o Método de Newton-Rapshon podemos notar que ele possui a desvantagem, de ter que obter 
 e calcular seu valor numérico a cada iteração.
Mas, existe uma forma de se contornar este problema, basta substituir a derivada 
 pelo quociente das diferenças:
onde 
 e 
 são duas aproximações para a raiz.
Neste caso, a função de iteração fica:
E, simplificando, temos:
, ou seja, 
Note que para aplicar este método são necessárias duas aproximações iniciais.
Interpretação Geométrica
Exemplos:
1) Aplicar o Método da Secante para encontrar a raiz da função 
, com 
, 
 e 
 e, sabendo-se que 
.
2) Aplicar o Método da Secante para encontrar a raiz da função 
, com 
, 
 e 
.
3) Resolva a equação x³-5=0 com x0 = 1,5 e x1 = 1,8 e precisão ( = 0,01pelo método da secante.
Como vimos, o Método da Secante é uma aproximação para o Método de Newton-Raphson, e por isso as condições de convergência são praticamente as mesmas, acrescentando apenas que este método pode divergir se 
. E, ainda podemos salientar que a ordem de convergência do Método da Secante não é quadrática como a do Método de Newton-Raphson, mas também não é linear. É possível provar que a ordem de convergência para este método é aproximadamente 
Alguns exercícios Resolvidos
1) Seja f(x) = 
, determine uma aproximação para o zero desta função no intervalo [1,2] com precisão ( = 0,5.
	x
	0
	1
	2
	3
	f(x)
	-5
	-0,8394
	0,737537
	1,483115
Raiz em (1,2)
1a iteração: x1 = (a+b)/2 = (1+2)/2 = 1,5
f(1) = -0,839397
f(1,5) = 0,109094		a= 1 e b = 1,5
f(2) = 0,737537
Critério de parada: |f(x1)|=0,109094 >(
2a iteração: x2 = (a+b)/2 = (1+1,5)/2 =1,25
f(1) = -0,839397
f(1,25) = -0,314490a= 1,25 e b = 1,5
f(1,5) = 0,109094
Critério de parada: |f(x2)|=0,314490 < (
Logo a solução será ( =1,25.
Pelo método de Newton:
F(x) = x1/2 – 5e-x
f’(x) = ½ x-1/2 -5e-x(-1) = 
1ª iteração: k = 0
f(1,5) = 0,109094
f’(1,5) = 1,523899
x1 = x0 – f(x0) / f’(x0) = 1,5 – 0,109094/1,523899 = 1,428411
Critério de Parada: |f(x1)| = |f(1,428411)| = 0,003286 <(
Logo, a raiz procurada é α= 1,428411.
2) Calcular uma aproximação da raiz de f(x) = x³+x²+x pelo método da Bissecção com ( = 0,2 no intervalo [-1;0,5].
a=-1( f(a) = f(-1) = (-1)³+(-1)²+(-1)=-1
b=0,5( f(b) = f(0,5) = (0,5)³+(0,5)²+0,5 =0,875 
1a iteração: x1 = (a+b)/2 = (-1+0,5)/2 = -0,25
f(-1) = -1
f(-0,25) =-0,203125		a= -0,25 e b = 0,5
f(0,5) = 0,875		Critério de parada: |f(x1)| >(
2a iteração: x1 = (a+b)/2 = (-0,25+0,5)/2 = 0,125
f(-0,25) = -0,203125
f(0,125) =0,142578		a= -0,25 e b = 0,125
f(0,5) = 0,875		Critério de parada: |f(x)| = < (
Logo, α = 0,125 é a solução procurada.
3) Obter a raiz quadrada do número 7, usando o método de Newton com precisão ( =0,01.
f(x) = x²-7	( f ’(x) =2x		2<
<3 ( x0 = 2,5 ou faz tabela
1ª iteração: k = 0	f(2,5) = -0,75; 	f’(2,5) = 5
x1 = x0 – f(x0) / f’(x0) = 2,5 –(-0,75)/5 = 2,65 
Critério de Parada:|f(x1)| = |f(2,65)|= 0,0225 > (
2ª iteração: k = 1	f(2,65) = 0,0225		f’(2,65) = 5,3
x2 = x1 – f(x1) / f’(x1) = 2,65 –0,0225/5,3 = 2,645755 
Critério de Parada: |f(x2)| = |f(2,645755)| = 0,000019 < (
Logo, a raiz procurada é α= 2,645755.
5) Determine uma raiz de f(x) = 2x-cos(x) no intervalo [0,(/4] pelo método de Newton, com x0 = (/8 e ( = 0,1.
f(x) = 2x-cos(x)
f’(x) =2+sen(x)	
1ª iteração: k = 0
f(0,392699) = -0,138481
f’(0,392699) = 2,382683
x1 = x0 – f(x0) / f’(x0) = 0,392699 -(-0,138481)/2,382683 = 0,450819
Critério de Parada|f(x1)| = 0,138481 >(
2ª iteração: k = 1
f(0,450819) = 0,001547
f’(0,450819) = 2,435703
x2 = x1 – f(x1) / f’(x1) = 0,450819-0,001547/2,435703= 0,450184
Critério de Parada: |f(x2)| = |f(0,450184)| = 0,000001 < (
Logo, a raiz procurada é α= 0,450184.
6) Determine uma raiz de f(x) = x5-6 pelo método de Newton com x0 = 1,5 e ( = 0,1.
f(x) = x5-6		f ’(x) =5x4	
1ª iteração: k = 0
f(1,5) = 1,55 -6= 1,593750
f’(1,5) = 5.1,54 = 25,3125
x1 = x0 – f(x0) / f’(x0) = 1,5 -1,593750/25,3125 = 1,437037
Critério de Parada: |f(x1)|=1,593750 > (
2ª iteração: k = 1
f(1,437037) = 0,128296
f’(2,65) = 21,322681
x2 = x1 – f(x1) / f’(x1) = 1,437037 –0,128296/21,322681 = 1,431020 
Critério de Parada: |f(x2)| = |f(1,431020)| = 0,001068 < (
Logo, a raiz procurada é α= 1,431020.
7) Determine uma raiz de f(x) = xlnx-1 pelo método de Newton com x0 = 1,5 e ( = 10-7.
Por termos o erro 10-7 devemos trabalhar com 8 casas decimais.
f(x) = xlnx-1
f’(x) = 1.ln(x)+x.1/x = lnx+1
f(1,5) = 1,5.ln(1,5)-1 = -0,39180234
f’(1,5) = ln(1,5)+1 = 1,40546511
x1 = x0 – f(x0) / f’(x0) = 1,5 – (-0,39180234)/1,40546511 = 1,77877059
Critério de Parada: |f(x1)| > (
2ª iteração: k = 1
f(1,77877059) = 0,02443391
f’(1,77877059) = 1,57592268
x2 = x1 – f(x1) / f’(x1) = 1,77877059-0,02443391/1,57592268 = 1,76326608
Critério de Parada: |f(x2)|>(
3ª iteração: k = 2
f(1,76326608) = 0,00006777
f’(1,76326608) = 1,56716782
x3 = x2 – f(x2) / f’(x2) = 1,76326608-0,00006777/1,56716782= 1,76322283
Critério de Parada: |f(x3)| = |f(1,76322283)| = 0,00000001 < (
Logo, a raiz procurada é α= 1,76322283.
8) A equação x³-2x-1=0 tem raízes reais.
Determine o intervalo onde está a raiz positiva e onde está a raiz negativa.
	X
	-5
	-4
	-3
	-2
	-1
	0
	1
	2
	3
	4
	5
	F(x)
	-116
	-57
	-22
	-5
	0
	-1
	-2
	3
	20
	55
	114
α=-1 é uma raiz negativa e a raiz positiva, está no intervalo [1,2].
Determine uma aproximação para a raiz em [1,2] com 3 iterações usando o método da Bisseção. Estime o erro.
f(1) = 1³-2.1-1 = -2 		f(2) = 2³-2*2-1 =3
k=1( x1 = (1+2)/2 = 1,5
f(1,5) = 1,5³-2.1,5-1 =-0,625
a= 1,5 e b = 2
k=2( x2 = (1,5+2)/2 = 1,75
f(1,75) = 1,75³-2.1,75-1=0,859375
a=1,5 e b = 1,75
 
k=3( x3 = (1,5+1,75)/2 = 1,625
f(1,625) =1,625³-2.1,625-1 =0,041016
a= 1,5 e b = 1,625
CP: |f(x3)| =0,041016 < (= 10-1.
Logo, α = 1,625 é a raiz procurada.
Usando a aproximação encontrada no item b, como x0, determine a raiz pelo método de Newton, com (=10-6.
f(x) = x³-2x-1
f’(x)=3x²-2
k=0(
f(x0)= f(1,593750) =1,593750³-2.1,593750-1= -0,139313
f’(x0) =f ’(1,593750) = 3.1,593750²-2 = 5,620117
x1 = x0 – f(x0)/f’(x0) = 1,593750- (-0,139313/5,620117) =1,618538
cp: |f(x1)| = |1,618538 ³-2.1,618538-1|=0,002952>10-6.
k=1(
f(x1)= f(1,618538) =1,618538³-2.1,618538-1= 0,002952
f ’(x1) =f ’(1,618538) = 3.1,618538²-2 = 5,858996
x2 = x1– f(x1)/f’(x1) =1,618538 - (0,002952/5,858996) =1,618034
cp: |f(x2)| = |1,618034 ³-2.1,618034-1|=0,0000001<10-6.
Logo, α = 1,618034 é a raiz aproximada.
Aplique o método da secante com x0 =1,5 e x1 = 1,7 comm 3 iterações.
f(x0)= f(1,5) =1,5³-2.1,5-1= -0,625
f(x1)= f(1,7) =1,7³-2.1,7-1= 0,513
k=1 ( x2 = 
=
=1,609842
f(x1)= f(1,7) =1,7³-2.1,7-1= 0,513
f(x2) =f(1,609842) = 1,609842³-2.1,609842 -1 = -0,047632
k=2 ( x3 = 
=
=1,617502
f(x3) =f(1,617502) = 1,617502³-2.1,617502 -1 = -0,003113
k=3 ( x4 = 
=
=1,618038
CP: |f(x4)| = |1,617502³-2.1,617502 -1| =0,000023 < 10-3.
Logo, ( =10-4 e α = 1,608038 é a raiz procurada.
Determine 3 funções de iteração, use ((x) =
 com 3 iterações para estimar a raiz, use x0 =1,5.
x³ -2x-1 = 0
x³ = 2x+1
x = 
(1(x) = 
x³ -2x-1 = 0
-2x = -x³+1
2x = x³-1
X = 
( (2 (x)= 
X³-2x-1=0
X(x²-2)-1=0
X(x²-2)= 1
X =
 
(3 (x)= 
Determinando a raiz, usando x0:
X1 =((x0) =((1,5)= 
=1,587401
X2 =((x1) =((1,587401)= 
=1,610196
X3 =((x2) = ((1,610196)=
=1,613036
Cp: f(1,613036) = -0,01168 ( |f(1,613036)|= 0,01168 < 10-1
α= 1,613036 é a raiz procurada, com precisão (=10-1.
1ª Lista de exercícios
1) Determine uma raiz de f(x) = 2x-cos(x) no intervalo [0,(/4] pelo método de Newton, com x0 = (/8 e ( = 0,1.
2) Determine uma raiz de f(x) = x5-6 pelo método de Newton com x0 = 1,5 e ( = 0,1.
3) Determine uma raiz de f(x) = xlnx-1 pelo método de Newton com x0 = 2,5 e ( = 10-7.
4) Encontre uma raiz aproximada para as funções a seguir utilizando o Método da Bissecção com precisão 
. Utilize a estimativa do número de intervalos para saber o número mínimo de iterações necessárias.
a) 
 e 
		b) 
 e 
.
5) Considere a função 
. No o intervalo 
 aplique o método da Bissecção até atingir a precisão 
. 
6) Usando o Método de Newton-Raphson encontre a raiz da função 
, com 
 e precisão
.
7) Dado o polinômio 
 tem suas raízes reais no intervalo (-1,1). Encontre pelo Método de Newton uma raiz aproximada com precisão 
 e 
.
8) Aplique o Método de Newton à equação 
 com 
. Justifique o que acontece.
9) Use o Método de Newton para obter a menor raiz positiva das equações a seguir com precisão 
. Para encontrar a menor raiz utilize os processos de localização de raízes (gráficos ou pelos teoremas).
a) 
	b) 
	c) 
10) Use o Método da Secante para obter a menor raiz positiva das equações a seguir com precisão 
. Para encontrar a menor raiz utilize os processos de localização de raízes (gráficos ou pelos teoremas).
a) 
	b) 
	c) 
11) Usando o Método da Secante encontre a raiz da função 
, com 
, 
 e precisão
.
12) A equação 2x³-7x²+2x+2=0 tem 3 raízes reais.
Determine o intervalo onde estão essas raízes. 
	X
	-3
	-2
	-1
	0
	1
	2
	3
	4
	5
	F(x)
	
	
	
	
	
	
	
	
	
Determine uma aproximação para a raiz em [0,1] com 5 iterações usando o método da Bisseção. Estime o erro.
Usando a aproximação encontrada no item b, como x0, determine a raiz pelo método de Newton, com (=10-6, (use7 casas decimais após a vírgula).
Aplique o método da secante com x0 =3 e x1 = 3,5 com 3 iterações.
13) Sabendo que a função tem uma raiz exata igual a , determine uma aproximação de com pelo menos 5 casas decimais significativas, utilizando o método da bissecção. Utilizar Ɛ < 0,01.
14) Seja a função :
Complete a tabela e encontre um intervalo onde se encontra uma raiz de f.
	x
	1
	2
	3
	4
	f(x) 
	
	
	
	
Partindo desse intervalo, utilizar o método de Newton para determinar o valor dessa raiz após 5 iterações.
A equação tem como raiz exata . Aplicar o método da falsa posição para determinar uma aproximação de . Utilizar Ɛ < 0,01.
Seja . Considerando o método do ponto fixo, verificar a convergência ou divergência da seguinte função de iteração de f(x) no intervalo
I = [0; 1[. ϕ(x) = . Se ϕ(x) for convergente no intervalo I, determinar uma aproximação da raiz de f nesse intervalo, utilizando Ɛ < 0,005.
Utilize o método de Newton para calcular o valor aproximado de . Considere Ɛ < 0,005
3. Interpolação Polinomial
Dado um conjunto de dados {xi,f(xi)} tal como na tabela abaixo:
	xi
	0
	1,5
	3,0
	4,5
	6,0
	f(xi)
	0,001
	0,016
	0,028
	0,046
	0,057
Como obter o valor de f(x) para um valor de x que não tenha sido medido, como por exemplo, x =2.0?
Quando se deseja saber o valor de f(x) para um x intermediário entre duas medidas, isto é, xi<x<xi+1, pode-se usar as técnicas da interpolação. 
A interpolação consiste em determinar uma função, que assume valores conhecidos em certos pontos (nós de interpolação), ou seja, aproximar essa função por uma outra função g(x), escolhida entre uma classe de funções definida a priori e que satisfaça algumas propriedades. A função g(x) é então usada em substituição à função f(x).
A necessidade de obter um valor intermediário que não consta de uma tabela ocorre comumente. Dados experimentais, tabelas estatísticas e de funções complexas são exemplos desta situação, e também:
Quando são conhecidos somente os valores numéricos da função para um conjunto de pontos e é necessário calcular o valor da função em um ponto não tabelado.
Quando a função em estudo tem uma expressão tal que operações como a diferenciação e a integração são difíceis.
Quando a função f(x) não é conhecida explicitamente.
A classe de funções escolhida para a interpolação é a priori arbitrária, e deve ser adequada às características que pretendemos que a função possua. 
Se a função a ser considerada é um polinômio, então teremos Interpolação Polinomial. Mas podemos escolher funções racionais, trigonométricas, dependendo das propriedades da função.
O problema geral da interpolação por meio de polinômios consiste em: Interpolar um ponto x a um conjunto de n +1 dados {xi,f(xi)}, ou seja, calcular o valor de f(x), sem conhecer a forma analítica de f(x) ou ajustar uma função analítica aos dados.Assim, obter um polinômio p(x) que passe por todos os pontos do conjunto de (n+1) dados {xi,f(xi)}, isto é: 
	p(x0)=f(x0)
	p(x1)=f(x1)
	…
	p(xn)=f(xn)
Obs: contagem começa em zero, portanto tem-se n+1 pontos na expressão, ou seja o polinômio tem grau n: pn(x) = a0+a1x+a2x²+... + anxn
 
O polinômio pn(x) é chamado de polinômio interpolador. Pode-se demonstrar que existe um único polinômio p(x) de grau menor ou igual a n que passa por todos os (n+1) pontos do conjunto {xi,f(xi)}. O conjunto de equações corresponde a um sistema linear de n+1 equações e n+1 variáveis.
	
Quais são as variáveis independentes? ai ou xi ?
A matriz dos coeficientes A é dada por 
é uma matriz de Vandermonde e, portanto desde que xi, i = 0,1,2,...n sejam pontos distintos, temos que o sistema linear acima admite solução única. Ou seja:
Teorema: Existe um único polinômio pn(x) de grau menor ou igual a n, tal que pn(xk) = f(xk), k = 0,1,2,...,n, desde que xk (xj, k(j.
Formas de se obter o polinômio
Uma das formas é resolver o sistema linear, porém, a determinação dos coeficientes do polinômio interpolador por meio da resolução de um sistema de equações lineares, apesar de ser conceitualmente simples, requer um certo esforço computacional.
Outras formas são: a forma de Lagrange e a forma de Newton 
Exemplo: Seja f(x) dada pela tabela abaixo. Encontre o polinômio interpolador e estime f(3).
	x
	0
	1
	2
	f(x)
	3
	6
	9
Nem sempre a resolução do sistema linear é simples e exata, pois a matriz dos coeficientes é uma matriz de Vandermonde, podendo ser mal condicionada.
Exemplo:
	x
	0,1
	0,2
	0,3
	0,4
	f(x)
	5
	13
	-4
	-8
Resolvendo o sistema linear por eliminação de Gauss e 3 casas decimais obtemos p3(x) = 6330x³-5050x²+1150x-66, onde vemos que para x = 0,4 , p3(0,4) = -10 ( -8 = f(0,4).
Forma de Lagrange
Seja x0,x1,x2,..., xn, n+1 pontos distintos e yi = f(xi), i = 0,1,2,..., n. Seja pn(x) o polinômio de grau menor ou igual a n, que interpola f em x0, x1, x2, ....xn. Podemos representar pn(x) na forma pn(x) = y0L0(x) +y1L1(x)+...+ynLn(x) onde os polinômios Lk(x) são de grau n. Para cada i, queremos que a condição pn(xi) = yi seja satisfeita. A forma mais simples de se satisfazer esta condição é impor:
 Lk (xi) = 
 e para isso, definimos Lk(x) por:
Lk (x)=
É fácil verificar que realmente Lk(xk) =1 e Lk(xi) =0 para k(i.	
Como o numerador de Lk(x) é um produto de n fatores da forma (x – xi) , temos que Lk(x) é um polinômio de grau n e, assim pn(x) é um polinômio de grau menor ou igual a n. Além disso, para x = xi temos: pn(x) =
Então a forma de Lagrange para o polinômio interpolador é pn(x) =
.
Exemplos: Dados os pontos abaixo, calcular o polinômio interpolador usando a fórmula de Lagrange.
	X
	-1
	0
	3
	F(x)
	15
	8
	-1
a)
P2(x) = y0.L0+y1.L1+ y2.L2
P2(x) = 15. 
 +8. 
-1. 
P2(x) = 
P2(x) = x²-6x+8
	X
	-1
	0
	3
	5
	F(x)
	15
	8
	-1
	10
b)
P3(x) = y0.L0+y1.L1+ y2.L2 + y3.L3
P3(x) = 15. 
+8. 
 -1. 
+10. 
P3(x) = 
 
P3(x) = 
x³+
x²-
x+8
Interpolação de Newton
Considerando os n+1 pontos (x0, f(x0)), ..., (xn,f(xn)) e o polinômio interpolador pn(x). Newton propôs de representar o polinômio pn(x) da forma:
Pn(x) = f[x0] + f[x0,x1] (x-x0)+ f[x0,x1, x2]. (x-x0).(x-x1)+….+ f[x0,x1, x2, …., xn]. (x-x0).(x-x1)…(x-xn-1)
Os coeficientes dk, k=0,...,n são diferenças divididas de ordem k entre os pontos (xj ,f(xj)), j=0,...,k
Operador diferenças divididas: dj
Sendo f(x) é uma função tabelada em x0,...,xn. Os operadores de diferenças divididas são definidos por:
Onde usamos, o seguinte esquema prático
	X
	Ordem 0
	Ordem 1
	Ordem 2
	...
	Ordem n
	x0
	f[x0]
	
	
	
	
	
	
	f[x0,x1]
	
	
	
	x1
	f[x1]
	
	f[x0,x1,x2]
	
	
	
	
	f[x1,x2]
	
	
	
	x2
	f[x2]
	
	f[x1,x2,x3]
	
	
	
	
	
	
	
		f[x0,...,xn]
	
	
	
	f[xn-2, xn-1, xn]
	
	
	....
	
	f[xn-1, xn]
	
	
	
	xn
	f[xn]
	
	
	
	
Ou mais completo
Exemplos: Dados os pontos abaixo, calcular o polinômio interpolador usando a fórmula de Newton.
	X
	-1
	0
	3
	F(x)
	15
	8
	-1
	X
	Ordem 0
	Ordem 1
	Ordem 2
	-1
	15
	
	
	
	
	-7
	
	0
	8
	
	1
	
	
	-3
	
	3
	-1
	
	
P2(x) = d0 + d1 (x-x0)+ d2 (x-x0).(x-x1)
P2(x) = 15-7(x+1)+1(x+1)(x-0)
P2(x) = 15-7x-7+x²+x
P2(x) = x²-6x+8
	X
	-1
	0
	3
	5
	f(x)
	15
	8
	-1
	10
b)
	X
	Ordem 0
	Ordem 1
	Ordem 2
	Ordem 3
	-1
	15
	
	
	
	
	
	-7
	
	
	0
	8
	
	1
	
	
	
	-3
	
	7/60
	3
	-1
	
	17/10
	
	
	
	11/2
	
	
	5
	10
	
	
	
P3(x) = d0 + d1 (x-x0)+ d2 (x-x0).(x-x1)+ d3 (x-x0).(x-x1)(x-x2)
P3(x) =x²-6x+8 +
(x+1)(x-0)(x-3)
P3(x) = x²-6x+8+
(x³-2x²-3x) ( P3(x) = 
x³+
x²-
x+8
Exercício: Preencha a tabela de diferença divididas. Determine os polinômios interpoladores de graus 1,2 e 3. 
	X
	-1
	01
	2
	3
	F(x)
	1
	1
	0
	-1
	-2
	x
	Ordem 0
	Ordem 1
	Ordem 2
	Ordem 3
	Ordem 4
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
Exercícios
Encontre a forma de Lagrange para o polinômio interpolador, considerando a função dada pela tabela:
	x
	-1
	0
	2
	f(x)
	4
	1
	-1
Para um tanque de água, são fornecidos valores de temperatura em função da profundidade conforme a tabela a seguir:
	Profundidade (m), x
	Temperatura (oC), y
	1.0
	66
	1.5
	52
	2.0
	18
	2.5
	11
	3.0
	10
Calcule a temperatura estimada para 1,75 m usando interpolação linear.
Calcule a temperatura estimada para 2,75 m usando interpolação quadrática.
Calcule a temperatura estimada para 3,25 m usando interpolação cúbica.
Qual a relação entre o número de pontos usados na interpolação e o grau do polinômio interpolador que pode ser calculado?
Se você tiver um conjunto de 5 dados {(x0,f(x0), (x1,f(x1), (x2,f(x2), (x3,f(x3), (x4,f(x4),}, e deseja fazer uma interpolação linear, isto é, encontrar uma reta que lhe permita obter o valor de f(x(), onde x1<x(<x2:
	a) Qual seria o grau do polinômio que você calcularia, isto é, quantos pontos você utilizaria?
	b) E quais pontos da tabela você usaria?
A seguinte tabela informa o número de carros que passam por um determinado pedágio em um determinado dia: 
	Horário
	10:00
	10:30
	11:00
	11:30
	12:00
	12:30
	Número (em mil)
	2.69
	1.64
	1.09
	1.04
	1.49
	2.44
a) Faça um gráfico de horário vs. número de carros para verificar qual a tendência da curva.
b) Estime o número de carros que passariam pelo pedágio às 11:10, usando a forma de Lagrange para encontrar um polinômio interpolador p(x) que estima o número de carros em função do tempo. Use uma reta como função interpoladora.
c) Agora, faça a mesma estimativa, mas utilizando uma parábola como polinômio interpolador.
Na fabricação de determinadas cerâmicas é muito importante saber as condições de temperatura em que o produto foi assado no forno. Como não é possível medir a temperatura do forno a todo instante, ela é medida em intervalos periódicos de tempo e esses dados são interpolados para o instante em que cada peça foi “queimada” a fim de se conhecer a temperatura do forno nesse instante. Em um dia de funcionamento do forno, os seguintes dados foram coletados:
	Horário
	7:00
	10:00
	13:00
	16:00
	19:00
	21:00
	Temperatura (102 oC)
	2.32
	2.51
	2.63
	2.55
	2.41
	2.28
Construa a tabela de diferenças divididas para esses pontos.
Estime a temperatura do forno ás 14:30 usando a forma de Newton para apenas dois pontos.
Faça essa estimativa novamente, desta vez usando 3 pontos.
4. Integração numérica
Em determinadas situações, integrais são difíceis, ou mesmo impossíveis de se resolver analiticamente. Por ex
emplo: o valor de f(x) é conhecido apenas em alguns pontos, num intervalo [a, b]. Como não se conhece a expressão analítica de f(x), não é possível calcular a integral.
Ideia básica da integração numérica: substituição da função f(x) por um polinômio que a aproxime razoavelmente no intervalo [a, b]. 
Integração numérica de uma função f(x) num intervalo [a,b] cálculo da área delimitada por essa função, recorrendo à interpolação polinomial, como, forma de obtenção de um polinômio – pn(x). 
Faremos uso de métodos numéricos quando:
- f(x) for uma função muito difícil de integrar, contrariamente a um polinômio;
- a única informação sobre f(x) ser um conjunto de pares ordenados. 
Regra dos Trapézios Simples 
Seja f uma função real dada sobre dois pontos (xo, f(xo)) e (x1,f(x1)). Para esta função podemos construir um polinômio interpolador de Lagrange de grau menor ou igual a 1, que aproxima a função f(x). 
Seja h = x1-x0. E seja p1(x) = y0.L0(x)+ y1.L1(x) com
e
	
Logo, 
.
Regra dos Trapézios Composta (Repetida)
Se o intervalo [a, b] relativamente pequeno, a 
aproximação do valor do integral é aceitável. Mas, se
 o intervalo [a, b] é de grande amplitude, a 
aproximação será defasada. Assim, subdividimos em n sub-intervalos, e em cada um a função é aproximada por uma função linear.  A amplitude dos sub-intervalos será h=(b-a)/n . 
Pois, os termos f(x0) e f(xn) não se repetem.
Exemplo: Estimar o valor de 
.
	X
	0,0
	0,5
	1,0
	1,5
	2,0
	2,5
	3,0
	3,5
	4,0
	Y=(1+X²)-1/2
	1,00000
	0,89445
	0,70711
	0,55475
	0,44722
	0,37138
	0,31623
	0,27473
	0,24254
Regra dos Trapézios Simples - 2 pontos (x0 = 0,0 e x1 = 4,0) R: 2,48508
Regra dos Trapézios Composta - 3 pontos (x0=0,0, x1 =2,0, x2 = 4,0) R: 2,1369 
Regra dos Trapézios Composta - 9 pontos R: 2,0936 
A aproximação para 9 pontos é melhor, dado que o valor real é 2,0947.
Exemplo: Calcule 
onde f(x) é dada pela tabela: 
	x
	-3
	-1
	1
	3
	4
	5
	7
	f(x)
	4
	2
	0
	-1
	2
	4
	9
h=2
	x
	-3
	-1
	1
	3
	f(x)
	4
	2
	0
	-1
h=1
	x
	3
	4
	5
	f(x)
	-1
	2
	4
h=2		
	
	x
	5
	7
	f(x)
	4
	9
. Portanto,
= 7+3,5+13 = 23,5.
Regra 1/3 de Simpson
Novamente, podemos usar a fórmula de Lagrange para estabelecer a fórmula de integração resultante da aproximação de f(x) por um polinômio interpolador, agora de grau menor ou igual a 2. Para isto é preciso que f(x) seja dada sobre 3 pontos distintos: (x0, f(x0)); (x1,f(x1); (x2 f(x2)). Supondo que x1-x0 = x2 -x1=h, temos que x2-x0 = 2h. 
	
	
Exercício: Fazendo x-x0 = zh e substituindo x1 = x0+h, podemos calcular esta integral usando o método de mudança de variáveis. Chegando ao resultado:
Regra 1/3 de Simpson Repetida
	
	
Aplica-se a regra de Simpson repetidas vezes no intervalo [a, b] = [x0, xm] com h = xi - xi-1 (i = 1, 2,…, m ), para m par. Ou seja, x0, x1,..., xm são pontos igualmente espaçados e cada parábola precisa de três pontos consecutivos.
Obs: Sempre que possível utilize a regra de 1/3 de Simpson, pois chegamos a uma melhor estimativa do que a regra dos Trapézios.
Exemplos:
	x
	5
	6
	7
	f(x)
	0
	-1
	1
1) Dada a tabela calcule 
 pela regra de 1/3 de Simpson.
	x
	-3
	-1
	1
	3
	5
	6
	7
	8
	f(x)
	9
	7
	6
	2
	0
	-1
	1
	4
2) Dada a tabela calcule 
 pela regra possível.
 3) Calcule uma estimativa para a integral
 usando m = 2 e depois m =10 (ou seja, com 2 subintervalos: 3 pontos e 10 subintervalos: 11 pontos.
Exercícios Resolvidos
1) Calcule para a integral
.
Usando o Cálculo integral.
Usando a regra do Trapézio com 2 subintervalos: 3 pontos. Estime o erro.
	x
	0
	0,5
	1
	f(x)
	1
	1,2214
	1,4918
n=2 (subintervalos)( h = b-a/n = 1-0/2 = 1/2=0,5 
Para preencher a tabela use f(x) =e0,4x (que vem na integral)
( 
=
Usando a regra de 1/3 de Simpson com 10 subintervalos: 11 pontos. Estime o erro.
	x
	0
	0,1
	0,2
	0,3
	0,4
	0,5
	0,6
	0,7
	0,8
	0,9
	1
	F(x)
	1,0000
	1,0408
	1,0833
	1,1275
	1,1735
	1,2214
	1,2712
	1,3231
	1,3771
	1,4333
	1,4918
n= 10 subintervalos( h = (b-a)/n = (1-0)/10 =0,1
=1,2296
2) Calcule 
 pela melhor estimativa usando a tabela
	x
	0
	1
	2
	4
	5
	6
	7
	8
	f(x)
	3
	5
	7
	8
	4
	10
	11
	15
h=1
	X
	0
	1
	2
	f(x)
	3
	5
	7
=
h=2
	x
	2
	4
	f(x)
	7
	8
h=1			
	x
	4
	5
	6
	7
	8
	f(x)
	8
	4
	10
	11
	15
=103/3 = 34,33
Portanto, 
= 10+15+34,33 = 59,33
	x
	2
	4
	6
	f(x)
	8
	12
	20
3) Seja f(x) dada por 
Calcule 
 utilizando a regra dos trapézios.
Calcule 
 utilizando a regra de 1/3 de Simpson.
= 
(8+4.12+20) =
=50,6667
Calcule dois polinômiosde grau 1 que interpolam f(x) sobre os 3 pontos dados.
Escolhe os 2 primeiros pontos: x0 = 2 e x1 = 4
P(x) = y0L0+y1L1
P(x) = 8.(-0,5x+2)+12.(0,5x-1) = -4x+16+6x-12=2x+4
Escolhe outros 2 pontos: x0 = 4 e x1 = 6
P(x) = y0L0+y1L1
P(x) = 12.(-0,5x+3)+20.(0,5x-2) = -6x+36+10x-40=4x-4
Calcule o polinômio de grau 2 que interpola f(x) sobre os 3 pontos dados.
3 pontos( grau 2 Assim, x0 = 2, x1 = 4 e x2= 6
P2(x) = y0.L0+y1.L1+ y2.L2
P2(x) = 8. 
+12. 
+20. 
=x²-10x+24+3(-x²+8x-12)+2,5(x²-6x+8)
=x²-3x²+2,5x²-10x+24x-15x+24-36+20
=0,5x²-x+8
Calcule 
 utilizando os polinômios encontrados nos itens a e c.
=
=
= x²+4x
+2x²-4x
=
=(4²+4.4)-(2²+4.2) + (2.6²-4.6)-(2.4²-4.4) = 52 (igual ao item a) 
=
=
=
=36-18+48 -4/3 +2-16=152/3 =50,6667 (igual ao item b)
4) Dada a tabela abaixo 
	x
	1
	2
	3
	4
	5
	f(x)
	-1
	3
	23
	71
	159
Determine o polinômio interpolador linear e estime o valor f(1,5).
P(x) = y0L0+y1L1
P(x) = -1.(-x+2)+3.(x-1) = x-2+3x-3=4x-5
F(1,5) = p(1,5) =4.1,5-5 = 1
Determine o polinômio interpolador quadrático e estime o valor f(1,25).
grau 2(3 pontos: Assim, x0 = 1, x1 = 2 e x2= 3
P2(x) = y0.L0+y1.L1+ y2.L2
P2(x) = -1. 
+3. 
+23. 
=-0,5x²+2,5x-3-3x²+12x-9+11,5x²-34,5x+23
=8x²-20x+11 ( P(1,25)=-1,5
Determine o polinômio interpolador cúbico e estime o valor f(1,45).
P3(x) = y0.L0+y1.L1+ y2.L2
P3(x) = -1. (
)+3.(
)+23.(
)+71.(
)
=
P(x)=2x³-4x²+2x-1( P(1,45)=-0,41275
5. Sistemas de equações Lineares:
	Definição 1.2: Um sistema linear S de m equações lineares com n incógnitas 
 é um conjunto de m equações lineares, cada uma delas com n incógnitas, consideradas simultaneamente. Isto é:
	Uma solução de S é uma n-upla 
 de números reais que é solução de cada uma das equações do sistema.
Métodos Iterativos para a Solução de Sistemas Lineares
Seja os Sistema Linear 
 onde:
 matriz de coeficientes 
 vetor de variáveis 
 vetor independente (constantes) 
Ideia Geral dos Métodos Iterativos
Converter o sistema de equações 
 em um processo iterativo 
, onde:
 matriz com dimensões 
 vetor com dimensões 
 função de iteração matricial
Esquema Iterativo Proposto
Partindo de uma vetor aproximação inicial 
, constrói-se uma sequência iterativa de vetores:
Forma Geral: 
Observação
Se a sequência de aproximação 
, 
, 
, ......, 
 é tal que 
, então 
 é a solução do sistema 
.
 
Teste de Parada
Como em todos os processos iterativos, necessita-se de um critério para a parada do processo. 
Máximo desvio absoluto: 
Máximo desvio relativo: 
Desta forma, dada uma precisão 
, o vetor 
 será escolhido como solução aproximada da solução exata, se 
, ou dependendo da escolha, 
.
Método Iterativo de Gauss-Jacobi
Considere o sistema linear: 
Supondo 
, isola-se o vetor 
 mediante a separação pela diagonal da matriz de coeficientes.
Assim, tem-se o sistema iterativo 
, onde:
 
Dado uma aproximação inicial 
, o Método de Gauss-Jacobi consiste em obter uma sequência 
, 
, 
, ......, 
, por meio da relação recursiva: 
 
Observe que o processo iterativo utiliza somente estimativas da iteração anterior.
Exemplo: Resolver o sistema de equações lineares, pelo Método de Gauss-Jacobi com solução inicial 
 e tolerância 
.
Separando-se os elementos diagonais, tem-se:
Solução para k=0 
Cálculo de 
:
Para k=1:
Para k=2:
 é solução com erro menor que 0,05.
Condições Suficientes para a Convergência do Método de Gauss-Jacobi
Teorema: Seja o sistema linear
 e seja 
Se 
, então o método G-J gera uma sequência 
 convergente para a solução do sistema dado, independentemente da escolha da aproximação inicial 
.
Observe que esta é uma condição suficiente, se for satisfeita o método converge, entretanto se não for satisfeita nada se pode afirmar.
Exemplo 1: Seja a matriz do exemplo dado anteriormente:
Tem-se a convergência garantida para qualquer vetor inicial.
Exemplo 2: Seja o sistema de equações lineares:
 
 
As condições de convergência do teorema não são satisfeitas, entretanto o Método de Gauss-Jacobi gera uma sequência convergente para a solução exata 
. Se as condições de suficiência não são satisfeitas, não significa que o método não possa convergir.
Exemplo 3: Considere o sistema linear:
As condições do teorema não são satisfeitas. Uma solução possível é permutar as equações. Seja no exemplo permutar a primeira equação com a segunda equação:
As condições passam a ser satisfeitas e a convergência é garantida para qualquer vetor inicial. Este tipo de procedimento nem sempre é possível.
Fórmula Matricial do Método Gauss-Jacobi 
Decompõe-se a matriz de coeficientes 
 em:
Onde:
L – Matriz Triangular Inferior
D – Matriz Diagonal
U – Matriz Triangular Superior
Método Iterativo de Gauss-Seidel
Assim como no Método de Gauss-Jacobi o sistema linear 
 é escrito na forma equivalente: 
Como no Método Gauss-Jacobi, é realizada uma separação diagonal, e o processo iterativo de atualização é seqüencial, componente por componente. A diferença é que, no momento de realizar-se a atualização das componentes do vetor numa determinada iteração, a formulação utiliza as componentes da iteração já atualizadas na iteração atual, com as restantes não atualizadas da iteração anterior. Por exemplo, ao se calcular a componente 
 da iteração (k+1), utiliza-se no cálculo as componentes já atualizadas 
 com as componentes ainda não atualizadas da iteração anterior 
.
Exemplo: Resolver o sistema linear utilizando o Método Iterativo de Gauss-Seidel, com 
 e tolerância 
.
O processo iterativo é dado por:
Para k=0 e 
 ( 
Cálculo de 
:
Para k=1 e 
:
Para k=2 e 
:
 é solução com erro menor que 0,05.
Fórmula Matricial do Método Gauss-Seidel
Decompõe-se a matriz de coeficientes 
 em:
. Onde: L – Matriz Triangular Inferior, D – Matriz Diagonal, U – Matriz Triangular Superior
�� EMBED Equation.3 
Interpretação Geométrica do Caso 
Considere o Sistema Linear:
O esquema iterativo utilizando o Método de Gauss-Seidel é dado por:
Para k=0 e 
:
Para k=1 e 
:
Para k=2 e 
:
A solução exata é dada por: 
.
Esse processo iterativo até à convergência pode ser interpretado geometricamente num gráfico com a componente 
 na abscissa e a componente 
na ordenada.
Observação 1: Verifica-se pelo gráfico acima que a sequência 
, 
, 
, ......, 
 está convergindo para a solução exata 
.
Observação 2: A sequência gerada pelo Método de Gauss-Seidel depende fortemente da posição das equações. Esta observação pode ser melhor entendida modificando a ordem das equações do exemplo anterior.
Considere o mesmo Sistema Linear anterior, porém modificando a ordem das equações:
O esquema iterativo utilizando o Método de Gauss-Seidel é dado por:
Para k=0 e 
:
Para k=1 e 
:
Estudo da Convergência do Método de Gauss-Seidel
Existem dois critérios de suficiência para a convergência do Método de Gauss-Seidel. O critério de linhas e o critério de Sassenfeld. O critério de linhas é o mesmo da Método de Gauss-Jacobi.
Critério de Linhas:Seja o sistema linear
, com A dimensão 
 e seja: 
Se 
, então o método Gauss-Seidel gera uma seqüência 
 convergente para a solução do sistema dado, independentemente da escolha da aproximação inicial 
. A matriz que satisfizer o critério de linhas é chamada de diagonal dominante estrita.Critério de Sassenfeld
Seja o sistema linear
, com A dimensão 
 e seja:
 e para 
:
Define-se 
.
Se 
, então o Método de Gauss-Seidel gera uma sequência convergente para a solução do sistema, qualquer que seja o vetor inicial. Além disso, quanto menor for o valor de 
 mais rápida é a convergência.
Exemplo: Verificar as condições de convergência do Método de Gauss-Seidel no sistema abaixo:
Critério de Linhas 
 não satisfaz.
Critério de Sassenfeld 
 não satisfaz.
Como a convergência do Método de Gauss-Seidel é fortemente dependente da posição das equações, pode-se trocar a posição das equações.
Tentativa 1: Troca-se a primeira equação pela terceira equação.
Critério de Linhas 
 não satisfaz.
Critério de Sassenfeld 
 não satisfaz.
Tentativa 2: Troca-se a primeira coluna pela terceira coluna na equação anterior.
Critério de Linhas
 satisfaz. 
 não satisfaz.
Critério de Sassenfeld
 satisfaz. 
 satisfaz. 
 satisfaz.
Com a última modificação o sistema passa a ser convergente para qualquer vetor inicial. Modificações desse tipo são puramente acadêmicas e são difíceis de serem realizadas em sistemas reais. Principalmente pelas dimensões dos problemas, resultando num grande esforço computacional, e das incertezas quanto a sua eficiência. 
Exemplo : Verifique a convergência do sistema abaixo pelo critério de linhas e Sassenfeld.
Critério de Linhas 
 Não satisfaz
Critério de Sassenfeld 
 Satisfaz
Exemplo : Verifique a convergência do sistema abaixo pelo critério de linhas e Sassenfeld.
Lista de exercícios
1) Dada aproximação inicial x(0) = (0,9 1,9 1,9) determine a solução do sistema de equações lineares, 
, usando ( = 0,01. Use 4 casas decimais. 
a) Pelo método de Gauss- Jacobi.
b) Pelo método de Gauss- Seidel.
2) Dada aproximação inicial x(0) = (0,7 0,7 0,7 0,7) determine a solução do sistema de equações lineares, 
, usando ( = 0,1. Use 4 casas decimais. Pelo método de Gauss- Seidel.
3) Dada aproximação inicial x(0) = (0,7 -0,5 0,6 -0,8) determine a solução do sistema de equações lineares, 
, usando ( = 0,05 ou 3 iterações (o que ocorrer primeiro). 
a) Pelo método de Gauss- Jacobi.
b) Pelo método de Gauss- Seidel.
4) Dada aproximação inicial x(0) = (0 0 0) determine a solução do sistema de equações lineares, 
, usando ( = 0,05 ou 3 iterações (o que ocorrer primeiro). Sendo A B C D os 4 últimos dígitos do seu RA. (Um dígito para cada letra).
a) Pelo método de Gauss- Jacobi.
b) Pelo método de Gauss- Seidel.
� EMBED Equation.3 ���
� EMBED Equation.3 ���
� EMBED Equation.3 ���
� EMBED Equation.3 ���
� EMBED Equation.3 ���
� EMBED Equation.3 ���
�PAGE �
�PAGE �42�
_1328979200.unknown
_1477199715.unknown
_1533109807.unknown
_1533109823.unknown
_1533109839.unknown
_1533109848.unknown
_1533109852.unknown
_1533109856.unknown
_1533110095.unknown
_1533110331.unknown
_1533126401.unknown
_1533110222.unknown
_1533110094.unknown
_1533109854.unknown
_1533109855.unknown
_1533109853.unknown
_1533109850.unknown
_1533109851.unknown
_1533109849.unknown
_1533109843.unknown
_1533109845.unknown
_1533109847.unknown
_1533109844.unknown
_1533109841.unknown
_1533109842.unknown
_1533109840.unknown
_1533109831.unknown
_1533109835.unknown
_1533109837.unknown
_1533109838.unknown
_1533109836.unknown
_1533109833.unknown
_1533109834.unknown
_1533109832.unknown
_1533109827.unknown
_1533109829.unknown
_1533109830.unknown
_1533109828.unknown
_1533109825.unknown
_1533109826.unknown
_1533109824.unknown
_1533109815.unknown
_1533109819.unknown
_1533109821.unknown
_1533109822.unknown
_1533109820.unknown
_1533109817.unknown
_1533109818.unknown
_1533109816.unknown
_1533109811.unknown
_1533109813.unknown
_1533109814.unknown
_1533109812.unknown
_1533109809.unknown
_1533109810.unknown
_1533109808.unknown
_1533109393.unknown
_1533109799.unknown
_1533109803.unknown
_1533109805.unknown
_1533109806.unknown
_1533109804.unknown
_1533109801.unknown
_1533109802.unknown
_1533109800.unknown
_1533109397.unknown
_1533109399.unknown
_1533109400.unknown
_1533109398.unknown
_1533109395.unknown
_1533109396.unknown
_1533109394.unknown
_1533109385.unknown
_1533109389.unknown
_1533109391.unknown
_1533109392.unknown
_1533109390.unknown
_1533109387.unknown
_1533109388.unknown
_1533109386.unknown
_1533107624.unknown
_1533108387.unknown
_1533108388.unknown
_1533108384.unknown
_1477199719.unknown
_1477199721.unknown
_1477199722.unknown
_1477199723.unknown
_1477199720.unknown
_1477199717.unknown
_1477199718.unknown
_1477199716.unknown
_1329580496.unknown
_1477199699.unknown
_1477199707.unknown
_1477199711.unknown
_1477199713.unknown
_1477199714.unknown
_1477199712.unknown
_1477199709.unknown
_1477199710.unknown
_1477199708.unknown
_1477199703.unknown
_1477199705.unknown
_1477199706.unknown
_1477199704.unknown
_1477199701.unknown
_1477199702.unknown
_1477199700.unknown
_1329581143.unknown
_1329581376.unknown
_1329583873.unknown
_1329583927.unknown
_1453462450.unknown
_1329581555.unknown
_1329581175.unknown
_1329581191.unknown
_1329581172.unknown
_1329581047.unknown
_1329581089.unknown
_1329581098.unknown
_1329581061.unknown
_1329580968.unknown
_1329581035.unknown
_1329580951.unknown
_1328980307.unknown
_1328981171.unknown
_1329058047.unknown
_1329058147.unknown
_1329058168.unknown
_1329058070.unknown
_1328981220.unknown
_1329057954.unknown
_1328981192.unknown
_1328980835.unknown
_1328981067.unknown
_1328980973.unknown
_1328981016.unknown
_1328980796.unknown
_1328980805.unknown
_1328980784.unknown
_1328980773.unknown
_1328979887.unknown
_1328979995.unknown
_1328980109.unknown
_1328980283.unknown
_1328980044.unknown
_1328979930.unknown
_1328979979.unknown
_1328979900.unknown
_1328979553.unknown
_1328979636.unknown
_1328979808.unknown
_1328979579.unknown
_1328979514.unknown
_1328979525.unknown
_1328979506.unknown
_1328376277.unknown
_1328964874.unknown
_1328967828.unknown
_1328968465.unknown
_1328978642.unknown
_1328979045.unknown
_1328979155.unknown
_1328978990.unknown
_1328978249.unknown
_1328978260.unknown
_1328978146.unknown
_1328978175.unknown
_1328978100.unknown
_1328967972.unknown
_1328968040.unknown
_1328968083.unknown
_1328968019.unknown
_1328967919.unknown
_1328967938.unknown
_1328967909.unknown
_1328967407.unknown
_1328967594.unknown
_1328967615.unknown
_1328967645.unknown
_1328967505.unknown
_1328967532.unknown
_1328967458.unknown
_1328967427.unknown
_1328966546.unknown
_1328966736.unknown
_1328967004.unknown
_1328967095.unknown
_1328966725.unknown
_1328966422.unknown
_1328966443.unknown
_1328965110.unknown
_1328963554.unknown
_1328964117.unknown
_1328964318.unknown
_1328964363.unknown
_1328964727.unknown
_1328964354.unknown
_1328964242.unknown
_1328964274.unknown
_1328964150.unknown
_1328963889.unknown
_1328964081.unknown
_1328963609.unknown
_1328963838.unknown
_1328963877.unknown
_1328963632.unknown
_1328963583.unknown
_1328377637.unknown
_1328963460.unknown
_1328963518.unknown
_1328963538.unknown
_1328963499.unknown
_1328963392.unknown
_1328963416.unknown
_1328446009.unknown
_1328963340.unknown
_1328446010.unknown
_1328446008.unknown
_1328377163.unknown
_1328377236.unknown
_1328377350.unknown
_1328377205.unknown
_1328376986.unknown
_1328377131.unknown_1328376379.unknown
_1280657324.unknown
_1328373735.unknown
_1328374944.unknown
_1328375263.unknown
_1328375692.unknown
_1328375916.unknown
_1328375961.unknown
_1328375725.unknown
_1328375495.unknown
_1328375663.unknown
_1328375302.unknown
_1328375232.unknown
_1328375234.unknown
_1328375235.unknown
_1328375233.unknown
_1328375046.unknown
_1328375077.unknown
_1328374945.unknown
_1328374626.unknown
_1328374715.unknown
_1328374783.unknown
_1328374810.unknown
_1328374738.unknown
_1328374677.unknown
_1328374701.unknown
_1328374658.unknown
_1328374172.unknown
_1328374453.unknown
_1328374491.unknown
_1328374348.unknown
_1328374062.unknown
_1328374167.unknown
_1328374035.unknown
_1328371646.unknown
_1328372473.unknown
_1328373314.unknown
_1328373504.unknown
_1328373710.unknown
_1328373422.unknown
_1328372544.unknown
_1328373309.unknown
_1328372505.unknown
_1328372173.unknown
_1328372325.unknown
_1328372188.unknown
_1328372298.unknown
_1328371687.unknown
_1328372111.unknown
_1328371663.unknown
_1298310605.unknown
_1298482395.unknown
_1298733365.unknown
_1328370885.unknown
_1328371061.unknown
_1328371474.unknown
_1328371520.unknown
_1328371426.unknown
_1328370767.unknown
_1298482709.unknown
_1298482876.unknown
_1298482915.unknown
_1298483033.unknown
_1298482844.unknown
_1298482564.unknown
_1298482670.unknown
_1298482549.unknown
_1298455589.unknown
_1298482250.unknown
_1298482289.unknown
_1298482237.unknown
_1298455545.unknown
_1298455561.unknown
_1298455527.unknown
_1280662919.unknown
_1282487326.unknown
_1285506694.unknown
_1285508833.unknown
_1285510773.unknown
_1286104116.unknown
_1286104557.unknown
_1286104924.unknown
_1286104979.unknown
_1286104923.unknown
_1286104269.unknown
_1286101957.unknown
_1286103227.unknown
_1286102722.unknown
_1286101909.unknown
_1285510631.unknown
_1285510729.unknown
_1285509784.unknown
_1285510386.unknown
_1285509202.unknown
_1285508748.unknown
_1282487972.unknown
_1285506606.unknown
_1285505707.unknown
_1282487940.unknown
_1282482057.unknown
_1282487080.unknown
_1282481423.unknown
_1280657523.unknown
_1280658676.unknown
_1280658799.unknown
_1280658652.unknown
_1280657445.unknown
_1280657518.unknown
_1280657350.unknown
_1280657444.unknown
_1271931084.unknown
_1279459791.unknown
_1280648511.unknown
_1280657214.unknown
_1280657270.unknown
_1280657304.unknown
_1280657286.unknown
_1280657242.unknown
_1280657258.unknown
_1280648582.unknown
_1280656333.unknown
_1280648533.unknown
_1280066934.unknown
_1280648442.unknown
_1280648457.unknown
_1280648433.unknown
_1279459854.unknown
_1279460007.unknown
_1279459799.unknown
_1279459184.unknown
_1279459726.unknown
_1279459735.unknown
_1279459679.unknown
_1279459707.unknown
_1279459246.unknown
_1271931682.unknown
_1271931747.unknown
_1271932398.unknown
_1271934402.unknown
_1279458924.unknown
_1271934173.unknown
_1271932052.unknown
_1271931702.unknown
_1271931163.unknown
_1271931199.unknown
_1271931114.unknown
_1111124697.unknown
_1111130367.unknown
_1271930714.unknown
_1271930938.unknown
_1271931004.unknown
_1271931018.unknown
_1271930953.unknown
_1271930806.unknown
_1271930929.unknown
_1271930789.unknown
_1111405036.unknown
_1111414206.unknown
_1111416020.unknown
_1113134326.unknown
_1113134810.unknown
_1113134823.unknown
_1113134789.unknown
_1111416304.unknown
_1111817501.unknown
_1111817704.unknown
_1113134311.unknown
_1111817587.unknown
_1111416404.unknown
_1111416907.unknown
_1111416226.unknown
_1111416259.unknown
_1111416178.unknown
_1111414921.unknown
_1111415938.unknown
_1111415964.unknown
_1111415161.unknown
_1111415177.unknown
_1111414493.unknown
_1111414824.unknown
_1111414315.unknown
_1111410502.unknown
_1111413643.unknown
_1111413766.unknown
_1111414165.unknown
_1111413662.unknown
_1111412473.unknown
_1111413609.unknown
_1111412470.unknown
_1111412471.unknown
_1111410710.unknown
_1111408429.unknown
_1111408541.unknown
_1111408557.unknown
_1111407415.unknown
_1111408144.unknown
_1111408401.unknown
_1111408409.unknown
_1111408383.unknown
_1111408021.unknown
_1111405666.unknown
_1111222282.unknown
_1111403213.unknown
_1111404521.unknown
_1111405001.unknown
_1111405027.unknown
_1111404578.unknown
_1111403281.unknown
_1111403772.unknown
_1111404040.unknown
_1111402133.unknown
_1111402988.unknown
_1111402514.unknown
_1111402619.unknown
_1111402213.unknown
_1111401543.unknown
_1111401604.unknown
_1111131736.unknown
_1111132265.unknown
_1111221188.unknown
_1111221198.unknown
_1111221999.unknown
_1111132550.unknown
_1111131900.unknown
_1111132173.unknown
_1111131668.unknown
_1111131053.unknown
_1111131096.unknown
_1111131179.unknown
_1111131253.unknown
_1111130558.unknown
_1111130629.unknown
_1111130781.unknown
_1111130441.unknown
_1111127379.unknown
_1111128634.unknown
_1111129832.unknown
_1111129911.unknown
_1111130145.unknown
_1111130046.unknown
_1111129885.unknown
_1111129153.unknown
_1111129269.unknown
_1111128789.unknown
_1111128948.unknown
_1111127668.unknown
_1111128068.unknown
_1111128455.unknown
_1111128022.unknown
_1111127564.unknown
_1111127623.unknown
_1111127415.unknown
_1111125844.unknown
_1111126367.unknown
_1111126746.unknown
_1111126786.unknown
_1111126558.unknown
_1111125892.unknown
_1111125281.unknown
_1111125601.unknown
_1111125615.unknown
_1111125633.unknown
_1111125371.unknown
_1111125351.unknown
_1111125210.unknown
_1111125244.unknown
_1111125136.unknown
_1111124728.unknown
_1111125109.unknown
_1111124343.unknown
_1111124380.unknown
_1111124516.unknown
_1111124630.unknown
_1111124388.unknown
_1111124372.unknown
_1111124252.unknown
_1111124304.unknown
_1111124328.unknown
_1111123997.unknown