%
% cpldosc_fhat.m
%

function q =  cpldosc_fhat(v, p)

    v1        = v(1);
    v2        = v(2);
    sigma1    = p(1);
    sigma2    = p(2);
    omega     = p(3);

    q = zeros(2,1);
    q(1) = (-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);
    q(2) = (-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