FRC Phase

function psi = frc_psi(rho,omega,gamma,lambda,f)

This function analytically computes the phase difference of the response with respect to the external forcing for a leading order non-autonomous approximation of the reduced dynamics.

[a,b] = frc_ab(rho,omega,gamma,lambda);

psi = atan2( (rho.*b * real(f) - a*imag(f)), (-a*real(f) - rho.*b*imag(f)) );

end