function out = wavepinninggloballocal
out{1} = @init;
out{2} = @fun_eval;
out{3} = [];
out{4} = [];
out{5} = [];
out{6} = [];
out{7} = [];
out{8} = [];
out{9} = [];

% --------------------------------------------------------------------------
function dydt = fun_eval(t,kmrgd,k0,gam,k,del,hill,T)
dydt=[(T-kmrgd(2))*(k0+gam*kmrgd(1)^hill/(k^hill+kmrgd(1)^hill))-del*kmrgd(1);
(T-kmrgd(2))*(k0+gam*kmrgd(2)^hill/(k^hill+kmrgd(2)^hill))-del*kmrgd(2);];

% --------------------------------------------------------------------------
function [tspan,y0,options] = init
handles = feval(wavepinninggloballocal);
y0=[0,0];
options = odeset('Jacobian',[],'JacobianP',[],'Hessians',[],'HessiansP',[]);
tspan = [0 10];

% --------------------------------------------------------------------------
function jac = jacobian(t,kmrgd,k0,gam,k,del,hill,T)
% --------------------------------------------------------------------------
function jacp = jacobianp(t,kmrgd,k0,gam,k,del,hill,T)
% --------------------------------------------------------------------------
function hess = hessians(t,kmrgd,k0,gam,k,del,hill,T)
% --------------------------------------------------------------------------
function hessp = hessiansp(t,kmrgd,k0,gam,k,del,hill,T)
%---------------------------------------------------------------------------
function tens3  = der3(t,kmrgd,k0,gam,k,del,hill,T)
%---------------------------------------------------------------------------
function tens4  = der4(t,kmrgd,k0,gam,k,del,hill,T)
%---------------------------------------------------------------------------
function tens5  = der5(t,kmrgd,k0,gam,k,del,hill,T)
