global nu sigma k kt L q h

nu = 0.1;
%nu = 0.04;
sigma = 1.2;
k = 0.4;
kt = 0.4;
L = 10^6;
q = 100;
h = 10;

[T,Y] = ode23('func_camp',[0,2500],[92.366,10,2]);  

figure(1);
hold on;
plot(T,Y(:,1),'b');
plot(T,Y(:,2),'r');
plot(T,Y(:,3),'g');
xlabel('t');
ylabel('\alpha,\beta,\gamma');

