function dy=unidirectional_solve(t,y)

(-g*a/(a*tau-1)*exp(-t*(a*tau-1)/tau)+g/(a*tau-1)*a)*exp(-1/tau*t)