@Carl Love If one of the goals is to avoid potentially quadratic memory use (ie. O(k^2)) then it might be a good idea to explicitly unassign the names which get the indexed assignments, before the loop begins.

That would add safeguard against the user re-running the code without doing a restart, and incurring such a cost.

It could also help the user avoid a runtime error which trying to assign to a "long list", if re-run with restart, in the case that k>100.

It might also be more clear to use different names: say tableA for the indexed assignments and listA for the final conversion.

@kegj You have ignored my query: where are the statement terminators between the statements in the for-do loop?

More specifically, why is there no statement terminator at the end of this statement?

Y[j]:=Y[j-1]+dY

@maple2015 Did you try print(NCP) after first doing your plotsetup call? Try it.

@nm Your screenshots each show a single execution group containing the first three commands.

@nm I was talking about `add` in particular, because you mentioned it specifically. I have a suspicion that, apart from a general note, the only things that usually get in such Compatibility sections are mention of new or changed options.

FWIW, there are details on changes to `cat` and `rand` on the ?updates,Maple2015,Language help page.

Of course, this too does not cover your Question in full. I don't know of any uniform way to find details on all changes.

@nm Make sure that the three first commands are in separate execution groups or document blocks.

There is a long-standing problem in that the Std GUI implementation of the interface command may not take effect if it follows a restart in the same execution group.

You'll know it has worked because the printout will be in plaintext mode, not typeset (or italic).

@Preben Alsholm My ~/.mapleinit initialization file contains code to read a .mapleinit file from currentdir(), if found. Sometimes I use this mechanism to customize libname, etc, according to directory.

Basically, I often have multiple projects on the go, and instead of having to change working directory and pass the -b option to the TTY interface I can just set up the details once in a local directory-specific initialization file.

...and I have such local initialization files print something to the terminal, to remind me whether or not they are active in the current session. I don't want that output pasted in my responses here.

Sure, I could just use the -s option to prevent the initialization file from running. But I also use the TTY interface history feature quite a bit, and I usually prefer not to quit the session since I can lose the connection between that project and that command history.

There are quite a few words that I habitually misspell. Fortunately, restart doesn't seem to be one of them.

@Markiyan Hirnyk I think that it's more automatic to not have to supply v>=3, which is why I wrote that Kitonum's answer is better. I realize that English may not be not your first language, but please try and pay closer attention.

@Markiyan Hirnyk Your suggestion is not better, since it is not more automatic. It still depends on the user utilizing special knowledge of the structure; in this case the position of the value.

Kitonum's answer is more automatic, and it is better.

@richas I have some customized procedures for compacting trig expressions, which can sometimes find shorter equivalents that what simplify(...,size) gives. But they can be very computationally expensive, with quickly diminishing returns on time and memory resources invested.

One not so expensive idea is that combine can be restricted to act on just some unfrozen subset of all the trig calls in expressions like yours. The point being that simplify(...,size) doesn't attempt such mathematical conversions as combining trig calls. And combine tries to merge as much as it can, without regard to size-simplification. By using frontend to call combine with specific subsets of the trig calls unfrozen, a shorter equivalent may be found. This can be a brute force attack, leaving unfrozen all members of the full power set of the trig terms.

Another, very expensive idea is to notice that subsets of a SUM DAG (a `+` of subterms) can sometimes be simplified selectively to obtain a shorter equivalent than what simplify(...,size) gives. For example suppose the original expression T is a sum with 4 addends T[i], i=1..4. It might turn out that T[1]+T[3] and T[2]+T[4] can be size-simplified, and then the resulting pair added and size-simplied, and that this is short. Testing all ways to accomplish this kind of towered addition gets very expensive as the number or original addends increases. Even with memoization this is expensive. I would not be surprised if there were much combinatorially better ways to build up candidate equivalents.

These two ideas can be combined, as a very expensive brute force attack.

But as I mentioned above, the results from simplify(...,size) produced just as good a common-subexpression optimized procedure for computing all of [X,Y,Z] together. Well, as good as for any short [A,B,C] that I've found so far. It you have to do 35 such examples then I suggest just emitting procedures of Matlab functions is in the fast and mostly automated way in my Answer.

It might also be worth noting that it could be that neither length not MmaTranslator:-Mma:-LeafCount might be the most appropriate metric for estimating the computational cost of evaluating the resulting expression at purely floating-point values such as might occur during a simluation. How should the cost of a float evaluation of a trig call be compared with that of an addition or a multiplication? Perhaps the trig calls should be given more weight, so that reducing the total number of trig calls becomes more significant.

The following A,B,C are shorter representations for X,Y,Z respectively than what simplify(...,size) produces.

But they don't appear to produce a better computation sequence using common subexpressions (as codegen[optimize] is used for, above).

A := (-sin(q6)*L11+d1*cos(q6)+L10+((L16*sin(q10)+d2*cos(q10)+L13)*cos(q8)
 +((L16*cos(q10)-d2*sin(q10)+L15)*sin(q9)+L14*cos(q9))*sin(q8)+L12)*cos(q6+q7)
 -((-L16*cos(q10)+d2*sin(q10)-L15)*cos(q9)+L14*sin(q9))*sin(q6+q7))*cos(q5)
 -(((-L16*cos(q10)+d2*sin(q10)-L15)*sin(q9)-L14*cos(q9))*cos(q8)
 +sin(q8)*(L16*sin(q10)+d2*cos(q10)+L13))*sin(q5):

B := (-L11*sin(q6)+d1*cos(q6)+L10+((L16*sin(q10)+d2*cos(q10)+L13)*cos(q8)
 +((L16*cos(q10)-sin(q10)*d2+L15)*sin(q9)+L14*cos(q9))*sin(q8)+L12)*cos(q6+q7)
 -(L14*sin(q9)-cos(q9)*(L16*cos(q10)-sin(q10)*d2+L15))*sin(q6+q7))*sin(q5)
 +(((-L16*cos(q10)+sin(q10)*d2-L15)*sin(q9)-L14*cos(q9))*cos(q8)
 +sin(q8)*(L16*sin(q10)+d2*cos(q10)+L13))*cos(q5)-L9:

C := ((L16*cos(q10)-sin(q10)*d2+L15)*cos(q9)-L14*sin(q9))*cos(q6+q7)-L8
 +(((sin(q10)*d2-L16*cos(q10)-L15)*sin(q9)-L14*cos(q9))*sin(q8)
 +(-L16*sin(q10)-d2*cos(q10)-L13)*cos(q8)-L12)*sin(q6+q7)-d1*sin(q6)-L11*cos(q6):

length(B), length(B), length(C);
                                646, 646, 406

@Melvin Brown sorry, I'm not sure I understand about how you are using 77 on ms-windows. On Linux/Unix/OSX g77 was as front end to gcc. What's the back end C compiler on windows? Is it gcc, because you're working within cygwin?

On Linux the -s switch just strips the symbol table. But it is -shared which builds a shared object. And option -fPic to get relocatable object files. In your Question I don't see how your code would build a shared dynamic object and library. In your original what was your executable called, that you ran? How was that linked to a .dll?

@pspathis You and JurgenWolf should both try reinstalling, possibly using a fresh download of the full Installer. And if you don't know how send a message to support@maplesoft.com

ps. pspathis, please don't spam by reposting your same message as comments on other topic/threads on this site.

@mzaman You didn't change the indexing from round brackets to square brackets.

FWIW I am seeing the following results in Maple 2015.0 on 64bit Linux,

int( tan(x)^(n-2)*sec(x)^2, x ); 

                                       (n - 1)
                                 tan(x)
                                 -------------
                                     n - 1

int( -csc(x)^2/cot(x)^(2/3), x );

                                         (1/3)
                                 3 cot(x)

acer