function f = omega(a,b,t,w0)
f = b*w0*exp(-b*t)./(b+a*w0*(1-exp(-b*t)));