odeoptions

Descrição

Esta função exibe interativamente um comando que deveria ser executado para se ajustar diversas opções para solucionadores de EDOs. A variável global %ODEOPTIONS ajusta as opções.

CUIDADO: a função ode verifica se esta variável existe e, neste caso, a utiliza. Para se utilizar valores padrões, deve-se limpar esta variável. Note que odeoptions não cria esta variável. Para criá-la, você deve executar a linha de comando exibida por odeoptions.

A variável %ODEOPTIONS é um vetor com os seguintes elementos:

[itask,tcrit,h0,hmax,hmin,jactyp,mxstep,maxordn,maxords,ixpr,ml,mu]

O valor padrão é:

[1,0,0,%inf,0,2,500,12,5,0,-1,-1]

O significado de cada elemento é descrito abaixo

itask 1 : computação normal nos tempos especificados 2 : computação nos tempos de malha (dados na primeira linha da saída de ode) 3 : um passo em um ponto interno da malha e retorna 4 : computação normal sem overshooting ("desvalorização excessiva") tcrit 5 : um passo, sem passagem de tcrit, e retorna

tcrit assume itask igual a 4 ou 5, descrito acima

h0 primeiro passo tentado

hmax tamanho máximo do passo

hmin tamanho mínimo do passo

jactype 0 : iterações funcionais, nenhum jacobiano usado ("adams" ou "stiff" apenas) 1 : jacobiano completo fornecido pelo usuário 2 : jacobiano completo internamente gerado 3 : jacobiano diagonal internamente gerado ("adams" ou "stiff" apenas) 4 : jacobiano com bandas fornecido pelo usuário (ver ml e mu abaixo) 5 : jacobiano com bandas internamente gerado (ver ml e mu abaixo)

maxordn máxima ordem não-stiff permitida, no máximo 12

maxords máxima ordem stiff permitida, no máximo 5

ixpr nível de impressão, 0 ou 1

ml,mu se jactype é igual a 4 ou 5, ml e mu são as semi-bandas inferior e superior do jacobiano com bandas: a banda são os i,j's com i-ml <= j <= ny-1. Se jactype é igual a 4 a função jacobiana deve retornar uma matriz J que é ml+mu+1 x ny (onde ny=dim de y em ydot=f(t,y)) tal que a coluna 1 de J é feita de mu zeros seguidos por df1/dy1, df2/dy1, df3/dy1, ... (1+ml entradas não-nulas possivelmente) a coluna 2 é feita de mu-1 zeros seguidos por df1/dx2, df2/dx2, etc

Ver Também

• ode — Solucionador de equações diferenciais ordinárias