y'(t)=alpha*y(t)/(b+2theta*t)-alpha*z(t)/(b+2theta*t);z'(t)=-alpha*y(t)/(b+2theta*(l-t))+alpha*z(t)/(b+2theta*(l-t))