This might be a change in newer versions of Maple, but currently a with command performed inside a procedure has no effect on other commands in the procedure. But the execution of such a procedure causes the with to be performed globally. The procedures that are called within a procedure depend on what name-bindings are in effect at the time the procedure was defined, not the time that it is executed.

@Markiyan Hirnyk At the bottom of Markelov's LiveJournal post to which you gave a link there a link to a 1988 article in Russian: http://kvant.mccme.ru/1988/06/tri_formuly_ramanudzhana.htm. I don't read Russian, but it looks like this article attributes the formula to a Ramanujan. I am not sure. Could you check that please?

@Markiyan Hirnyk Thank you Markiyan for that stunningly beautiful identity. Where did you get that? It is reminiscent of Ramanujan.

The identity can indeed be proved in Maple. I feel that the proof warrants a separate post, so I am branching off. Put any replies there.

Actually, the theorem is more general than that: n can be a product of a power of 2 and any number of distinct Fermat primes.

But a proof of "A and B" is also a proof of "A" alone, so I don't think that it is fair to nitpick when someone cites a specific case of a more-general theorem. I was just using the case that I needed to make my point.

All of your questions are answered in detail in my Answer in this thread.

All of your questions are answered in detail in my Answer in this thread.

Both of your questions are answered in detail in my Answer.

Preben, your modifications work because they both return a two-member list even when they are passed symbolic arguments. See my much-more-detailed Answer.

Preben, your modifications work because they both return a two-member list even when they are passed symbolic arguments. See my much-more-detailed Answer.

@Sergio Parreiras I am eager to hear how any of your explorations with Maple IDE go.

Here's another basic debugging command to add to my list above:

infolevel[all]:= 2;

The infolevel and printlevel can be adjusted to various levels.

From running your code with infolevel set, I see that some large polynomials in β are being factored, perhaps one of your input polynomials. I'd be interested in running the code on the individual factors to see if one of them generates the bug.

@Sergio Parreiras I am eager to hear how any of your explorations with Maple IDE go.

Here's another basic debugging command to add to my list above:

infolevel[all]:= 2;

The infolevel and printlevel can be adjusted to various levels.

From running your code with infolevel set, I see that some large polynomials in β are being factored, perhaps one of your input polynomials. I'd be interested in running the code on the individual factors to see if one of them generates the bug.

Is G really an Array (with a capital A), or is it a list? What about K? And what are their sizes? Is it 14 and 4?

I'd need to see the rest of the code, where you define BC, XY, and SOL. How about uploading a worksheet?

I'll sign up for this beta testing if at least one person on this list with a reputation higher than mine also signs up. So, if you're signing up, please let me know. Send me email if you want to keep it private.

Or, if any such person has a specific reason why they are not signing up (like, more specific than I'm too busy), please let me know.

I don't understand this line in your code:

B := select(A->abs(A[4])B[5,1]^2+B[5,2]^3+B[5,3]^4; #check

Could you either explain it or correct it?