clear all
% set up global variables for use in the DE routine.
global r1 rm1 r2 rm2 r3

t0  = 0;     %
t1  = 5;    %
a0  = 10;     %
b0  = 10;   %
c0  = 10;
% initial state for use in DE solver.

x0  = [a0;b0;c0];


r1    = 1;
rm1   = 1;
r2    = 1;
rm2   = 20;
r3    = 1;



[t,xt] = ode45('GendronFirstMatLabProjectFunc',[t0,t1],x0);
plot(t,xt),legend('a','b','c')