ode_discrete solucionador de equações diferenciais ordinárias, simulação de tempo discreto

## Seqüência de Chamamento

    y = ode("discrete",y0,k0,kvect,f)

## Parâmetros

y0

  vetor de reais ou matriz (condições iniciais).

t0

  escalar real (tempo inicial).

f

  função externa i.e. função, string ou lista.

k0

  inteiro (tempo inicial).

kvect

  vetor de inteiros.

## Descrição

Com esta sintaxe (primeiro argumento igual a `"discrete"`) `ode`, computa recursivamente `y(k+1)=f(k,y(k))` de um estado inicial `y(k0)` e retorna `y(k)` para `k` em `kvect`. `kvect(1)` deve ser maior do que ou igual a `k0`.

Outros argumentos e opções são os mesmos que para `ode`, veja a página de ajuda de ode.

## Exemplos

    y1 = [1;2;3];  deff("yp = a_function(k,y)", "yp = A*y + B*u(k)")
    A = diag([0.2,0.5,0.9]);  B = [1;1;1];  u = 1:10;  n = 5;
    y = ode("discrete", y1, 1, 1:n, a_function);
    y(:,2) - (A*y1+B*u(1))

    // Agora, y é avaliado em [y3,y5,y7,y9]
    y = ode("discrete", y1, 1, 3:2:9, a_function)

## Ver Também

- ode