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.

So, changing your colons to semicolons, we see that the solution set is a one-dimensional subset of C^3, but its intersection with R^3 is zero-dimensional.

So, changing your colons to semicolons, we see that the solution set is a one-dimensional subset of C^3, but its intersection with R^3 is zero-dimensional.

That's tricky as there are no parameters. I guess that the difference from regular solve is that it uses the SolveTools:-SemiAlgebraic algorithm, which in turn uses the massive RegularChains package.

It is also interesting that the solution is returned as a listlist.

That's tricky as there are no parameters. I guess that the difference from regular solve is that it uses the SolveTools:-SemiAlgebraic algorithm, which in turn uses the massive RegularChains package.

It is also interesting that the solution is returned as a listlist.

@GOODLUCK Why do your equations contain, repeatedly, expressions such as A/A, ctilde^s/ctilde^s?

@ysf Could you provide more details, please? Does it put the entire program in a single long line in the file? Could you upload your Maple code, please?

@ysf Could you provide more details, please? Does it put the entire program in a single long line in the file? Could you upload your Maple code, please?