Note: I am using
2*x+cos(x) as an example below. I'm asking this question for any function where
Solve doesn't yield a closed form result. My specific interest at the moment is getting Mathematica to spit out the (well known) power series expansion solution of Kepler's Equation of the Center, just as a test.
2*x+cos(x) is monotonic increasing (derivative is
2-sin(x)), and thus invertible. There's no simple form for the inverse:
Solve[2*x+Cos[x] == y, x] Solve::nsmet: This system cannot be solved with the methods available to Solve.
Question: can I get a power series approximation to the solution? Either using Mathematica builtin functions or some sort of voodoo?