@nm 

There is nothing wrong with having text or any other data structure in a matrix. Indeed, certain commands in Maple require an Array, Matrix, or Vector of plots. I understand your point about "lists of lists," but note that a Maple Matrix is different from a Maple list of lists.

There is nothing wrong with having keys in a set.

@mzaman 

My guess that using the timelimit wouldn't be worthwhile was correct. But I include the technique here to show how to use the command. In the worksheet below, I call it with the minimal amount of time, .01 seconds. I tried it at 5 seconds, but the program ran for 40 minutes without completing any factorizations.

I also improved the code for selecting the kp2 closest to 56 bits. In the run below, you can see that p is 511 bits and q is 512 bits.

Download WuSun.mw

@gdoug 

What version of Maple are you using? The above works for me, even with a single vector, in Maple 16 through 2015.

Also, did you cut and paste the code exactly as I had it? Or did you retype it? If you retyped it, please post the retyped version here.

@Yiannis Galidakis 

Those are wonderful animations. Could you post the Maple code to generate them?

And your Reply above would make an excellent Post. Please consider Posting it.

What makes you think that the precision is currently set at 100 digits? The default is 10. To control the number of digits, set the Digits environment variable. For example,

Digits:= 15:

However, I doubt that reducing the number of digits will do much to speedup your problem, even if your precision is currently 100.

How about posting your code?

@J4James 

What "issue" do you mean? And what is the significance of the plot that you posted? It has three curves. The sum from 10001 to infinity (the part that I truncated) is much too small to make any difference in the plot!

Could you please provide the complete code, such as the definitions of bc_cond, slope_cond, etc.? You can upload a Maple worksheet by using the green uparrow, which is the last item on the second row of the toolbar in the MaplePrimes editor.

By "use procedure programming structure in order to get the same results," do you mean that you want to use something like Newton's method?

@GPY 

It does invoke the function; it just doesn't show you the results. The results of statements nested two or more levels deep are by default not displayed. If you want to see it, then you can wrap the invocation in print(...), as in print(cubesum(cubeprod)).

And, yes, the check iroot(cubeprod,3) = i is redundant.

@nm 

Your screen shot shows that you are using Adode Acrobat 9 Pro, which is a PDF reader and writer, not just a PDF reader, such as Adobe Postscript. And, I don't think that Acrobat is freely available.

@Chia 

Your eval statement contains the following list of variables:

That list should contain only the boolean variables, and there should be nv of them.

@Yiannis Galidakis 

You do have a reputation of 10, which is enough to vote up Answers. Also, it does not require any reputation to select an Answer as the best Answer to one's own Question.

@GPY 

Yes, I noticed that you started a new thread for this question. I am about to Answer your Question in that thread. So please end discussion of it here in this thread.

@GPY 

It's an interesting problem. Yes, I believe that it would work. Note that you'd be finding solutions to x^3+y^3 = n, not x^3+y^3 = n^3. The command analogous to isqrt is iroot. There's some question as to whether it's worth using an alternative to x3:= x^3.

@Yiannis Galidakis 

What you want to do with animations can't be done as a frame-by-frame animation. It's a severe limitation of Maple animation. There will always be one fixed set of axes that applies to every frame of the animation.

However, you can use the viewpoint option of plot3d to "get closer to" or "fly by" or even "fly through" a surface. This essentially creates an animation with a single frame viewed from multiple positions. See ?plot3d,viewpoint.