t = linspace (0, 500, 1000);
x0 = [ 4; 1.1; 4 ];
t0 = cputime();
y = lsode ("oregonatorM", x0, t);
timeUsed = cputime()-t0


figure(1);
plot(t,y)
grid on

figure(2);
plot(t,y(:,2))
grid on

t0 = cputime();
y = lsode ("oregonator", x0, t);
timeUsed = cputime()-t0
