Campus Livre
ode113
Solver de Adams-Bashforth-Moulton de ordem variável para equações diferenciais não rígidas.
Categoria: Mathematics > Numerical Integration and Differential Equations
Sintaxe
[t, y] = ode113(odefun, tspan, y0)Integra y' = odefun(t, y) em tspan a partir da condição inicial y0.
[t, y] = ode113(odefun, tspan, y0, options)Aplica também opções de odeset.
sol = ode113(___)Retorna uma estrutura para uso com deval.
Descrição
Um método de múltiplos passos reaproveita valores anteriores da derivada, de modo que cada passo precisa de poucas avaliações de odefun. Isso torna ode113 eficiente quando a função é cara de avaliar e a solução é suave, e ela lida bem com tolerâncias rigorosas.
Não serve para problemas rígidos (use ode15s) nem para soluções com descontinuidades abruptas, em que um método de passo único como ode45 é melhor.
Argumentos de entrada
odefun- Handle
dydt = odefun(t, y)que retorna um vetor coluna. tspan[t0 tf], ou um vetor mais longo de instantes em que reportar a solução.y0- Vetor de valores iniciais.
options- Estrutura de
odeset(tolerâncias, eventos etc.).
Argumentos de saída
t- Vetor coluna de instantes.
y- Array com uma linha por instante e uma coluna por componente.
Exemplo
Decaimento exponencial
[t, y] = ode113(@(t, y) -y, [0 1], 1);
fprintf('%.4f %.4f\n', y(end), exp(-1))Saída
0.3679 0.3679