(q(t)={(((((Subscript[q, I])-(Subscript[q, E])))/([1-exp(-α)]))[exp(-((α t)/(Subscript[t, f])))-exp(-α)]+(Subscript[q, E])t≤(Subscript[t, f])(Subscript[q, E])exp[-(OverHat[β])(t-(Subscript[t, f]))]t>(Subscript[t, f]),)Eq.   1)