FRC_ab

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

This function computes the autonomous reduced dynamics in polar coordinates as specified in SSMTool 1.

a = rho * real(lambda);
b = rho * (imag(lambda) - omega);

for j = 1:length(gamma)
    a = a + real(gamma(j))* rho.^(2*j+1);
    b = b + imag(gamma(j)) * (rho.^(2*j+1));
end

end