`bvp5c`

Resolve problemas de valor de contorno com um método de quarta ordem que controla o erro verdadeiro.

Categoria: Mathematics > Numerical Integration and Differential Equations

Sintaxe

sol = bvp5c(odefun, bcfun, solinit)

Resolve o sistema y' = odefun(x, y) com condições de contorno bcfun(y(a), y(b)) = 0.

sol = bvp5c(odefun, bcfun, solinit, options)

Aplica também opções de bvpset.

Descrição

bvp5c tem o mesmo padrão de chamada de bvp4c, então muitas vezes uma pode ser trocada pela outra. A diferença está no controle do erro: bvp4c limita o resíduo do polinômio de colocação, enquanto bvp5c limita o erro na própria solução, o que pode importar quando se exige tolerância apertada.

Se um solver trava em um problema difícil, experimente o outro e melhore o palpite inicial.

Argumentos de entrada

odefun, bcfun: Handles de função para as equações e para os resíduos de contorno, como em bvp4c.
solinit: Malha inicial e palpite de bvpinit.
options: Estrutura de bvpset.

Argumentos de saída

sol: Estrutura com os campos x, y, yp e stats, utilizável com deval.

Exemplo

Mesmo problema de bvp4c

odefun = @(x, y) [y(2); -y(1)];
bcfun = @(ya, yb) [ya(1); yb(1) - 1];
solinit = bvpinit(linspace(0, 1, 5), [0; 1]);
sol = bvp5c(odefun, bcfun, solinit);
y = deval(sol, 0.5);
fprintf('%.4f\n', y(1))

Saída

0.5697

Veja também

bvp4c bvpinit bvpset deval