Enunciado

Em cada um dos problemas de 1 a 6, determine um valor aproximado da solução em t = 0,4 e t = 0,5 usando o método especificado. Para os valores iniciais, use o método de Runge-Kutta; veja os Problemas de 1 a 6 da Seção 8.3. Compare os resultados dos vários métodos entre si e com a solução exata (se disponível).

(a) Use o método de previsão ecorreção de quarta ordem comh=0,1. Use a fórmula de correção uma vez em cada passo.

(b) Use o método de Adam-Moulton de quarta ordem com h = 0,1.

(c) Use o método inverso de diferenciação de quarta ordem com h = 0,1.

1. y ' = 5 t - 3 y   ,   y 0 = 2

Passo 1

a) Vamos começar usando a notação f n = f ( t n , y n ) .

Desse modo, a fórmula do método preditor vai ficar assim...

y n + 1 = y n + h 24 55 f n - 59 f n - 1 + 37 f n - 2 - 9 f n - 3

Já no caso da fórmula de correção com f n + 1 = f ( t n + 1 , y n + 1 ) ... Temos que:

y n + 1 = y n + h 24 ( 9 f n + 1 + 19 f n - 5 f n - 1 + f n - 2 )

Usando os valores iniciais gerados pelo método de R u n g e - K u t t a . No caso peguemos os valores da figura abaixo do código que fizemos na S e ç ã o   3 para essa mesma E D O com passo h = 0.1 ;

Ou seja, pegando a seguinte tabela abaixo...

Para o método preditor com t n = 0.40 ...

y n + 1 = y n + h 24 55 f n - 59 f n - 1 + 37 f n - 2 - 9 f n - 3

Sendo que...

f n = 5 ⋅ 0.3 - 3 1.12686 = - 1.68461

f n - 1 = 5 ⋅ 0.2 - 3 1.33362 = - 2.46447

f n - 2 = 5 ⋅ 0.1 - 3 1.62231 = - 3.321098

f n - 3 = 5 ⋅ 0 - 3 2 = - 4.24264

Obtemos que...

y n + 1 = 1.12686 + 0.1 24 [ 55 ⋅ ( - 1.6846 ) - 59 ⋅ ( - 2.4644 ) + 37 ⋅ ( - 3.3211 ) - 9 ⋅ ( - 4.246 ) ]

y n + 1 = 0.993751

Passo 2

No caso do método corretor, temos que para t n = 0.4 ...

y n + 1 = y n + h 24 ( 9 f n + 1 + 19 f n - 5 f n - 1 + f n - 2 )

Com...

f n + 1 = 5 ⋅ 0.4 - 3 0.993751 = - 0.990612

f n = 5 ⋅ 0.3 - 3 1.12686 = - 1.68461

f n - 1 = 5 ⋅ 0.2 - 3 1.33362 = - 2.46447

f n - 2 = 5 ⋅ 0.1 - 3 1.62231 = - 3.321098

Chegamos ao seguinte resultado:

y n + 1 = 1.12686 + 0.1 24 [ 9 ⋅ - 0.9906 + 19 ⋅ - 1.68461 - 5 ⋅ - 2.46447 - 3.321098 ]

y n + 1 = 0.993852

Seguindo o mesmo procedimento para outros t n = 0.5 , com n = 5 .... Chegamos até a seguinte tabela abaixo...

Passo 3

b) Presta atenção! O método de Adam-Moulton tem a seguinte cara abaixo...

Com f n + 1 = f ( t n + 1 , y n + 1 ) , temos que...

y n + 1 = y n + h 24 ( 9 f n + 1 + 19 f n - 5 f n - 1 + f n - 2 )

Até aí tá bem parecido com o método corretor. O diferencial é que...

f n + 1 = 5 t n + 1 - 3 y n + 1

Uma vez que a E D O é linear, nós vamos fazer...

y n + 1 = y n + h 24 45 t n + 1 - 27 y n + 1 + 19 f n - 5 f n - 1 + f n - 2

E com isso, para t n = 0.4 ... Sabendo dos valores abaixo...

y n = 1.12686

y n + 1 = 0.993852

h = 0.1

t n + 1 = 0.5

f n = 5 ⋅ 0.3 - 3 1.12686 = - 1.68461

f n - 1 = 5 ⋅ 0.2 - 3 1.33362 = - 2.46447

f n - 2 = 5 ⋅ 0.1 - 3 1.62231 = - 3.321098

Alcançamos o resultado para y ( 0.4 ) :

y n + 1 = y n + h 24 45 t n + 1 - 27 y n + 1 + 19 f n - 5 f n - 1 + f n - 2

y n + 1 = 1.1268 + 0.1 24 45 ⋅ 0.5 - 27 0.9938 + 19 ⋅ - 1.6846 - 5 ⋅ - 2.4644 - - 3.3211

y n + 1 = 0.925746

Seguindo o mesmo procedimento para t n = 0.5 ...Obtemos a seguinte tabela abaixo

Passo 4

c) No caso do método inverso de diferenciação de quarta ordem com h   =   0.1 ... Bora fazer;

y n + 1 = 1 25 48 y n - 36 y n - 1 + 16 y n - 2 - 3 y n - 3 + 12 h f n + 1

Só que neste método, o diferencial assim como no método anterior, é que...

f n + 1 = 5 t n + 1 - 3 y n + 1

O que nos faz ficar com...

y n + 1 = 1 25 48 y n - 36 y n - 1 + 16 y n - 2 - 3 y n - 3 + 12 h ( 5 t n + 1 - 3 y n + 1 )

E como já sabemos todos os valores de cada variável aí, fazendo as iterações, para t n = 0.4 e t n = 0.5 ... Obtemos a tabela abaixo:

Resposta

a)

b)

c)

Ver Outros Exercícios desse livro

Exercícios de Livros Relacionados

Em cada um dos problemas de 13 a 16, encontre a solução do problema de valor inicial dado. Faça o gráfico da solução e descreva como ela se comporta quando x → 0 16. x 2 y ″ +3 x y ' +5 y=0 y ( 1 )=1
Em cada um dos problemas de 17 a 34, encontre todos os pontos singulares da equação dada e determine se cada um deles é regular ou irregular.21. 1–x 2 2 y ″+x ( 1–x ) y '+( 1+x ) y=0
Em cada um dos problemas de 13 a 16, encontre a solução do problema de valor inicial dado. Faça o gráfico da solução e descreva como ela se comporta quando x → 0 14. 4 x 2 y ″ +8 x y ' +17 y=0 y ( 1 )
Em cada um dos problemas de 1 a 12, determine a solução geral da equação diferencial dada, válida em qualquer intervalo que não inclui o ponto singular1. x 2 y ″+4 x y '+2 y=0
Em cada um dos problemas de 1 a 12, determine a solução geral da equação diferencial dada, válida em qualquer intervalo que não inclui o ponto singular6. ( x - 1 ) 2 y ″ + 8 ( x - 1 ) y ' + 12 y=0