function dy = func_camp(t,y)

global nu sigma k kt L q h lambda

a = y(1);
g = y(2);

dalpha = nu - sigma*phi(a,g);
dgamma = lambda*phi(a,g) - k*g;

dy = [dalpha;dgamma];
