FEAT: in lamb1, add two functions to solve two-type Rayleigh equations

[Dengda98]

void grt_rayleigh1_roots(real_t nu, cplx_t y3[3]);
void grt_rayleigh2_roots(real_t m, cplx_t y3[3]);