Carl Love

Carl Love

26488 Reputation

25 Badges

12 years, 261 days
Himself
Wayland, Massachusetts, United States
My name was formerly Carl Devore.

MaplePrimes Activity


These are answers submitted by Carl Love

In your first expression, TPP, you have square brackets used algebraically; in other words, they're used as a substitute for parentheses. While that is indeed standard in mathematical typesetting, and I agree that it improves readibility, it's not alllowed in Maple (nor in any other computer language that I'm aware of). The reason is that with the severe limitation of the number of special characters available on a standard keyboard, they need to all have their own meanings.

If I replace those square brackets with paremtheses, the error goes away.

In relatively recent versions of Maple, and only using 1-D input, a for-do loop (or, indeed, any do loop) can be used as an embedded expression whose output is an expression sequence:

A:= select(isprime, {$3..101});
    A := {3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 
      59, 61, 67, 71, 73, 79, 83, 89, 97, 101}

B:= {for a in A do if a mod 3 = 1 then a fi od};
    B := {7, 13, 19, 31, 37, 43, 61, 67, 73, 79, 97}

This is not just a syntactic nicety; it's often the most efficient way to solve set-building problems. However, in this case, acer's select works efficiently also, and it's what I would use.

The difference between what you perceive as the correct and incorrect results is only the constant 3/4. Being different by a constant is not incorrect for an indefinite integral.

Here's a solution that allows a decision to be made based on the number of arguments and avoids any awkwardness that may come from having arguments that don't match named parameters due to there not being any named parameters. An overload is a list of (usually anonymous) procedures all assigned to a single procedure name such that for any call to the procedure, Maple will choose the first member of the list whose declared parameters match the arguments in the call.

The purpose of the (called the end-of-parameters marker) in the following example is to force those procedures to balk at (rather than the default, to ignore) extra arguments. To advance to the next procedure in the list, it's necessary that the arguments not match the declared parameters; I think that no other error will cause it to advance.

There seemed to be no good reason to limit this to real constants (type realcons), so I used type algebraic. But it would work just as well if you changed those to realcons..

P:= overload([
    proc(a::algebraic, b::algebraic, $) option overload; 
        a/b 
    end proc,

    proc(a::algebraic, b::algebraic, c::algebraic, $) option overload;
        a + b*c
    end proc,
    
    #Fall-through case: You can choose to return an error, to return unevaluated (as shown here), 
    #or to leave off this case entirely:
    ()-> 'procname'(args)
]):

P(3,5);
                               3/5
P(7,1,1,1);
                         P(7, 1, 1, 1)
P(3,5,7);
                               38

 

It's just an undocumented feature. There are many commands in Maple that work on a wider range of types of input than they are documented to work on. If they documented it, they'd never be able to change it in the future. On the other hand, it's often convenient to treat a listlist as a Matrix.

Just about any procedure can be made into an operator-form procedure by giving it a name beginning with (see ?&). This is most natural with procedures of 2 arguments, but it's not limited to them. For example,

`&<--`:= (a, b::{list,set}(`=`))-> eval['recurse'](a,b):
(x+3) &<-- [x= 4];

                               
7

&-operators have a very high precedence among infix operators (see ?operators,precedence), so you need to be wary of correct usage of parentheses and other grouping symbols.

Regarding a comparison of eval['recurse'] and MmaTranslator:-Mma:-ReplaceRepeated: The code for both of these is short and easy to read. The actual procedure name for eval[recurse] is `simpl/eval` (it's actually an appliable module rather than a procedure, but that's an insignificant technicality here). So,

showstat(`simpl/eval`);
showstat(MmaTranslator:-Mma:-ReplaceRepeated);
 

The concept is weirdly named. I'd call it the UNbinding number. (I think some mathematician had an over-sexualized imagination when naming it the bondage number.) Anyway, for your example 7-regular graph with 48 automorphisms, I can very easily find a set of 6 edges whose removal increases the domination number from 4. It's more difficult to prove that no set of 5 edges suffice (binomial(84,5) ~ 31 million). By brute-force all-subsets search, I've proved that no 4-subset of edges will work. Thus, the "bondage number" is 5 or 6.

There are only 300 (exactly) minimal dominating sets. If you go back to the thread on those, I've posted new code that returns the set of 300. For each candidate set of edges, it's much easier computationally to check whether any of these 300 is still a dominating set after the edges are removed than it is to recompute the domination number from scratch. I'll post some code here along those lines.

Regarding your comment on symmetry (primarily from your linked StackExchange thread): My intuition is that the symmetry will have a profound effect on the computational effort required by various algorithms to compute the bondage number. In my new Answer in the domination-number thread, I've shown how to partition the 300 minimal dominating sets into 12 equivalence classes under automorphism (two sets being considered equivalent if there's an automorphism mapping one to the other). This partitioning takes less than one second, including the time required to find the 48 automorphisms and convert them to vertex permutations. My intuition is that these equivalence classes could somehow be used to reduce the computation needed to get the bondage number, but I haven't been able to make any further developments on that idea. And it may be an average-case reduction rather than a worst-case reduction.

From your example, it looks like perhaps by "similar" you mean "identical".  If that is the case, then Maple already recognizes identical subexpressions and makes use of that in its internal storage. It's only the displayed form of expressions that don't show this. 

Here's a procedure that will show you the common subexpressions that Maple has recognized:

CommonSubExprs:= e->
    select(se-> numboccur(expr, se) > 1, indets(e, {`+`, `*`, `^`, 'function'}))
:

The idea can be easily extended to more types than I've used above. I was just guessing that  `+`, `*`, `^`, function would cover what you need.

Any of the expressions returned by the procedure can be substituted by a variable with the ordinary subs command.

You just need to use the given ordering as an index of the vertex list:

R:= IsChordal(s, eliminationordering):
if [R][1] then op(3,s)[R[2]] fi;

If by Gamma you mean the generalized factorial function, then you need to spell it GAMMA (all uppercase). The cause of the error is that it doesn't know what Gamma means,

I improved my algorithm from above (made it faster) and added extensive explanation of how it works. Considering the code you've posted recently, you desperately need to study this explanation.

Your example graph has a lot of symmetry, 48 automorphisms. In the worksheet below, I show some techniques for analyzing the symmetry, including dividing the 300 minimal dominating sets into 12 equivalence classes under the automorphisms. The equivalence relation is specifically that two minimal dominating sets are equivalent if there is an automorphism of mapping one to the other. 
 

Minimal dominating sets of a graph

 

restart
:

#Author: Carl Love <carl.j.love@gmail.com>, the Glorious 1st of June, 2023
#
#Returns all minimal dominating sets of the input simple graph.
#
                                         #Reference #s for expository comments:
MinDominatingSets:= proc(G::Graph, {any::truefalse:= false})                # 1
local                         
    n:= nops(op(3,G)), V:= {$n},                                            # 2
    A:= [seq](op(4,G) union~ `{}`~(V)),                                     # 3
    N:= table([()= {}]), N0,                                                # 4
    R:= <()>,                                                               # 5
    s, Ns, v                                                                # 6
;
    do                                                                      # 7
        N0:= eval(N); N:= evaln(N);                                         # 8
        for s,Ns in eval(N0) do                                             # 9
            for v from `if`(s::posint, s, `if`(s=(), 0, s[-1])) + 1 to n do #10
                if (N[s,v]:= Ns union A[v]) = V then                        #11
                    if any then return {s,v} else R,= {s,v} fi              #12
                fi                    
            od
        od;
        if numelems(R) <> 0 then return {seq}(R) fi                         #13
    od    
end proc
:

Expository comments on the code and algorithm:

 

#0: Note that the procedure does not call any command from GraphTheory. There's usually no need to do that. Indeed, this procedure doesn't use any commands from any package.

 

#1: If the keyword option any is passed, then the procedure returns only the first minimal dominating set that it finds.

 

#2: One key to efficient graph theory algorithms is to never refer to vertices by their labels. Instead, they are always referred to by their index into the list of vertices. The list of vertices is op(3, G), the same as GraphTheory:-Vertices(G). So, V here is simply the set {$1..n} where n is the number of vertices. All further references to vertices, including within edges, will use these indices.

 

#3: The second key to efficient graph theory algorithms (this is very important) is what I call the adjacency array. I don't know if it has a formal name. It's returned by op(4, G). It's a sparse representation of the adjacency matrix. It's a one-dimensional array A indexed by the vertices such that A[k] is the set of vertices adjacent to k.

 

For the purpose of calculating dominating sets, it's convenient to consider vertices as being adjacent to themselves; so, the union~ adds each k to the set A[k].

 

#4: N is a table that is the key data structure of this algorithm. The kth iteration of the main loop (#7) constructs all k-susbsets of V and their neighborhoods. The k-subsets (stored as sequences) are the indices of  N, and the neighborhoods are the entries. They are constructed from the (k-1)-subsets and their neighborhoods, which are moved to table N0 at at the beginning of loop (#7).

 

#5: R is a array to hold the dominating sets to be returned.

 

#6: On the kth iteration of the main loop (#7), s will be a (k-1)-subset of vertices as a sequence, Ns will the neighborhood of s, and v be a vertex to be added to s to make a k-subset.

 

#7: Everything that needs to be said about this line is covered in #4 and #6 above. Although I've referred to the kth iteration of this loop several times, the code itself does not use k.

 

#8: Move the (k-1)-subsets and their neighborhoods to table N0, then clear and garbage-collect N in preparation for the k-subsets.

 

#9: When a for loop has two control variables, the first is the index into the container, and the second is the container's element at that index. So, in this case, s runs over all indices of N0, and Ns = N0[s].

 

#10: Add to s all vertices v in such a way that each subset of V is constructed at most once. Specifically, we only add vertices v that are greater than all vertices is s.

 

#11: The index of N is the new k-subset (stored as a sequence); the entry is its neighborhood. If the neighborhood is all of V, then the index is a dominating set.

 

#12: If we have a dominating set and keyword any is set, then return it; otherwise save the dominating set for later return. The R,= is an efficient way to append a new element to an Array without indexing.

 

#13: If any dominating sets have been stored in R, then return them. They're necessarily all dominating sets of size k and all minimal.

 

Note that a return statement can be used anywhere in a procedure. You do not need to exit loops before using return.

 

#Example: A 7-regular graph with 24 vertices
#
g:= (GT:= GraphTheory):-Graph(
    [$24],
    {
        {1, 2}, {1, 3}, {1, 4}, {1, 5}, {1, 7}, {1, 9}, {1, 11},
        {2, 3}, {2, 5}, {2, 6}, {2, 7}, {2, 12}, {2, 14},
        {3, 4}, {3, 7}, {3, 8}, {3, 9}, {3, 16},
        {4, 5}, {4, 9}, {4, 10}, {4, 11}, {4, 17},
        {5, 6}, {5, 11}, {5, 12}, {5, 19},
        {6, 7}, {6, 12}, {6, 13}, {6, 14}, {6, 21},
        {7, 8}, {7, 14}, {7, 15},  
        {8, 9}, {8, 14}, {8, 15}, {8, 16}, {8, 22},
        {9, 10}, {9, 16}, {9, 17},
        {10, 11}, {10, 17}, {10, 18}, {10, 19}, {10, 23},
        {11, 12}, {11, 18}, {11, 19},  {12, 13}, {12, 19}, {12, 20},
        {13, 14}, {13, 19}, {13, 20}, {13, 21}, {13, 24},
        {14, 15}, {14, 21},  {15, 16}, {15, 21}, {15, 22}, {15, 24},
        {16, 17}, {16, 22}, {16, 23},  {17, 18}, {17, 22}, {17, 23},
        {18, 19}, {18, 20}, {18, 23}, {18, 24},  {19, 20},
        {20, 21}, {20, 23}, {20, 24},  {21, 22}, {21, 24},
        {22, 23}, {22, 24}, {23, 24}
    }
);

GRAPHLN(undirected, unweighted, [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24], Array(1..24, {(1) = {2, 3, 4, 5, 7, 9, 11}, (2) = {1, 3, 5, 6, 7, 12, 14}, (3) = {1, 2, 4, 7, 8, 9, 16}, (4) = {1, 3, 5, 9, 10, 11, 17}, (5) = {1, 2, 4, 6, 11, 12, 19}, (6) = {2, 5, 7, 12, 13, 14, 21}, (7) = {1, 2, 3, 6, 8, 14, 15}, (8) = {3, 7, 9, 14, 15, 16, 22}, (9) = {1, 3, 4, 8, 10, 16, 17}, (10) = {4, 9, 11, 17, 18, 19, 23}, (11) = {1, 4, 5, 10, 12, 18, 19}, (12) = {2, 5, 6, 11, 13, 19, 20}, (13) = {6, 12, 14, 19, 20, 21, 24}, (14) = {2, 6, 7, 8, 13, 15, 21}, (15) = {7, 8, 14, 16, 21, 22, 24}, (16) = {3, 8, 9, 15, 17, 22, 23}, (17) = {4, 9, 10, 16, 18, 22, 23}, (18) = {10, 11, 17, 19, 20, 23, 24}, (19) = {5, 10, 11, 12, 13, 18, 20}, (20) = {12, 13, 18, 19, 21, 23, 24}, (21) = {6, 13, 14, 15, 20, 22, 24}, (22) = {8, 15, 16, 17, 21, 23, 24}, (23) = {10, 16, 17, 18, 20, 22, 24}, (24) = {13, 15, 18, 20, 21, 22, 23}}), `GRAPHLN/table/1`, 0)

#Get 1 minimal dominating set:
ds:= CodeTools:-Usage(MinDominatingSets(g, any));

memory used=1.01MiB, alloc change=0 bytes, cpu time=0ns, real time=13.00ms, gc time=0ns

{1, 11, 21, 22}

#Get them all:
ds:= CodeTools:-Usage(MinDominatingSets(g));

memory used=13.89MiB, alloc change=37.00MiB, cpu time=110.00ms, real time=80.00ms, gc time=62.50ms

{{1, 2, 19, 22}, {1, 5, 14, 23}, {1, 6, 8, 18}, {1, 6, 15, 18}, {1, 6, 16, 18}, {1, 6, 18, 22}, {1, 6, 19, 22}, {1, 7, 13, 17}, {1, 7, 13, 23}, {1, 7, 17, 20}, {1, 7, 19, 22}, {1, 7, 20, 23}, {1, 8, 10, 13}, {1, 8, 13, 17}, {1, 8, 13, 18}, {1, 8, 13, 23}, {1, 9, 13, 24}, {1, 9, 20, 21}, {1, 10, 12, 15}, {1, 10, 13, 15}, {1, 10, 13, 16}, {1, 10, 13, 22}, {1, 11, 13, 16}, {1, 11, 13, 22}, {1, 11, 14, 23}, {1, 11, 16, 21}, {1, 11, 21, 22}, {1, 12, 14, 23}, {1, 12, 15, 17}, {1, 12, 15, 18}, {1, 12, 15, 23}, {1, 13, 14, 17}, {1, 13, 14, 23}, {1, 13, 15, 17}, {1, 13, 15, 18}, {1, 13, 15, 23}, {1, 13, 16, 17}, {1, 13, 16, 18}, {1, 13, 16, 19}, {1, 13, 16, 23}, {1, 13, 17, 22}, {1, 13, 18, 22}, {1, 13, 19, 22}, {1, 13, 22, 23}, {1, 14, 17, 20}, {1, 14, 19, 22}, {1, 14, 19, 23}, {1, 14, 20, 23}, {1, 16, 19, 21}, {1, 19, 21, 22}, {2, 3, 10, 21}, {2, 3, 10, 24}, {2, 3, 18, 21}, {2, 3, 18, 24}, {2, 3, 19, 22}, {2, 4, 8, 20}, {2, 4, 15, 20}, {2, 4, 16, 20}, {2, 4, 19, 22}, {2, 4, 20, 22}, {2, 5, 9, 24}, {2, 8, 10, 13}, {2, 8, 10, 20}, {2, 8, 10, 21}, {2, 8, 10, 24}, {2, 9, 10, 21}, {2, 9, 10, 24}, {2, 9, 11, 24}, {2, 9, 12, 24}, {2, 9, 18, 21}, {2, 9, 18, 24}, {2, 9, 19, 22}, {2, 9, 19, 24}, {2, 10, 12, 15}, {2, 10, 12, 22}, {2, 10, 13, 15}, {2, 10, 13, 16}, {2, 10, 13, 22}, {2, 10, 14, 23}, {2, 10, 15, 19}, {2, 10, 15, 20}, {2, 10, 15, 21}, {2, 10, 15, 24}, {2, 10, 16, 20}, {2, 10, 16, 21}, {2, 10, 16, 24}, {2, 10, 19, 22}, {2, 10, 20, 22}, {2, 10, 21, 22}, {2, 10, 22, 24}, {2, 11, 13, 16}, {2, 11, 16, 20}, {2, 11, 16, 21}, {2, 11, 16, 24}, {2, 12, 15, 17}, {2, 12, 17, 22}, {2, 14, 17, 18}, {2, 15, 17, 19}, {2, 17, 19, 22}, {3, 4, 13, 24}, {3, 4, 20, 21}, {3, 5, 10, 21}, {3, 5, 14, 23}, {3, 5, 17, 21}, {3, 5, 18, 21}, {3, 5, 21, 23}, {3, 6, 8, 18}, {3, 6, 10, 21}, {3, 6, 10, 24}, {3, 6, 15, 18}, {3, 6, 16, 18}, {3, 6, 18, 21}, {3, 6, 18, 22}, {3, 6, 18, 24}, {3, 6, 19, 22}, {3, 7, 19, 22}, {3, 10, 11, 21}, {3, 10, 12, 15}, {3, 10, 12, 21}, {3, 10, 19, 21}, {3, 11, 13, 16}, {3, 11, 13, 22}, {3, 11, 14, 23}, {3, 11, 16, 21}, {3, 11, 17, 21}, {3, 11, 18, 21}, {3, 11, 21, 22}, {3, 11, 21, 23}, {3, 12, 14, 23}, {3, 12, 15, 17}, {3, 12, 15, 18}, {3, 12, 15, 23}, {3, 12, 17, 21}, {3, 12, 18, 21}, {3, 12, 21, 23}, {3, 13, 16, 19}, {3, 13, 19, 22}, {3, 14, 19, 22}, {3, 14, 19, 23}, {3, 16, 19, 21}, {3, 17, 19, 21}, {3, 18, 19, 21}, {3, 19, 21, 22}, {3, 19, 21, 23}, {4, 5, 8, 20}, {4, 5, 8, 24}, {4, 5, 14, 23}, {4, 5, 15, 20}, {4, 5, 15, 24}, {4, 6, 8, 18}, {4, 6, 8, 20}, {4, 6, 15, 18}, {4, 6, 15, 20}, {4, 6, 16, 18}, {4, 6, 16, 20}, {4, 6, 18, 22}, {4, 6, 19, 22}, {4, 6, 20, 22}, {4, 7, 8, 20}, {4, 7, 13, 17}, {4, 7, 13, 23}, {4, 7, 15, 20}, {4, 7, 16, 20}, {4, 7, 17, 20}, {4, 7, 19, 22}, {4, 7, 20, 22}, {4, 7, 20, 23}, {4, 8, 12, 20}, {4, 8, 12, 24}, {4, 8, 14, 20}, {4, 10, 12, 15}, {4, 11, 14, 23}, {4, 12, 14, 23}, {4, 12, 15, 17}, {4, 12, 15, 18}, {4, 12, 15, 20}, {4, 12, 15, 23}, {4, 12, 15, 24}, {4, 13, 14, 17}, {4, 13, 14, 23}, {4, 14, 15, 20}, {4, 14, 16, 20}, {4, 14, 17, 20}, {4, 14, 19, 22}, {4, 14, 19, 23}, {4, 14, 20, 22}, {4, 14, 20, 23}, {5, 6, 8, 18}, {5, 6, 8, 23}, {5, 6, 9, 24}, {5, 6, 16, 18}, {5, 6, 16, 23}, {5, 7, 9, 24}, {5, 7, 13, 17}, {5, 7, 17, 20}, {5, 7, 17, 21}, {5, 7, 17, 24}, {5, 8, 9, 20}, {5, 8, 9, 24}, {5, 8, 10, 13}, {5, 8, 10, 20}, {5, 8, 10, 21}, {5, 8, 10, 24}, {5, 8, 13, 17}, {5, 8, 13, 18}, {5, 8, 13, 23}, {5, 8, 14, 18}, {5, 8, 14, 23}, {5, 8, 17, 20}, {5, 8, 17, 21}, {5, 8, 17, 24}, {5, 8, 18, 20}, {5, 8, 18, 21}, {5, 8, 18, 24}, {5, 8, 19, 22}, {5, 8, 20, 23}, {5, 8, 21, 23}, {5, 8, 23, 24}, {5, 9, 14, 23}, {5, 9, 14, 24}, {5, 9, 15, 20}, {5, 9, 15, 24}, {5, 14, 16, 18}, {5, 14, 16, 23}, {5, 15, 16, 19}, {6, 7, 10, 23}, {6, 7, 17, 18}, {6, 8, 9, 18}, {6, 8, 11, 18}, {6, 8, 11, 23}, {6, 9, 10, 21}, {6, 9, 10, 24}, {6, 9, 11, 24}, {6, 9, 12, 24}, {6, 9, 15, 18}, {6, 9, 16, 18}, {6, 9, 18, 21}, {6, 9, 18, 22}, {6, 9, 18, 24}, {6, 9, 19, 22}, {6, 9, 19, 24}, {6, 11, 13, 16}, {6, 11, 16, 18}, {6, 11, 16, 20}, {6, 11, 16, 21}, {6, 11, 16, 23}, {6, 11, 16, 24}, {7, 9, 11, 24}, {7, 9, 12, 24}, {7, 9, 19, 22}, {7, 9, 19, 24}, {7, 10, 12, 15}, {7, 10, 12, 22}, {7, 10, 15, 19}, {7, 10, 19, 22}, {7, 11, 13, 16}, {7, 11, 13, 17}, {7, 11, 16, 20}, {7, 11, 16, 21}, {7, 11, 16, 24}, {7, 11, 17, 20}, {7, 11, 17, 21}, {7, 11, 17, 24}, {7, 12, 13, 17}, {7, 12, 15, 17}, {7, 12, 17, 20}, {7, 12, 17, 21}, {7, 12, 17, 22}, {7, 12, 17, 24}, {7, 13, 17, 19}, {7, 15, 17, 19}, {7, 17, 19, 20}, {7, 17, 19, 21}, {7, 17, 19, 22}, {7, 17, 19, 24}, {8, 9, 12, 20}, {8, 9, 12, 24}, {8, 11, 12, 22}, {8, 11, 14, 18}, {8, 11, 14, 23}, {9, 10, 12, 15}, {9, 11, 14, 23}, {9, 11, 14, 24}, {9, 12, 14, 23}, {9, 12, 14, 24}, {9, 12, 15, 17}, {9, 12, 15, 18}, {9, 12, 15, 20}, {9, 12, 15, 23}, {9, 12, 15, 24}, {9, 14, 19, 22}, {9, 14, 19, 23}, {9, 14, 19, 24}, {11, 12, 15, 16}, {11, 13, 14, 16}, {11, 14, 16, 18}, {11, 14, 16, 20}, {11, 14, 16, 21}, {11, 14, 16, 23}, {11, 14, 16, 24}}

nops(%);

300

#Check symmetries (automorphisms):
#
GpT:= GroupTheory:
P:= CodeTools:-Usage(
    curry~(GpT:-PermApply, [GpT:-Elements(GT:-AutomorphismGroup(g))[]])
):
nops(P);

memory used=2.21MiB, alloc change=0 bytes, cpu time=187.00ms, real time=440.00ms, gc time=0ns

48

#There are 48 automorphisms; the graph has many symmetries.
#Partition the 300 dominating sets into
#equivalence classes under automorphism: 
#
EC:= CodeTools:-Usage(
    [do ec:= {(`{}`~@op@P~)(ds[1])[]} until (ds minus= ec) = {}]
);

memory used=1.12MiB, alloc change=0 bytes, cpu time=31.00ms, real time=21.00ms, gc time=0ns

[{{1, 2, 19, 22}, {1, 5, 14, 23}, {1, 13, 16, 17}, {1, 13, 22, 23}, {2, 5, 9, 24}, {2, 10, 15, 21}, {2, 10, 22, 24}, {3, 6, 8, 18}, {3, 7, 19, 22}, {3, 10, 11, 21}, {3, 18, 19, 21}, {4, 7, 8, 20}, {4, 10, 12, 15}, {4, 11, 14, 23}, {4, 14, 15, 20}, {5, 8, 18, 20}, {5, 8, 23, 24}, {6, 9, 12, 24}, {6, 9, 16, 18}, {6, 11, 13, 16}, {7, 12, 13, 17}, {7, 17, 19, 20}, {9, 12, 15, 17}, {11, 14, 16, 21}}, {{1, 6, 8, 18}, {1, 6, 19, 22}, {1, 8, 13, 17}, {1, 10, 12, 15}, {1, 10, 13, 16}, {1, 12, 14, 23}, {1, 13, 15, 23}, {1, 13, 18, 22}, {2, 4, 8, 20}, {2, 4, 19, 22}, {2, 8, 10, 21}, {2, 9, 11, 24}, {2, 10, 13, 15}, {2, 10, 16, 24}, {2, 10, 20, 22}, {2, 11, 13, 16}, {3, 5, 10, 21}, {3, 5, 14, 23}, {3, 6, 15, 18}, {3, 11, 17, 21}, {3, 12, 15, 17}, {3, 12, 18, 21}, {3, 14, 19, 22}, {3, 19, 21, 23}, {4, 6, 15, 20}, {4, 6, 16, 18}, {4, 7, 16, 20}, {4, 12, 15, 18}, {4, 14, 19, 23}, {4, 14, 20, 22}, {5, 7, 9, 24}, {5, 7, 13, 17}, {5, 8, 10, 20}, {5, 8, 13, 18}, {5, 8, 17, 24}, {5, 8, 21, 23}, {6, 9, 18, 22}, {6, 9, 19, 24}, {6, 11, 16, 20}, {7, 9, 19, 22}, {7, 11, 16, 21}, {7, 11, 17, 20}, {7, 12, 17, 21}, {7, 17, 19, 24}, {9, 11, 14, 23}, {9, 12, 14, 24}, {9, 12, 15, 23}, {11, 14, 16, 24}}, {{1, 6, 15, 18}, {1, 8, 10, 13}, {1, 12, 15, 18}, {1, 13, 15, 18}, {2, 4, 16, 20}, {2, 8, 10, 13}, {2, 10, 16, 20}, {2, 11, 16, 20}, {3, 5, 17, 21}, {3, 12, 14, 23}, {3, 12, 15, 23}, {3, 12, 21, 23}, {4, 6, 18, 22}, {4, 6, 19, 22}, {4, 6, 20, 22}, {5, 7, 17, 21}, {5, 8, 10, 13}, {5, 8, 17, 21}, {7, 9, 11, 24}, {7, 11, 16, 24}, {7, 11, 17, 24}, {9, 14, 19, 22}, {9, 14, 19, 23}, {9, 14, 19, 24}}, {{1, 6, 16, 18}, {1, 8, 13, 23}, {1, 10, 13, 22}, {1, 12, 15, 17}, {1, 13, 15, 17}, {1, 13, 16, 18}, {1, 14, 19, 22}, {1, 14, 19, 23}, {2, 4, 15, 20}, {2, 8, 10, 24}, {2, 9, 19, 22}, {2, 9, 19, 24}, {2, 10, 13, 22}, {2, 10, 15, 20}, {2, 10, 16, 21}, {2, 11, 16, 21}, {3, 5, 18, 21}, {3, 6, 18, 22}, {3, 6, 19, 22}, {3, 10, 12, 15}, {3, 10, 12, 21}, {3, 11, 14, 23}, {3, 11, 21, 23}, {3, 17, 19, 21}, {4, 6, 8, 18}, {4, 6, 8, 20}, {4, 7, 19, 22}, {4, 7, 20, 22}, {4, 12, 14, 23}, {4, 12, 15, 23}, {4, 14, 16, 20}, {5, 7, 17, 20}, {5, 8, 10, 24}, {5, 8, 13, 23}, {5, 8, 17, 20}, {5, 8, 18, 21}, {5, 9, 14, 23}, {5, 9, 14, 24}, {6, 9, 11, 24}, {6, 9, 15, 18}, {6, 11, 16, 24}, {7, 9, 12, 24}, {7, 11, 13, 16}, {7, 11, 13, 17}, {7, 12, 17, 24}, {7, 17, 19, 21}, {9, 12, 15, 18}, {11, 14, 16, 20}}, {{1, 6, 18, 22}, {1, 8, 13, 18}, {1, 10, 13, 15}, {1, 12, 15, 23}, {2, 4, 20, 22}, {2, 8, 10, 20}, {2, 10, 13, 16}, {2, 11, 16, 24}, {3, 5, 21, 23}, {3, 12, 15, 18}, {3, 12, 17, 21}, {3, 14, 19, 23}, {4, 6, 15, 18}, {4, 6, 16, 20}, {4, 14, 19, 22}, {5, 7, 17, 24}, {5, 8, 10, 21}, {5, 8, 13, 17}, {6, 9, 19, 22}, {7, 9, 19, 24}, {7, 11, 16, 20}, {7, 11, 17, 21}, {9, 11, 14, 24}, {9, 12, 14, 23}}, {{1, 7, 13, 17}, {1, 11, 13, 16}, {1, 13, 14, 23}, {1, 13, 19, 22}, {2, 3, 10, 21}, {2, 9, 10, 24}, {2, 10, 12, 15}, {2, 10, 19, 22}, {3, 6, 18, 21}, {3, 11, 16, 21}, {3, 19, 21, 22}, {4, 5, 8, 20}, {4, 7, 17, 20}, {4, 12, 15, 20}, {4, 14, 20, 23}, {5, 6, 8, 18}, {5, 8, 9, 24}, {5, 8, 14, 23}, {6, 9, 18, 24}, {6, 11, 16, 18}, {7, 12, 15, 17}, {7, 17, 19, 22}, {9, 12, 15, 24}, {11, 14, 16, 23}}, {{1, 7, 13, 23}, {1, 7, 17, 20}, {1, 11, 13, 22}, {1, 11, 16, 21}, {1, 13, 14, 17}, {1, 13, 16, 19}, {1, 14, 20, 23}, {1, 19, 21, 22}, {2, 3, 10, 24}, {2, 3, 18, 21}, {2, 9, 10, 21}, {2, 9, 18, 24}, {2, 10, 12, 22}, {2, 10, 15, 19}, {2, 12, 15, 17}, {2, 17, 19, 22}, {3, 6, 10, 21}, {3, 6, 18, 24}, {3, 11, 13, 16}, {3, 11, 21, 22}, {3, 13, 19, 22}, {3, 16, 19, 21}, {4, 5, 8, 24}, {4, 5, 15, 20}, {4, 7, 13, 17}, {4, 7, 20, 23}, {4, 8, 12, 20}, {4, 12, 15, 24}, {4, 13, 14, 23}, {4, 14, 17, 20}, {5, 6, 8, 23}, {5, 6, 16, 18}, {5, 8, 9, 20}, {5, 8, 14, 18}, {5, 9, 15, 24}, {5, 14, 16, 23}, {6, 8, 11, 18}, {6, 9, 10, 24}, {6, 9, 18, 21}, {6, 11, 16, 23}, {7, 10, 12, 15}, {7, 10, 19, 22}, {7, 12, 17, 22}, {7, 15, 17, 19}, {8, 9, 12, 24}, {8, 11, 14, 23}, {9, 12, 15, 20}, {11, 14, 16, 18}}, {{1, 7, 19, 22}, {1, 11, 14, 23}, {1, 13, 16, 23}, {1, 13, 17, 22}, {2, 3, 19, 22}, {2, 9, 12, 24}, {2, 10, 15, 24}, {2, 10, 21, 22}, {3, 6, 16, 18}, {3, 10, 19, 21}, {3, 11, 18, 21}, {4, 5, 14, 23}, {4, 7, 15, 20}, {4, 8, 14, 20}, {4, 12, 15, 17}, {5, 6, 9, 24}, {5, 8, 18, 24}, {5, 8, 20, 23}, {6, 8, 9, 18}, {6, 11, 16, 21}, {7, 12, 17, 20}, {7, 13, 17, 19}, {9, 10, 12, 15}, {11, 13, 14, 16}}, {{1, 7, 20, 23}, {1, 11, 21, 22}, {2, 3, 18, 24}, {2, 12, 17, 22}, {3, 13, 16, 19}, {4, 5, 15, 24}, {4, 13, 14, 17}, {5, 6, 16, 23}, {6, 9, 10, 21}, {7, 10, 15, 19}, {8, 9, 12, 20}, {8, 11, 14, 18}}, {{1, 9, 13, 24}, {2, 10, 14, 23}, {3, 4, 20, 21}, {5, 8, 19, 22}, {6, 7, 17, 18}, {11, 12, 15, 16}}, {{1, 9, 20, 21}, {2, 14, 17, 18}, {3, 4, 13, 24}, {5, 15, 16, 19}, {6, 7, 10, 23}, {8, 11, 12, 22}}, {{1, 14, 17, 20}, {1, 16, 19, 21}, {2, 9, 18, 21}, {2, 15, 17, 19}, {3, 6, 10, 24}, {3, 11, 13, 22}, {4, 7, 13, 23}, {4, 8, 12, 24}, {5, 9, 15, 20}, {5, 14, 16, 18}, {6, 8, 11, 23}, {7, 10, 12, 22}}]

nops~(EC); add(%);

[24, 48, 24, 48, 24, 24, 48, 24, 12, 6, 6, 12]

300

#Check global symmetries:
#    1. Given any pair [u,v] of vertices, is there an automorphism of g
#    that takes u to v?
#
evalb(nops((`{}`@op@P)~({$24})) = 1);

true

#    Yes, there is.
#

#    2. Given any pair [u,v] of vertices, is there an automorphism of g
#    that exchanges them?
#
pairs:= combinat:-choose([$24], 2):
P2:= subsindets((P~)~(pairs), listlist(posint), `[]`~@op):
andseq(pairs[k][[2,1]] in P2[k], k= 1..nops(pairs));

false

#    There is not.
#

#    3. Given any 2 pairs {u1,v1}, {u2,v2} of vertices, is there an
#    automorphism of g taking the 1st pair to the 2nd?
#
pairs:= combinat:-choose({$24}, 2):
evalb(nops(
    subsindets(subsindets((P~)~(pairs), set(list), `{}`~@op), list, `{}`@op)
) = 1);

false

#    There is not.
#

#    4. Given any two edges of g, is there an automorphism taking the 1st
#    to the 2nd?
#
edges:= GT:-Edges(g):
evalb(nops(
    subsindets(subsindets((P~)~(edges), set(list), `{}`~@op), list, `{}`@op)
) = 1);

false

#    There is not.
#

#Classify the edges by their number of automorphisms:
#
ListTools:-Classify(
    nops@`{}`@op@`{}`~@op@P~, edges
);

table( [( 12 ) = {{1, 9}, {2, 14}, {3, 4}, {5, 19}, {6, 7}, {8, 22}, {10, 23}, {11, 12}, {13, 24}, {15, 16}, {17, 18}, {20, 21}}, ( 24 ) = {{1, 2}, {1, 3}, {1, 4}, {1, 5}, {1, 7}, {1, 11}, {2, 3}, {2, 5}, {2, 6}, {2, 7}, {2, 12}, {3, 7}, {3, 8}, {3, 9}, {3, 16}, {4, 5}, {4, 9}, {4, 10}, {4, 11}, {4, 17}, {5, 6}, {5, 11}, {5, 12}, {6, 12}, {6, 13}, {6, 14}, {6, 21}, {7, 8}, {7, 14}, {7, 15}, {8, 9}, {8, 14}, {8, 15}, {8, 16}, {9, 10}, {9, 16}, {9, 17}, {10, 11}, {10, 17}, {10, 18}, {10, 19}, {11, 18}, {11, 19}, {12, 13}, {12, 19}, {12, 20}, {13, 14}, {13, 19}, {13, 20}, {13, 21}, {14, 15}, {14, 21}, {15, 21}, {15, 22}, {15, 24}, {16, 17}, {16, 22}, {16, 23}, {17, 22}, {17, 23}, {18, 19}, {18, 20}, {18, 23}, {18, 24}, {19, 20}, {20, 23}, {20, 24}, {21, 22}, {21, 24}, {22, 23}, {22, 24}, {23, 24}} ] )

nops~(%);

table( [( 12 ) = 12, ( 24 ) = 72 ] )

#So there are 12 edges with 12 automorphic images and 72 with 24 images.

 

Download DominationSet.mw

It can be done simply with subs by substituting the reciprocal of the desired denominator for the reciprocal of the current denominator, like this:

subs(1/rhs(eq_Dx)= 1/Dx, eq_K1_m4);

Let's say that you want the indices from 1 to 2, 1 to 3, 1 to 4, and 1 to 5. Then do

S:= Array(1..2, 1..3, 1..4, 1..5);

I'm guessing that @mmcdara is joking with his "non-answer", but his Reply got me thinking, and I came up with another "non-answer" that may be useful for you (so, it's actually an Answer).

If you have a device (such as a phone) with a camera and the Android or iOS operating system, then download the free Maple Calculator app to it. The app can translate photos of typeset or handwritten mathematical formulae into Maple code. I'm not aware of that photo-translation capability being available in Maple itself.

Tap the camera 📷 icon, then frame and photograph your formula. After a few seconds, you'll probably get an error message like "Sorry, this took too long to compute", which doesn't matter because you're not interested in computing it anyway. Tap the "upload to Cloud" icon. When it's done, it'll tell you the filename in the MapleCloud. In this case, the filename is the exact Maple sum command that you need.

Of course, the intended purpose of this photo-to-cloud facility is to let you download that Cloud file to Maple proper. But, for your immediate purpose, the file title has everything you need, so you don't even need to open Maple.

The error is because the target label is not in a procedure or loop. Not being in a procedure or loop is also called being "at top level".

There is no good reason for you to ever use goto. In addition to being widely frowned upon in most languages, the Maple implementation is very inefficient. 

There used to be a help page for it. But using the command is so discouraged that the help page has been removed.

First 11 12 13 14 15 16 17 Last Page 13 of 382