def METHOD_NAME(u,tt,mu,A,dhhdt,eta,t0,t1):
   tau = eta/mu
   W2 = dhhdt*tau*((tt/tau) - (u/alpha_fun(mu,A)) + (u/alpha_fun(mu,A))*(1 - t0*beta_fun(u,mu,eta,A))*math.exp(-beta_fun(u,mu,eta,A)*(tt - t0))*heaviside(tt-t0,-1)- ((tt/tau) - (u/alpha_fun(mu,A)) + (u/alpha_fun(mu,A))*(1 - t1*beta_fun(u,mu,eta,A))*math.exp(-beta_fun(u,mu,eta,A)*(tt - t1)))*heaviside(tt-t1,-1));
   return W2