Campus Livre
ode89
Solver de Runge-Kutta de oitava e nona ordem para equações não rígidas com tolerâncias muito apertadas.
Categoria: Mathematics > Numerical Integration and Differential Equations
Sintaxe
[t, y] = ode89(odefun, tspan, y0)Integra y' = odefun(t, y).
[t, y] = ode89(odefun, tspan, y0, options)Aplica também opções de odeset.
Descrição
ode89 é o solver explícito de maior ordem da família. Compensa em soluções muito suaves integradas por longos intervalos com tolerâncias rigorosas, em que atinge uma dada precisão com o menor número de avaliações da função.
Como todo método explícito, dá passos minúsculos em problemas rígidos; use ode15s nesses casos.
Argumentos de entrada
odefun- Handle
dydt = odefun(t, y). tspan[t0 tf]ou um vetor de instantes de saída.y0- Vetor de valores iniciais.
Argumentos de saída
t, y- Instantes e solução, com o mesmo formato de
ode45.
Exemplo
Tolerância rigorosa
opts = odeset('RelTol', 1e-12, 'AbsTol', 1e-14);
[t, y] = ode89(@(t, y) -y, [0 1], 1, opts);
fprintf('%.10f %.10f\n', y(end), exp(-1))Saída
0.3678794412 0.3678794412