// Equacoes: Caso de Pendulo Simples
// 
g = 1.0;
l = 1.0;
k = g/l;
deff('[dy]=exemp(t,y)',[
    'dy(1) = y(2);';
    'dy(2) = - k * sin(y(1));'])
// Abaixo temos que os dados serao gerados desde 0 com passo 0.01 ate 3.0
t=[0:0.01:12.0];
// Vetor das condicoes iniciais
y0 = [1; 0];
// Instante inicial
t0 = 0.0;
y=ode(y0,t0,t,exemp);
// Parada e preparacao da geracao de grafico
clf;
plot(t, y);