`bvpxtend`

Estende a solução de um problema de valor de contorno para uma nova extremidade do intervalo.

Categoria: Mathematics > Numerical Integration and Differential Equations

Sintaxe

solext = bvpxtend(sol, xnew, ynew)

Retorna um palpite no intervalo estendido até xnew, a partir de sol e do valor ynew no novo ponto.

Descrição

É a base da continuação: resolve-se um problema fácil em um intervalo curto, estende-se o palpite e resolve-se de novo em um intervalo maior, repetindo até chegar ao intervalo desejado. Ajuda em problemas difíceis de iniciar com um palpite grosseiro.

A extensão é apenas um palpite; ainda é preciso chamar o solver com as novas condições de contorno.

Argumentos de entrada

sol: Estrutura de solução de bvp4c ou bvp5c.
xnew: Novo ponto extremo, fora do intervalo atual.
ynew: Palpite para a solução em xnew (opcional).

Argumentos de saída

solext: Estrutura utilizável como argumento solinit de uma nova chamada do solver no intervalo mais longo.

Exemplo

Estendendo o intervalo

odefun = @(x, y) [y(2); -y(1)];
bcfun = @(ya, yb) [ya(1); yb(1) - 1];
sol = bvp4c(odefun, bcfun, bvpinit(linspace(0, 1, 5), [0; 1]));
solext = bvpxtend(sol, 2, [1; 0]);
solext.x(end)

Saída

ans = 2

Veja também

bvp4c bvp5c bvpinit odextend