function q1 = cpldosc_fhat_q1(v1,v2,p)
    sigma1    = p(1);
    sigma2    = p(2);
    omega     = p(3);

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

   return