@Sphericalmoments Procedures for both the random bipartite tournament and all bipartite tournaments of a given size would be very easy to write, so I wouldn't quite call it a "package". I'd expect about 4 lines of code for the former and 7 for the latter.

I'd write them right now, but I had a minor muscle injury in my back 12 days ago that nonetheless is causing me incapacitating pain. But maybe I can do it later today.

Would you please define bipartite tournament? To my knowledge, a tournament is a directed graph constructed from an undirected complete graph by assigning a direction to each edge. I'm having trouble reconciling that with it being bipartite.

I wrote the above before you posted the code. I'll likely be able to figure it out from the code. But please attach your worksheet using the green uparrow on the editor's toolbar. It's much easier to read your code that way.

@AT Siacara As I said in my Answer below, a procedure made with the command proc cannot make indirect references to its parameters, and unapply should be used for such procedures. For example, you have

cstr5:= proc(NN, MM) FS - G1 end proc

where FS and/or G1 presumably contain references to NN and MM. That's what I mean by an indirect reference. It should be changed to 

cstr5:= unapply(FS - G1, [NN, MM])

I agree with you that the root cause is likely the same as with that other Question, but I haven't been able to figure out that root cause.

@ThomasLev If the order of subtraction of any one of the three square roots were reversed, then the integral would be a real number of the same magnitude. You should double check that you got those orders correct.

Can you give a brief defiintion or example of a "variation tab"?

(To avoid typographic confusion, I'm using L instead of l.) Starting with E being any positive number and using your definitions for R, L, and the interval of r, we have that r^2 < 4*E^2, so the integrand is guaranteed imaginary over the entire interval. For example, evaluating at the interval's endpoints, we get

restart:
R:= E/10:  L:= R/100:
eval~(r^2-4*E^2, r=~ [R-L, R+L]);

which is clearly negative.
 

@mmcdara If you're going to sort the output of combinat:-randperm, then you'd might as well just use combinat:-randcomb instead.

Sorry, I just accidentally separated this Reply of yours into a new Question. I can't undo that. But you can copy this back as a Reply (or, better, as an Answer).

You will make substantial progress simply by entering the code exactly as shown in the paper. People here are not likely to help you if you won't even do that. If you encounter some stumbling block, just attach your worksheet with error messages.

@Thomas Richard Your example makes it even weirder then: Sometimes assuming {} is okay, and sometimes it needs to be changed to assuming [].

I recommend experimenting with calling prefix-form `assuming`(expr, Alist).

@panke The plot3d command that you showed is missing multiplication signs and contains \pi, which is nonsense to Maple, although I understand perfectly well what you mean. It's not possible that this worked for you in some previous version of Maple or Windows.

However, of course, if Maple can't understand what you typed, it should just give you an error message instead of crashing. 

@KIRAN SAJJAN As far as I can see (I'm just reading on my phone at the moment, not looking in Maple), your system doesn't have parameters named Ha, S, etc. So how could those plots be made for your system?

By the way, I am very familiar with those plots because I wrote the worksheet that created them, and I can see in your Question that you're trying to adapt my worksheet.

@nm The OP is getting that strangely unsimplified result from a context menu, Tutors, rather than straightforward command entry.

I haven't looked at this in Maple yet, just speaking off-the-cuff: While sqrt(9) and 9^(1/2) mean the same thing mathematically, they are distinct Maple commands: sqrt(9) invokes the executable code sqrt (I think it's an appliable module) and returns 3, but 9^(1/2) remains unevaluated. This surely has something to do with the OP's problem, but I don't know whether the user can control this directly when using Tutors.

@tomleslie I just wrote a much-longer explanation as a Reply to my Answer above. The quoted help above is true to an extent, but it lacks any nuance. The last highligthed sentence applies across versions, not within a given session. Maple has gone to great lengths to ensure that re-executing code within one version, but possibly on different operating systems, will produce the same results (well, great lengths within the kernel; I can't say about the GUI). Any deviation is a reportable bug.

It's a bit specious for them to say "unordered sequence". As you say, it needs to be put in order at some point. That this happens the instant the set is created is easy to verify with addressof, disassemble, etc., and it's also mentioned in the Programming Guide.