clc; clear all; close all; 
format short
%parámetros de simulación
ti=0; h=0.0025; tf = 100; ts=ti:h:tf;  
opciones=odeset('RelTol',1e-3,'InitialStep', 2.5e-3, 'MaxStep',2.5e-3);
[t,x]=ode45('mra',ts,[0; 0],opciones);
plot(t,x(:,1),t,x(:,2))