## 490 Reputation

7 years, 115 days
i would also not like to ask, but if not ask, what should i do?

## @Kitonum  i edited your function f...

i edited your function for all elements in matrix, but not all cases considered, such as transpose or partly transpose

with(LinearAlgebra);
IsEquivalent:=proc(obj1::Matrix, obj2::Matrix, var::list)
local P, S;
uses combinat, LinearAlgebra;
P:=permute(var);
#print(P):
S:=[seq(subs(zip(`=`,var,p),obj2), p=P)];
#print(S):
SS := choose(S,2);
result := false:
for i from 1 to nops(SS) do
if `or`(seq(Equal(obj1,s), s=SS[i])) then
result := true:
end if:
od:
return result:
end proc:
elementss := proc(AA)
local elements1:
elements1 := [op(indets(AA))];
elements1 := [op(elements1), op(indets(AA,float))];
elements1 := [op(elements1),op(indets(AA,constant))];
return elements1:
end proc:
A1 := Matrix([[777,x,x],[77,7,7],[2,y,y]]);
B1 := Matrix([[777,y,y],[77,7,7],[2,x,x]]);
IsEquivalent(A1, B1, elementss(A1)); # correct
A1 := Matrix([[777,x,x],[77,7,7],[2,y,y]]);
B1 := Matrix([[777,y,y],[77,7,7],[2,x,y]]);
IsEquivalent(A1, B1, elementss(A1)); # correct

A1 := Matrix([[x,x,x],[77,7,7],[z,y,y]]);
B1 := Matrix([[z,y,y],[77,7,7],[x,x,x]]);
IsEquivalent(A1, B1, elementss(A1)); # incorrect

A1 := Matrix([[x,x,x],[77,7,7],[z,y,y]]);
B1 := Matrix([[z,y,y],[77,7,7],[x,x,y]]);
IsEquivalent(A1, B1, elementss(A1)); # correct

A1 := Matrix([[777,x,x],[77,7,7],[2,y,y]]);
B1 := Transpose(A1);
IsEquivalent(A1, B1, elementss(A1)); # incorrect

A1 := Matrix([[777,x,x],[77,7,7],[2,y,y]]);
B1 := Matrix([[777,x,2],[77,7,7],[x,y,y]]);
IsEquivalent(A1, B1, elementss(A1)); # incorrect

A1 := Matrix([[777,x,x],[77,7,7],[2,y,y]]);
B1 := Matrix([[777,77,2],[x,7,7],[x,y,y]]);
IsEquivalent(A1, B1, elementss(A1)); # incorrect

## @Kitonum  i got error when input a...

i got error when input a list having 3 elements x, y, z

it expects 2 arguments

how to edit the function for more than 2 elements?

## @Kitonum  is there any clues to wr...

is there any clues to write the extension?

what things need to consider?

## @Kitonum    in maple 12 , ...

in maple 12 , it return not true, can it change to run in maple 12?

second question is about two operations, such as symmetric flip and then rotate

for example:

IsEquivalent(Matrix([[777,x,x],[77,7,7],[2,y,y]]), Matrix([[777,y,y],[77,7,x],[2,7,x]]), [x,y]);

## @Kitonum can this express in print ...

can this express in print or printf to show the operation it recognised ?

because expect to count how many same operations in series.

I have some understanding not clear is that

if there are two or three operations to achieve between two matrix , is it one equivalent can be considered as one operations for two or more operations?

## @Joe Riel  i mean self modifying m...

i mean self modifying maple code application

and need a database to record the change

it write maple code itself

## @Kitonum  Not use procedure  ...

Not use procedure

and run directly , it can show animation now

## @Kitonum do not know why i run in m...

do not know why i run in maple 2015 , no result.

## @Kitonum  i run your code in maple...

i run your code in maple 12 and maple 2015 , no graph and animation was shown.

## @Kitonum assume rotation matrix is&...

assume rotation matrix is (Matrix([[cos(t),sin(t)],[sin(t),cos(t)]])

how to animate a particle rotate around another particle?

how to plot the track of x and y coordinate with the change of t from -Pi to Pi?

would like to draw this track in 2D diagram but there are 3 variable?

how to set an initial point for x and y and then change t to see how x and y change in 2D coordinate?

## @vv  i understand now, by the way...

i understand now,

by the way,

do you know how to predict a future matrix

belong to which one of 3 matrix?

actually eigenvector has very small difference in

Decimal value , can linear regression help to predict the exact eigenvector values?

i find book mention magnetron , neucleon these words, but I use Time series only to find Pauli equation. Is it the end after find Pauli equation ?

how to use Pauli equation for Time series?

## @asa12  i succeed to do and comple...

i succeed to do and completed today

## @Kitonum  @Kitonum    ...

with(LinearAlgebra):
A := Matrix(2, 2, {(1, 1) = (1/6)*sqrt(3)+(1/2)*I, (1, 2) = (1/6)*sqrt(3)-(1/2)*I, (2, 1) = (1/6)*sqrt(3)-(1/2)*I, (2, 2) = (1/6)*sqrt(3)+(1/2)*I}):
with(MTM):
factor(Re(A));
A3 := Matrix([[0,0],[0,0]]):
for i from 1 to 2 do
for j from 1 to 2 do
A3[i,j] := Im(A[i,j]);
od:
od:
counterd := 0:
for i from 1 to 2 do
for j from 1 to 2 do
if A2[i,j] mod ((1/6)*sqrt(3)) = True then
counterd := counterd + 1:
end if:
od:
od:
counterd;
if counterd = 4 then
for i from 1 to 2 do
for j from 1 to 2 do
A2[i,j] := A2[i,j]/((1/6)*sqrt(3)):
od:
od:
end if:
A2;
A2 := Matrix([[0,0],[0,0]]):
for i from 1 to 2 do
for j from 1 to 2 do
A2[i,j] := Im(A[i,j]);
od:
od:
counterd := 0:
for i from 1 to 2 do
for j from 1 to 2 do
if A2[i,j] mod (1/2) = True then
counterd := counterd + 1:
end if:
od:
od:
counterd;
if counterd = 4 then
for i from 1 to 2 do
for j from 1 to 2 do
A2[i,j] := A2[i,j]/(1/2):
od:
od:
end if:
Error, invalid input: modp received 1/2, which is not valid for its 2nd argument, m;

i do it in another method
mod has error in maple 12

how to mod in maple 12 for fraction ?

## @Kitonum  moreover , is it possibl...

moreover , is it possible sqrt(3)/6 and 1/2 combine into one term and the rest is one matrix result?

## @Kitonum A := Matrix(2, 2, {(1, 1) ...

Your doing is correct.

A := Matrix(2, 2, {(1, 1) = (1/6)*sqrt(3)+(1/2)*I, (1, 2) = (1/6)*sqrt(3)-(1/2)*I, (2, 1) = (1/6)*sqrt(3)-(1/2)*I, (2, 2) = (1/6)*sqrt(3)+(1/2)*I}):
A1:=(sqrt(3)/6)%*(Matrix(2,(i,j)->Re(A[i,j]))/(sqrt(3)/6))%+((1/2)%*(Matrix(2,(i,j)->A[i,j]-Re(A[i,j]))/(1/2)));

Error, (in rtable/Product) invalid arguments
in maple 12

how to run in maple 12?

 5 6 7 8 9 10 11 Last Page 7 of 44
﻿