a bug or a feature.

[Han]

(11-12-2015 03:05 PM)ji3m Wrote:  I am having great difficulty doing cas matrix operations in a program.

The test case below shows one of the problems

EXPORT hsbat2(M)
BEGIN
LOCAL I,P,Q;
DEBUG;
I := inv(M); <- the inverse seems ok. Both M and I are lists
P := I*M; <- a problem here since '*' does silly list multiply
Q := simplify(P); <- since P is wrong this doesnt do much.
{M,I,P,Q}
END;

How are you calling the program? It runs just fine to me.
M5:=[[1,2,3],[−1,0,1],[0,0,1]];
hsbat2(M5) returns

Code:

{
  [[1,2,3],
   [−1,0,1],
   [0,0,1]],

  [[0,−1,1],
   [0.5,0.5,−2],
   [0,0,1]],

  [[1,0,0],
   [0,1,0],
   [0,0,1]],

  [[1,0,0],
   [0,1,0],
   [0,0,1]]
}

Why are you using lists? Can you show us your example that fails?