@Markiyan Hirnyk 

All powers can be removed from an expression ex by evalindets(ex, `^`, x-> op(1,x)).

@Markiyan Hirnyk 

All powers can be removed from an expression ex by evalindets(ex, `^`, x-> op(1,x)).

@yewslyewsl Do you have reason to believe that there is a solution? You have two equations

Fn(n,m) = 0, Fm(n,m) = 0.

I can't find a shred of evidence that there is any point (n,m) that makes either one equal to zero, let alone one that makes them both equal zero. The evidence suggests that Fm is strictly positive and Fn is strictly negative. Mind you, I haven't proved anything; I've just looked at it empirically.

It looks like you may have tried to upload a file, but the link in your post is not clickable.

F:= subs(b= 1-a, eval(f)):
limit((F(x)-F(1))/(x-1), x= 1, right);
                               -1
limit((F(x)-F(1))/(x-1), x= 1, left);
                              2 a

F:= subs(b= 1-a, eval(f)):
limit((F(x)-F(1))/(x-1), x= 1, right);
                               -1
limit((F(x)-F(1))/(x-1), x= 1, left);
                              2 a

@Alex Joannou What version of Maple are you using? I think that Iterator requires Maple 17.

@Alex Joannou What version of Maple are you using? I think that Iterator requires Maple 17.

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.