%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% DOPRI 5/4 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% DOPRI 5/4
abstol = 10e-6;
reltol = 10e-6;
options = odeset('AbsTol',abstol,'RelTol',reltol);  % Options to ODE45
h = 1e-3;
x0 = [2;2];
tspan = [0,50];
a = 1;
b = 1;
[T,X] = DOPRI54(@PreyPredator,tspan, x0, h, abstol, reltol, a, b);

figure;
plot(X(:,1),X(:,2), 'b-x');
hold on
% Comparing with ode45
[T1, X1] = ode45(@PreyPredator, tspan, x0, options, a, b);

plot(X1(:,1), X1(:,2), 'r-x');
legend('DOPRI54','ODE45');
title('Embeded Dopri 5/4 Prey-Predator: a = 1, b = 1, x0 = [2;2].');
xlabel('X1');
ylabel('X2');