NAME:
  kepler1
 PURPOSE:   (one line only)
  Solve Kepler's Equation (small eccentricity)
 DESCRIPTION:
  Solve Kepler's Equation. Method by S. Mikkola (1987) Celestial
     Mechanics, 40 , 329-334. 
    result from Mikkola then used as starting value for
       Newton-Raphson iteration to extend the applicability of this
       function to higher eccentricities, Courtesy of Joern Wilms.
 CATEGORY:
  Celestial Mechanics
 CALLING SEQUENCE:
  kepler1,m,ecc,eccanom
 INPUTS:
  m    - Mean anomaly (radians, can be an array)
  ecc  - Eccentricity
 OPTIONAL INPUT PARAMETERS:
 KEYWORD INPUT PARAMETERS:
  THRESH -  stopping criterion for the Newton Raphson iteration; the
            iteration stops once abs(E-Eold)