function df_db = domega_db(a,b,t,w0)
df_db = ((b+a*w0*(1-exp(-b*t))).*(-b*w0*t.*exp(-b*t)+w0*exp(-b*t))-b*w0*exp(-b*t).*(1+a*w0*t.*exp(-b*t)))./((b+a*w0*(1-exp(-b*t))).^2);
