Description: Mean longitude of the Moon minus mean longitude of the ascending node.
The source code for the Fortran routine FAF03 is available: