function df_da = domega_da(a,b,t,w0)
df_da = (-b*w0*exp(b*t)*w0.*(1-exp(b*t)))


./((b+a*w0*(1-exp(b*t))).^2);