Calculo Numerico
Calculo Numerico
Calculo Numerico
yk +1 = yk + hφ(xk , yk ), ∀k = 0, 1, . . . ,
yk +1 = yk + hf(xk , yk ), ∀k = 0, 1, . . . ,
Exemplo 1
Use o método de Euler para obter uma aproximação numérica da
solução do sistema de equações diferenciais
(
y 0 = z, y (0) = 1
0
z =y +e x z(0) = 0
Para k = 0, temos
y1 y0 z0 1 0 1
= +h x = + 0.1 = .
z1 z0 y0 + e 0 0 1+1 0.2
Para k = 1, temos
y2 y1 z1 1 0.2 1.02
= +h = + 0.1 = .
z2 z1 y1 + ex1 0.2 1 + e0.2 0.4105
Exemplo 2
Use o método de Heun para obter uma aproximação numérica da
solução do sistema de equações diferenciais
(
y 0 = z, y (0) = 1
0
z =y +e x z(0) = 0
Logo,
y1 y h 1 0.1 0 0.2
= 0 + (k1 + k2 ) = + +
z1 z0 2 0 2 2 2.1052
1.01000
=
0.20526
h
yk +1 = yk + (k1 + 2k2 + 2k3 + k4 ) ,
6
em que
k1 = f(xk , yk ),
k2 = f(xk + h/2, yk + k1 h/2),
k3 = f(xk + h/2, yk + k2 h/2),
k4 = f(xk + h, yk + k3 h).
Logo,
h 1.0102
y1 = y0 + (k1 + 2k2 + 2k3 + k4 ) = .
6 0.2055
Logo,
h 1.0415
y2 = y1 + (k1 + 2k2 + 2k3 + k4 ) = .
6 0.4241
e a condição inicial
p0 = 80 e q0 = 30.
Equivalentemente, temos
(
p0 = 0.25p − 0.01pq, p(t0 ) = 80,
q 0 = −1.00q + 0.01pq, q(t0 ) = 30.
140
h=1.00
h=0.50
130 h=0.25
Exata
120
110
Presa
100
90
80
70
60
10 15 20 25 30 35 40 45
Marcos Eduardo Valle MS211 - Cálculo
PredadorNumérico 15 / 20
Método de Heun
140
h=1.00
h=0.50
130 h=0.25
Exata
120
110
Presa
100
90
80
70
60
10 15 20 25 30 35 40 45
Marcos Eduardo Valle MS211 - Cálculo
PredadorNumérico 16 / 20
Método de Runge-Kutta ordem 4
140
h=1.00
h=0.50
130 h=0.25
Exata
120
110
Presa
100
90
80
70
60
10 15 20 25 30 35 40 45
Marcos Eduardo Valle MS211 - Cálculo
PredadorNumérico 17 / 20
Equações Diferenciais de Ordem Superior
Uma equação diferencial de ordem m
u (m) = g(x, u, u 0 , u 00 , . . . , u (m−1) ),
pode ser escrita como um sistema com m equações diferenciais de
ordem 1:
y0 = f (x, y),
tomando
y1 (x) = u(x), y2 (x) = u 0 (x), ..., ym (x) = u (m−1) .
Especificamente, temos
y10 = y2 ,
y 0
= y3 ,
2
..
.
0
ym = g(x, y1 , y2 , . . . , ym ).
Marcos Eduardo Valle MS211 - Cálculo Numérico 18 / 20
Exemplo 4
Escreva o problema de valor inicial
(
u 00 − u = ex ,
u(0) = 1 e u 0 (0) = 0,