Is there a way to combine Sympy\'s solver with Astropy? For example, if I\'d like to solve Kepler\'s third law for p for given values of all other variables:
