%this routine solves a system of odes

clear
global mu lambda kappa eps m m0

%specify parameter values
mu = 1;
lambda = 1;
kappa = 1;
eps = 0.1;
m0 = 0.001;
m = 2;



%specify the output points
tspan = [0:.1:500];

%specify initial data

y0 = [.01,0,.02];



[T,Y] = ode23('lac_deRHS',tspan, y0);  

plot(T,Y)