...
So I'm glad that issue was resolved. IF these results look good to you, one
One final thing I want to do is change the eccentricity for the plane of the ecliptic solution (I arbitrarily set it to 0.3=ecc) to see if there is a lower angular rate solution. decided to do was optimize for the smallest angular rate given some fixed inclination = 157.5 deg (23.5 deg giving the same solution), by varying the value of a first then varying the value of eccentricity. This is because the orbit solution for a = 12000 km, e = 0.3, inc = 23 deg was an unoptimized solution; when I wrote my own code for it, I was now able to optimize for the a and ecc to give the smallest angular rate. The plots for these results are shown below. (These solutions are implemented with the constraint that the perigee value should be at least 400 km above Earth).
The two solutions to take from this are:
eccentricity solution: (the ang_rate is off by a factor of 1000 because of the way units work in PoliAstro, the actual value is 87.22 arcsec/sec)
semi-major solution:
I must say that the two solutions look very similar to each other, with very similar eccentricities and semi major axes; this is expected, as I just parameterize one of the values. So, it will just be said that the solution is:
a = 14600, ecc = 0.53, inc = 157.5 or inc = 23.5