Increase NMAX to 1000
Calculate the precession
Add the tidal force by gas disk.
Add gas disk gravity