`ode23tb`

Solver TR-BDF2 para equações diferenciais rígidas em tolerâncias grosseiras.

Categoria: Mathematics > Numerical Integration and Differential Equations

Sintaxe

[t, y] = ode23tb(odefun, tspan, y0)

Integra y' = odefun(t, y).

[t, y] = ode23tb(odefun, tspan, y0, options)

Aplica também opções de odeset.

Descrição

Cada passo combina um estágio de regra do trapézio com um estágio de diferenciação para trás de segunda ordem, usando a mesma matriz de iteração para ambos. O método é implícito e L-estável, e portanto amortece bem as componentes rígidas.

Como ode23s, pode superar ode15s quando as tolerâncias são folgadas.

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

Decaimento rígido

f = @(t, y) -100*(y - cos(t));
[t, y] = ode23tb(f, [0 1], 0);
fprintf('%.3f\n', y(end))

Saída

0.549

Veja também

ode15s ode23t ode23s odeset