function q2 = cpldosc_fhat_q2(v1,v2,p)
    sigma1    = p(1);
    sigma2    = p(2);
    omega     = p(3);

    q2 = (-v2+exp(-40.0*v1+(8.0/15.0)).*tanh(v2*sigma2)-4.0*omega-v2.*exp(-40.0*v1+(8.0/15.0))+tanh(v2*sigma2)-v2*omega)./(exp(-40.0*v1+(8.0/15.0))+1.0);

   return
   