solucion.sce

// Solución del modelo de Volterra
exec('data.sci', -1); // Carga de la data
Y= ode(Y_0, t_0, t, Volterra);
impares= [1:2:length(Y)]; // Índices impares
pares= [2:2:length(Y)]; // Índices pares
// Las x de la solución están los impares de Y
x= Y(impares);
// Las y de la solución están los pares de Y
y= Y(pares);
// Gráfico de x vs t
// Gráfico de y vs t
plot(t', x', "r", t', y', "b");
xtitle('Modelo de Volterra', 't', 'presa/depredador');
legend(['presas', 'depredadores']);