// Equacoes: Caso Presa-Predador
// 
alfa1 = 1.0;
alfa2 = 0.1;
beta1 = 1.0;
beta2 = 0.1;

deff('[dy]=exemp(t,y)',[
    'dy(1) = (alfa1 -  alfa2 * y(2)) * y(1);';
    'dy(2) = (-beta1 + beta2 * y(1)) * y(2);'])
// Abaixo temos que os dados serao gerados desde 0 com passo 0.01 ate 3.0
t=[0:0.01:30.0];
// Vetor das condicoes iniciais
y0 = [1; 1];
// Instante inicial
t0 = 0.0;
y=ode(y0,t0,t,exemp);
// Parada e preparacao da geracao de grafico
clf;
plot(t, y);