function dydt=equ2(t,y,flag,N5,N6)
Taus=0.7520;
bbas=2*10^(-25);
bbes=3*10^(-25);
as=0.00035;
h=6.63*10^(-34);c=3*10^8;A=1.39*10^(-11);Lp=790*10^(-9);Ls=2*10^(-6);
n=1.45;
%P0=2*h*(c/n)^1.32/(Ls^3);
P0=0;
%N5=1*10^23;
%N6=4*10^23;
dydt=-(Taus*(bbes*N6-bbas*N5)-as)*y+Taus*bbes*N6*P0;