FRC_PSI

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

Compute relative phase between forcing and the respnose frequency of a steady state

[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