function dy = func_camp(t,y)

global nu sigma k kt L q h

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

dalpha = nu - sigma*phi(a,g);
dbeta = q*sigma*phi(a,g) - kt*b;
dgamma = kt*b/h - k*g;

dy = [dalpha;dbeta;dgamma];
