Campus Livre
ode23
Solver de Runge-Kutta de baixa ordem para equações diferenciais não rígidas.
Categoria: Mathematics > Numerical Integration and Differential Equations
Sintaxe
[t, y] = ode23(odefun, tspan, y0)Integra y' = odefun(t, y) em tspan a partir de y0.
[t, y] = ode23(odefun, tspan, y0, options)Aplica também opções de odeset.
sol = ode23(___)Retorna uma estrutura para uso com deval.
Descrição
ode23 combina uma fórmula de Runge-Kutta de segunda e outra de terceira ordem para estimar o erro e adaptar o passo. É mais barato por passo que ode45, mas precisa de mais passos para uma dada precisão, sendo a melhor escolha para tolerâncias grosseiras ou problemas levemente rígidos.
Quando tspan tem apenas duas entradas, o solver retorna os pontos que escolheu; com mais entradas, retorna a solução exatamente nesses instantes.
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.
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] = ode23(@(t, y) -y, [0 1], 1);
fprintf('%.4f %.4f\n', y(end), exp(-1))Saída
0.3679 0.3679