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);
