ActiveUser

475 Reputation

10 Badges

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

MaplePrimes Activity


These are replies submitted by ActiveUser

@tomleslie 

very good function

it seems not need accumulate function

@vv 

my notebook is faster to run web now.

i attached my code,

but the result different from table of the book photo

at orbit size 3 , only two number of each, but book photo show all are 3 numbers

restart:
with(combinat):
with(group):
FixPoints := (a,n) -> {seq(i,i=1..n)} minus {op(map(op,a))}:
a := [[seq(i,i=1..6)],[seq(i,i=7..9)]]:
b := a:
aa := []:
countlist := []:
for i to 6 do
tt := b;
if tt = [] then
 tt := permgroup(9, {[]}):
end if:
count := 1:
orbitlist := []:
continuerun := 0:
for j from 1 to 7 do
if continuerun = 0 then
if b = [] then
if orbit(tt,j) in orbitlist then
 continuerun := 1:
else
 orbitlist := [op(orbitlist), orbit(tt,j)]:
end if:
else
if orbit({tt},j) in orbitlist then
 continuerun := 1:
else
 orbitlist := [op(orbitlist), orbit({tt},j)]:
end if:
end if:
end if:
od:
if nops(orbitlist) in countlist then
 rr := 0:
else
 countlist := [op(countlist), nops(orbitlist)]:
 lprint('a'^i=b, "orbit size=",nops(orbitlist)):
end if:
aa := [op(aa), b]:
b := mulperms(b,a);
od:

 

@vv 

I have network security concern to not to use mobile carrier network

how about the next two questions?

start from a^2 , result different from book photo?

how to calculate the orbit result at orbit size 3 at bottom table of the book photo?

@vv 

sorry for missing one photo

my notebook load maple prime web very slow, I can only use iphone

start from a^2 , result different from book above photo

and aa0 can not be calculated with orbit , and how to calculate the orbit result at orbit size 3 at bottom table of the book photo

 

@vv 

i understand now

a^4 repeated for a^2

so not choose a^4

and if calculate larger than orbit size 

it will repeated too

but I do not understand , start from a^2

the result is different from the book

and aa0 , can not calculate with orbit

return error

moreover when aa[3]= a^3

but how calculate the result having 3 numbers 036, 147 etc at orbit size 3

@vv 

why orbit size 1 is 678

how to calculate this?

a^6 = a^0

how to calculate (018) .. etc at the bottom table

i use fixpoints and orbit function not work

 

and why it choose orbit size 1,2,3,6 only?

@Kitonum 

but i need a while to do a further combination 

from this result I need to choose 51 different rows

and search

@tomleslie 

I need 1 and 0 in all combination as described above

i expect larger , 10 one and 75 zero

@Mariusz Iwaniuk 

my application is quantum matrix

one of invariant value sometimes is not exactly 1.0 

i am thinking why invariant is not always invariant.

i searched that is caused by superposition 

but convolution of list and sum result is the same as distributive law

so, how to use convolution to make invariant always invariant?

@Mariusz Iwaniuk 

 

restart:
conv := proc(a, b)
c := []:
for k from 1 to nops(a)+nops(b)-1 do
 result := 0:
 for i from 1 to k do
  if k-i+1 <= nops(b) and i <= nops(a) then
  result := result + a[i]*b[k-i+1]:
  end if:
 od:
 c := [op(c), result]:
od:
return c:
end proc:
with(LinearAlgebra):
M := Matrix([[1,2,3],[4,5,6],[7,8,9]]):
X := Matrix([[x1],[x2],[x3]]);
MatrixMatrixMultiply(M,X);
something should be convsolve(conv(M, X) = a value);
row1 := conv([seq(M[1][a],a=1..nops(M[1]))], [seq(Transpose(X)[1][a],a=1..nops(Transpose(X)))]):
sum(row1[mm], mm=1..nops(row1));
row2 := conv([seq(M[2][a],a=1..nops(M[2]))], [seq(Transpose(X)[1][a],a=1..nops(Transpose(X)))]);
sum(row2[mm], mm=1..nops(row2));
row3 := conv([seq(M[3][a],a=1..nops(M[3]))], [seq(Transpose(X)[1][a],a=1..nops(Transpose(X)))]);
sum(row3[mm], mm=1..nops(row3));
solve([sum(row1[mm], mm=1..nops(row1)), sum(row2[mm], mm=1..nops(row2)), sum(row3[mm], mm=1..nops(row3))], [x1,x2,x3]);
 
 
i do not know the definition of convolution of matrix with x to form a system of equations
i guessed shown above
if above is correct, then does it mean that solve can also solve this system?
 
how to replace distributive rules with convolution in maple?
if this can not be done, can i do it when only not use bracket ( ) to do multiply?
 
is there any example to show only convolution is correct when distributive is wrong?
 

@tomleslie 

is there any one know maple 12 use which method to expand the size of display graph?

@tomleslie 

 

what is yr ?

is there an simple example using dualaxisplot 

vertical line is only x=1, x=2, x=3 etc for example around 12 lines 

because i collect these lines with

L := [op(L), 1]

L := [op(L), 2]

etc

@acer 

firstly I enabled maple security, then =0 not work then I disable security, it work but can not changed yesterday

after tried today, it work in window maple

i just curious why this command is applied globally but not session applied and it has affected after closed and open again

then I try again window discover =2 , display %1 = imaginary value before display matrix

then set =3 to solve this issue in maple 

I use =0 mainly for copying many matrix in window 

because =2 or =3 can not copy matrix

@stefanv 

but i really find invariant in 4x4

why is it pointless in more dimension?

there is no papers mention how to construct a general invariant matrix

and there is not only one kind of matrix invariant. Is there book teaching different kind of general module to find numerical invariant matrix?

@mmcdara 

i do not know the formula 

is there example code to show what happen?

If cube is not matrix

why do people solve cube game?

 

1 2 3 4 5 6 7 Last Page 3 of 42