# 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