// Equacoes: exemplo de massa-mola
// Duas molas, uma massa
k = [1.0; 1.0]; // Constantes das molas
// Se os valores de k forem ambos diferentes de zero: Molas presas nas paredes
// Se um valor de k for nulo: Este lado esta' "solto"
deff('[dy]=exemp(t,y)',[
    'dy(1) = y(2);';
    'dy(2) = - (k(1) + k(2)) * y(1);'])
// 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; 0];
// Respectivamente: condicao para a velocidade e posicao
// Instante inicial
t0 = 0.0;
// Chamada da funcao ode() sem especificar o metodo a ser utilizado
y=ode(y0,t0,t,exemp);
// Parada e preparacao da geracao de grafico
clf;
plot(t, y);