Fixed:
1. The treatment of an internal variable bsi in tsolqu was wrong. It only affects the spin tracking anyway.

2. Unnecessary circular dependence of tfloor.f on tffs.f.

Added:
1. New optical functions GMX, GMY, GMZ to represent the gamma functions (1+alpha^2)/beta for matching, OpticsPlot and DRAW. They can be displayed by DISP GA and DISP Z. GMX and GMY are useful to estimate the natural chromaticities of a beam line.

Changes:
1. The usage of an internal function akang has been changed.

2. Some internal parameters p_.._QUAD and p_,,_MULT are changed.

3. Removed an unused arg poldiv of tmulbs.

4. SADHelp on SKn.

5. Rearrange some trigonometric/hyperbolic functions in tracking for efficiency. Also more use of hypot, p2h, etc.

6. Removed some computed GOTOs.

7. The number of iterations in teigen is increased to avoid inconvergence in some cases. Also merged tbal.f, thess.f, tqr.f into teigen.f as members of module eigen.
8. Inserting some sleeps in make, to avoid i/o error in the case of lxplus with eos.