You need more Digits in order to get an accurate answer to this. The inaccuracy of your answer is obscured because you checked the residuals in the evalfed form of the equation rather than in the original. See my solution.

You need more Digits in order to get an accurate answer to this. The inaccuracy of your answer is obscured because you checked the residuals in the evalfed form of the equation rather than in the original. See my solution.

@mjc Are you sure that you used a forward slash "/" instead of a backward slash "\" in the libname setting?

@mjc Are you sure that you used a forward slash "/" instead of a backward slash "\" in the libname setting?

@mjc What happens if you try to load the package with(Orbitals);

@mjc What happens if you try to load the package with(Orbitals);

I had an Answer with that rule, but I deleted it because of this weirdness, which I couldn't come up with any way to prevent:

applyrule(r, (a-d)*ee);

Note that this result does not contain d.

I had an Answer with that rule, but I deleted it because of this weirdness, which I couldn't come up with any way to prevent:

applyrule(r, (a-d)*ee);

Note that this result does not contain d.

@brian bovril I just used the Tools => Options menu, nothing from the Format => Styles menu. The 12-point Courier New bold brown is the default font for character style Maple Input.

What happens if you select Maple Input from the pull-down style selector, which is the leftmost item in the bottom row of the toolbar?

Are you creating your worksheets in Document mode or Worksheet mode? You might need to use Worksheet mode.

@brian bovril I just used the Tools => Options menu, nothing from the Format => Styles menu. The 12-point Courier New bold brown is the default font for character style Maple Input.

What happens if you select Maple Input from the pull-down style selector, which is the leftmost item in the bottom row of the toolbar?

Are you creating your worksheets in Document mode or Worksheet mode? You might need to use Worksheet mode.

@Markiyan Hirnyk The solution set over the complex numbers has one free parameter (to be expected since there are three unknowns and only two equations) whereas the solution set over the real numbers has zero free parameters.

@Markiyan Hirnyk The solution set over the complex numbers has one free parameter (to be expected since there are three unknowns and only two equations) whereas the solution set over the real numbers has zero free parameters.

@Markiyan Hirnyk That's because Alejandro's solution seemed to illustrate a general technique that could be applied to many problems. Indeed, solve(sys, vars, parametric, real) may be a complete replacement for the totally unreliable RealDomain:-solve when sys is polynomial.

Pagan's solution, on the other hand, seemed highly specialized because of the reason that I mentioned in my comment to that solution.

@Markiyan Hirnyk That's because Alejandro's solution seemed to illustrate a general technique that could be applied to many problems. Indeed, solve(sys, vars, parametric, real) may be a complete replacement for the totally unreliable RealDomain:-solve when sys is polynomial.

Pagan's solution, on the other hand, seemed highly specialized because of the reason that I mentioned in my comment to that solution.

@EWI If you remove the inequalities, solve will quickly return nine solution sets for the equations.