Simple solution

Preben Alsholm — Thu, 28 Feb 2013 07:27:28 Z

Since Ts and Ths are lists of matrices and Q is a matrix, you could do
Qb:=[seq(Ts[i]^(-1).Q.Ths[i],i=1..nops(Ts))];
to get a list of the matrices you need.

Another way:

Qb:=Ts^~(-1).~[Q$nops(Ts)].~Ths;

I find the first one easier to understand.

radians and zip

Carl Love — Thu, 28 Feb 2013 08:23:25 Z

Before answering your main question, I want to point out that angles need to be in radians in Maple. You can enter them in degrees and then convert to radians, like this

Angles:= [30, 60, 90];
Angles:= Angles *~ Pi/180:

Angles:= map(convert, Angles *~ degrees, radians):

Second, in my previous post, I said that you'd probably want to use more than 4 digits in actual practice, but you kept it at 4. You should change all occurences of evalf[4] to evalf[15] or simply evalf (for reasons that I won't go into here, 15 digits usually gives the optimal balance of speed and accuracy).

Third, you should probably use different names for the two Constructor procedures: the one that constructs the Ts and the one that constructs the Ths.

Now for your main questions. You have two lists of eight matrices each, Ts and Ths, and a single matrix Q, and you want to compute a list of matrices Qb such that Qb[k] = Ts[k]^(-1).Q.Ths[k] for k from 1 to 8. Here are two ways do that:

Qb:= [seq](Ts[k]^(-1).Q.Ths[k], k= 1..8):

Qb:= zip((X,Y)-> X^(-1).Q.Y, Ts, Ths):

Finally, you want to construct a 6x6 matrix from 3x3 matrices A, B, and D, with pattern

A B
B D.

The command for this is

M:= < < A | B > , < B | D > >:

Let me know if that is all clear to you!

MaplePrimes - answers and comments on Question, Automating matrix calculations

Simple solution

radians and zip