RukeとLuNaYuの日記
I know the truth.
I know whole.
And I...know you.
平凡な大学生活の日記です。時折まじめな長文を書く病気になります。興味がなければ読み飛ばしてください。
201708 << 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 >> 201710
スポンサーサイト (--/--/--(--) --:--:--)
上記の広告は1ヶ月以上更新のないブログに表示されています。
新しい記事を書く事で広告が消せます。
オクターヴ (2006/02/07(火) 02:34:13)
一応ここまでの計算のoctaveスクリプトを載せておく。
% prepare generators
p1=zeros(60);
% (56 57 58 59 60)
p1(56,60)=1;
p1(57,56)=1;
p1(58,57)=1;
p1(59,58)=1;
p1(60,59)=1;
% (42 45 48 51 54)
p1(42,54)=1;
p1(45,42)=1;
p1(48,45)=1;
p1(51,48)=1;
p1(54,51)=1;
% (41 44 47 50 53)
p1(41,53)=1;
p1(44,41)=1;
p1(47,44)=1;
p1(50,47)=1;
p1(53,50)=1;
% (43 46 49 52 55)
p1(43,55)=1;
p1(46,43)=1;
p1(49,46)=1;
p1(52,49)=1;
p1(55,52)=1;
% (22 26 30 34 38)
p1(22,38)=1;
p1(26,22)=1;
p1(30,26)=1;
p1(34,30)=1;
p1(38,34)=1;
% (23 27 31 35 39)
p1(23,39)=1;
p1(27,23)=1;
p1(31,27)=1;
p1(35,31)=1;
p1(39,35)=1;
% (21 25 29 33 37)
p1(21,37)=1;
p1(25,21)=1;
p1(29,25)=1;
p1(33,29)=1;
p1(37,33)=1;
% (24 28 32 36 40)
p1(24,40)=1;
p1(28,24)=1;
p1(32,28)=1;
p1(36,32)=1;
p1(40,36)=1;
% (7,10,13,16,19);
p1(7,19)=1;
p1(10,7)=1;
p1(13,10)=1;
p1(16,13)=1;
p1(19,16)=1;
% (8 11 14 17 20)
p1(8,20)=1;
p1(11,8)=1;
p1(14,11)=1;
p1(17,14)=1;
p1(20,17)=1;
% (6 9 12 15 18)
p1(6,18)=1;
p1(9,6)=1;
p1(12,9)=1;
p1(15,12)=1;
p1(18,15)=1;
% (1 2 3 4 5)
p1(1,5)=1;
p1(2,1)=1;
p1(3,2)=1;
p1(4,3)=1;
p1(5,4)=1;

p2=zeros(60);
% (30 31 49 48 47)
p2(30,47)=1;
p2(31,30)=1;
p2(49,31)=1;
p2(48,49)=1;
p2(47,48)=1;
% (29 32 50 58 46)
p2(29,46)=1;
p2(32,29)=1;
p2(50,32)=1;
p2(58,50)=1;
p2(46,58)=1;
% (13 33 51 57 27)
p2(13,27)=1;
p2(33,13)=1;
p2(51,33)=1;
p2(57,51)=1;
p2(27,57)=1;
% (14 34 59 45 28)
p2(14,28)=1;
p2(34,14)=1;
p2(59,34)=1;
p2(45,59)=1;
p2(28,45)=1;
% (11 15 35 60 44)
p2(11,44)=1;
p2(15,11)=1;
p2(35,15)=1;
p2(60,35)=1;
p2(44,60)=1;
% (12 16 52 56 26)
p2(12,26)=1;
p2(16,12)=1;
p2(52,16)=1;
p2(56,52)=1;
p2(26,56)=1;
% (3 17 53 42 25)
p2(3,25)=1;
p2(17,3)=1;
p2(53,17)=1;
p2(42,53)=1;
p2(25,42)=1;
% (4 36 54 43 10)
p2(4,10)=1;
p2(36,4)=1;
p2(54,36)=1;
p2(43,54)=1;
p2(10,43)=1;
% (2 18 38 41 24)
p2(2,24)=1;
p2(18,2)=1;
p2(38,18)=1;
p2(41,38)=1;
p2(24,41)=1;
% (5 37 55 23 9)
p2(5,9)=1;
p2(37,5)=1;
p2(55,37)=1;
p2(23,55)=1;
p2(9,23)=1;
% (1 19 39 22 8)
p2(1,8)=1;
p2(19,1)=1;
p2(39,19)=1;
p2(22,39)=1;
p2(8,22)=1;
% (6 20 40 21 7)
p2(6,7)=1;
p2(20,6)=1;
p2(40,20)=1;
p2(21,40)=1;
p2(7,21)=1;

s=zeros(60);
% (58)
s(58,58)=1;
% (57,59)
s(57,59)=1;
s(59,57)=1;
% (56,60)
s(56,60)=1;
s(60,56)=1;
% (48)
s(48,48)=1;
% (47,49)
s(47,49)=1;
s(49,47)=1;
% (46,50)
s(46,50)=1;
s(50,46)=1;
% (45,51)
s(45,51)=1;
s(51,45)=1;
% (44,52)
s(44,52)=1;
s(52,44)=1;
% (43,53)
s(43,53)=1;
s(53,43)=1;
% (42,54)
s(42,54)=1;
s(54,42)=1;
% (41,55)
s(41,55)=1;
s(55,41)=1;
% (30,31)
s(30,31)=1;
s(31,30)=1;
% (29,32)
s(29,32)=1;
s(32,29)=1;
% (28,33)
s(28,33)=1;
s(33,28)=1;
% (27,34)
s(27,34)=1;
s(34,27)=1;
% (26,35)
s(26,35)=1;
s(35,26)=1;
% (25,36)
s(25,36)=1;
s(36,25)=1;
% (24,37)
s(24,37)=1;
s(37,24)=1;
% (23,38)
s(23,38)=1;
s(38,23)=1;
% (22,39)
s(22,39)=1;
s(39,22)=1;
% (21,40)
s(21,40)=1;
s(40,21)=1;
% (13,14)
s(13,14)=1;
s(14,13)=1;
% (12,15)
s(12,15)=1;
s(15,12)=1;
% (11,16)
s(11,16)=1;
s(16,11)=1;
% (10,17)
s(10,17)=1;
s(17,10)=1;
% (9,18)
s(9,18)=1;
s(18,9)=1;
% (8,19)
s(8,19)=1;
s(19,8)=1;
% (7,20)
s(7,20)=1;
s(20,7)=1;
% (6)
s(6,6)=1;
% (3,4)
s(3,4)=1;
s(4,3)=1;
% (2,5)
s(2,5)=1;
s(5,2)=1;
% (1)
s(1,1)=1;

character_1=trace(ones(60))

character_12C5=trace(p1)

character_12C5_2=trace(p1^2)

character_20C3=trace(p1*p2)

character_15C2=trace(p2*p1*p2)

i=s*p1*p2^4*p1*p2^3*p1^4;
character_i=trace(i)

character_12S10=trace(i*p1)

character_12S10_3=trace(i*p1^2)

character_20S6=trace(i*p1*p2)

character_15sigma=trace(i*p1*p2*p1)

sum_1=eye(60);

sum_12C5=p1\
+p2+p1^1*(p2)*p1^4+p1^2*(p2)*p1^8 +p1^3*(p2)*p1^12+p1^4*(p2)*p1^16\
+p1^4\
+p2^4+p1^1*(p2^4)*p1^4+p1^2*(p2^4)*p1^8 +p1^3*(p2^4)*p1^12+p1^4*(p2^4)*p1^16;

sum_12C5_2=p1^2\
+p2^2+p1^1*(p2^2)*p1^4+p1^2*(p2^2)*p1^8 +p1^3*(p2^2)*p1^12+p1^4*(p2^2)*p1^16\
+p1^8\
+p2^8+p1^1*(p2^8)*p1^4+p1^2*(p2^8)*p1^8 +p1^3*(p2^8)*p1^12+p1^4*(p2^8)*p1^16;

sum_20C3=p1*p2+p1^1*(p1*p2)*p1^4+p1^2*(p1*p2)*p1^8 +p1^3*(p1*p2)*p1^12+p1^4*(p1*p2)*p1^16\
+p2^4*p1^4+p1^1*(p2^4*p1^4)*p1^4+p1^2*(p2^4*p1^4)*p1^8 +p1^3*(p2^4*p1^4)*p1^12+p1^4*(p2^4*p1^4)*p1^16\
+(p2^4*p1*p2^2)+p1^1*(p2^4*p1*p2^2)*p1^4 +p1^2*(p2^4*p1*p2^2)*p1^8 +p1^3*(p2^4*p1*p2^2)*p1^12+p1^4*(p2^4*p1*p2^2)*p1^16\
+(p2^8*p1^4*p2^1)+p1^1*(p2^8*p1^4*p2^1)*p1^4 +p1^2*(p2^8*p1^4*p2^1)*p1^8+p1^3*(p2^8*p1^4*p2^1)*p1^12 +p1^4*(p2^8*p1^4*p2^1)*p1^16;

sum_15C2=(p2*p1*p2)+p1*(p2*p1*p2)*p1^4 +p1^2*(p2*p1*p2)*p1^8 +p1^3*(p2*p1*p2)*p1^12 +p1^4*(p2*p1*p2)*p1^16\
+(p1*p2^2)+p1*(p1*p2^2)*p1^4+p1^2*(p1*p2^2)*p1^8 +p1^3*(p1*p2^2)*p1^12+p1^4*(p1*p2^2)*p1^16\
+(p2^4*p1*p2^3)+p1*(p2^4*p1*p2^3)*p1^4+p1^2*(p2^4*p1*p2^3)*p1^8 +p1^3*(p2^4*p1*p2^3)*p1^12+p1^4*(p2^4*p1*p2^3)*p1^16;

sum_i=i*sum_1;
sum_12S10=i*sum_12C5;
sum_12S10_3=i*sum_12C5_2;
sum_20S6=i*sum_20C3;
sum_15sigma=i*sum_15C2;

source "hamiltonian_octave.m"

P_Ag=(1/120)*(1*sum_1+1*sum_12C5+1*sum_12C5_2 +1*sum_20C3+1*sum_15C2 +1*sum_i+1*sum_12S10 +1*sum_12S10_3+1*sum_20S6 +1*sum_15sigma);
P_T1g=(3/120)*(3*sum_1+(1+sqrt(5))/2*sum_12C5 +(1-sqrt(5))/2*sum_12C5_2+0*sum_20C3 -1*sum_15C2+3*sum_i+(1-sqrt(5))/2*sum_12S10 +(1+sqrt(5))/2*sum_12S10_3 +0*sum_20S6-1*sum_15sigma);
P_T2g=(3/120)*(3*sum_1+(1-sqrt(5))/2*sum_12C5 +(1+sqrt(5))/2*sum_12C5_2+0*sum_20C3 -1*sum_15C2+3*sum_i+(1+sqrt(5))/2*sum_12S10 +(1-sqrt(5))/2*sum_12S10_3 +0*sum_20S6-1*sum_15sigma);

state=[rand(1,60)]';
(H*P_Ag*state)./(P_Ag*state)
(H*P_T1g*state)./(P_T1g*state)
(H*P_T2g*state)./(P_T2g*state)
スポンサーサイト
コメント
この記事へのコメント
コメントを投稿する

管理者にだけ表示を許可する
トラックバック
この記事のトラックバックURL
この記事へのトラックバック
(C)Copyright 2003-2007 by Ruke All rights reserved. Powered By FC2. VALID HTML? VALID CSS?
上記広告は1ヶ月以上更新のないブログに表示されています。新しい記事を書くことで広告を消せます。