The Stark problem applied to astrodynamics is a two-body model plus a constant perturbation acceleration.
This coefficient files and a basic Fortran wrapper for solving the Stark or Kepler problems using the F&G series approach along with a generalized Sundman transformation: LINK
Details on the F&G solution are found here:
Pellegrini, E., Russell, R.P., Vittaldev, V., “F and G Taylor Series Solutions to the Stark and Kepler Problems with Sundman Transformations,” Celestial Mechanics and Dynamical Astronomy, 118, No. 4, 2014, pp 355-378, DOI 10.1007/s10569-014-9538-7.
Also see:
Lantoine, G., Russell, R. P., “Complete, Closed-Form Solutions of the Stark Problem,” Celestial Mechanics and Dynamical Astronomy, Vol. 109, No. 4, 2011, pp. 333-366, DOI 10.1007/s10569-010-9331-1.
Hatten, N., Russell, R.P., “Comparison of Three Stark Problem Solution Techniques for the Bounded Case,” Celestial Mechanics and Dynamical Astronomy, 121, No. 1, 2015, pp 39-60, DOI 10.1007/s10569-014-9586-z.