[求助]这个对角化程序对不对啊

a1=0;
a2=0;
a3=0;
b1=500;
b2=500;
b3=500;
b4=500;
c1=2000;
c2=2000;
c3=2000;
e1=3600;
e2=3600;
e3=5400;
e4=5400;
f1=2000;
f2=200;


B=zeros(45);

B(1,1)=a2-8*b2;
B(2,2)=a2-8*b2;
B(3,3)=a2-8*b2;

B(4,4)=a2+8*b2+4*c2;B(4,7)=2*6^(1/2)*b3+6^(1/2)*c3;

B(5,5)=a2+2*c1;B(5,8)=-2*3^(1/2)*b3;
B(6,6)=a2+2*c1;B(6,9)=2*3^(1/2)*b3;
B(7,7)=a1+10*b1+5*c1;
B(8,8)=a1+b1+2*c1;
B(9,9)=a1+b1+2*c1;

B(10,10)=a1-5*b1;B(10,22)=6*b3;
B(11,11)=a1-5*b1;B(11,23)=6*b3;
B(12,12)=a1-5*b1;B(12,24)=6*b3;
B(13,13)=a1-5*b1;B(13,25)=6*b3;
B(14,14)=a1-5*b1;B(14,26)=6*b3;
B(15,15)=a1-5*b1;B(15,27)=6*b3;
B(16,16)=a1-5*b1;B(16,28)=6*b3;
B(17,17)=a1-5*b1;B(17,29)=6*b3;
B(18,18)=a1-5*b1;B(18,30)=6*b3;
B(19,19)=a1+b1+2*c1;B(19,43)=2*3^(1/2)*b1;
B(20,20)=a1+b1+2*c1;B(20,44)=2*3^(1/2)*b1;
B(21,21)=a1+b1+2*c1;B(21,45)=2*3^(1/2)*b1;
B(22,22)=a3+4*b3;
B(23,23)=a3+4*b3;
B(24,24)=a3+4*b3;
B(25,25)=a3+4*b3;
B(26,26)=a3+4*b3;
B(27,27)=a3+4*b3;
B(28,28)=a3+4*b3;
B(29,29)=a3+4*b3;
B(30,30)=a3+4*b3;
B(31,31)=a3+4*b3+2*c3;
B(32,32)=a3+4*b3+2*c3;
B(33,33)=a3+4*b3+2*c3;
B(34,34)=a3-8*b3;
B(35,35)=a3-8*b3;
B(36,36)=a3-8*b3;
B(37,37)=a3-8*b3;
B(38,38)=a3-8*b3;
B(39,39)=a3-8*b3;
B(40,40)=a3-8*b3;
B(41,41)=a3-8*b3;
B(42,42)=a3-8*b3;
B(44,44)=a3+2*c3;
B(43,43)=a3+2*c3;
B(45,45)=a3+2*c3;

C=zeros(45);

C(1,7)=-1/3*3^(1/2)*i*(-f2);
C(1,8)=-1/3*3^(1/2)*i*(f2);
C(1,9)=-1/3*3^(1/2)*i*(f2);
C(1,10)=(-f2);
C(1,12)=(-f2);
C(1,13)=-1/2*2^(1/2)*(f2);
C(1,15)=-1/2*2^(1/2)*(f2);
C(1,19)=-1/6*6^(1/2)*i*(-f2);
C(1,20)=-1/6*6^(1/2)*i*(2*f2);
C(1,21)=-1/6*6^(1/2)*i*(f2);
C(1,23)=-1/2*2^(1/2)*(f2);
C(1,24)=-1/2*2^(1/2)*(f2);
C(1,26)=-1/2*(f2);
C(1,27)=-1/2*(f2);
C(1,32)=-1/2*(-f2);
C(1,33)=-1/2*(-f2);
C(1,35)=-1/2*2^(1/2)*i*(f2);
C(1,36)=1/2*2^(1/2)*i*(f2);
C(1,38)=-1/2*i*(f2);
C(1,39)=1/2*i*(f2);
C(1,44)=1/2*i*(f2);
C(1,45)=-1/2*i*(f2);

C(2,7)=-1/6*6^(1/2)*i*(-f2);
C(2,8)=-1/6*6^(1/2)*i*(f2);
C(2,9)=-1/6*6^(1/2)*i*(f2);
C(2,10)=-1/2*2^(1/2)*(-f2);
C(2,12)=-1/2*2^(1/2)*(f2);
C(2,13)=-1/2*(f2);
C(2,15)=-1/2*(f2);
C(2,19)=-1/12*12^(1/2)*i*(-f2);
C(2,20)=-1/12*12^(1/2)*i*(2*f2);
C(2,21)=-1/12*12^(1/2)*i*(f2);
C(2,23)=-1/2*(f2);
C(2,24)=-1/2*(f2);
C(2,26)=-1/2*4^(1/2)*(f2);
C(2,27)=-1/2*4^(1/2)*(f2);
C(2,32)=-1/2*4^(1/2)*(-f2);
C(2,33)=1/2*4^(1/2)*(f2);
C(2,35)=-1/2*i*(f2);
C(2,36)=1/2*i*(f2);
C(2,38)=-1/4*2^(1/2)*i*(f2);
C(2,39)=1/4*2^(1/2)*i*(f2);
C(2,44)=-1/4*2^(1/2)*i*(-f2);
C(2,45)=-1/4*2^(1/2)*i*(f2);

C(4,7)=-1/6*6^(1/2)*i*(-i)*(-f2);
C(4,8)=-1/6*6^(1/2)*i*(-i)*(f2);
C(4,9)=-1/6*6^(1/2)*i*(-i)*(f2);
C(4,10)=-1/2*2^(1/2)*(-i)*(-f2);
C(4,12)=-1/2*2^(1/2)*(-i)*(f2);
C(4,13)=-1/2*(-i)*(f2);
C(4,15)=-1/2*(-i)*(f2);
C(4,19)=-1/12*12^(1/2)*i*(-i)*(-f2);
C(4,20)=-1/12*12^(1/2)*i*(-i)*(2*f2);
C(4,21)=-1/12*12^(1/2)*i*(-i)*(f2);
C(4,23)=-1/2*(-i)*(f2);
C(4,24)=-1/2*(-i)*(f2);
C(4,26)=-1/2*4^(1/2)*(-i)*(f2);
C(4,27)=-1/2*4^(1/2)*(-i)*(f2);
C(4,32)=-1/2*4^(1/2)*(-i)*(-f2);
C(4,33)=1/2*4^(1/2)*(-i)*(f2);
C(4,35)=-1/2*i*(-i)*(f2);
C(4,36)=1/2*i*(-i)*(f2);
C(4,38)=-1/4*2^(1/2)*i*(-i)*(f2);
C(4,39)=1/4*2^(1/2)*i*(-i)*(f2);
C(4,44)=-1/4*2^(1/2)*i*(-i)*(-f2);
C(4,45)=-1/4*2^(1/2)*i*(-i)*(f2);

C(6,7)=1/3*3^(1/2)*i*(i)*(-f2);
C(6,8)=1/3*3^(1/2)*i*(i)*(f2);
C(6,9)=1/3*3^(1/2)*i*(i)*(f2);
C(6,10)=-(i)*(-f2);
C(6,12)=-(i)*(-f2);
C(6,13)=1/2*2^(1/2)*(i)*(f2);
C(6,15)=1/2*2^(1/2)*(i)*(f2);
C(6,19)=1/6*6^(1/2)*i*(i)*(-f2);
C(6,20)=1/6*6^(1/2)*i*(i)*(2*f2);
C(6,21)=1/6*6^(1/2)*i*(i)*(f2);
C(6,23)=1/2*2^(1/2)*(i)*(f2);
C(6,24)=1/2*2^(1/2)*(i)*(f2);
C(6,26)=1/2*(i)*(f2);
C(6,27)=1/2*(i)*(f2);
C(6,32)=1/2*(i)*(-f2);
C(6,33)=1/2*(i)*(-f2);
C(6,35)=1/2*2^(1/2)*i*(i)*(f2);
C(6,36)=-1/2*2^(1/2)*i*(i)*(f2);
C(6,38)=1/2*i*(i)*(f2);
C(6,39)=-1/2*i*(i)*(f2);
C(6,44)=-1/2*i*(i)*(f2);
C(6,45)=1/2*i*(i)*(f2);

C(7,7)=1/3*(f1);
C(7,8)=1/3*(-f1);
C(7,9)=1/3*(-f1);
C(7,10)=-1/3*3^(1/2)*i*(-f1);
C(7,12)=-1/3*3^(1/2)*i*(-f1);
C(7,13)=-1/6*6^(1/2)*i*(-f1);
C(7,15)=-1/6*6^(1/2)*i*(-f1);
C(7,19)=1/18*18^(1/2)*(f1);
C(7,20)=1/18*18^(1/2)*(-2*f1);
C(7,21)=1/18*18^(1/2)*(-f1);
C(7,22)=-1/6*6^(1/2)*i*(-f2);
C(7,23)=-1/6*6^(1/2)*i*(-f2-f1);
C(7,24)=1/6*6^(1/2)*i*(-f2+f1);
C(7,25)=-1/6*3^(1/2)*i*(-f2);
C(7,26)=-1/6*3^(1/2)*i*(-f2-f1);
C(7,27)=1/6*3^(1/2)*i*(-f2+f1);
C(7,31)=-1/6*3^(1/2)*i*(-f2);
C(7,32)=-1/6*3^(1/2)*i*(-f2+f1);
C(7,33)=1/6*3^(1/2)*i*(-f2-f1);
C(7,34)=1/6*6^(1/2)*(-f2);
C(7,35)=1/6*6^(1/2)*(-f2-f1);
C(7,36)=1/6*6^(1/2)*(-f2-f1);
C(7,37)=1/6*3^(1/2)*(-f2);
C(7,38)=1/6*3^(1/2)*(-f2-f1);
C(7,39)=-1/6*3^(1/2)*(-f2-f1);
C(7,43)=1/6*3^(1/2)*(-f2);
C(7,44)=1/6*3^(1/2)*(-f2+f1);
C(7,45)=-1/6*3^(1/2)*(-f2+f1);


C(8,8)=1/3*(f1);
C(8,9)=1/3*(-1)*(-f1);
C(8,10)=-1/3*3^(1/2)*i*(-1)*(-f1);
C(8,12)=-1/3*3^(1/2)*i*(-1)*(-f1);
C(8,13)=-1/6*6^(1/2)*i*(-1)*(-f1);
C(8,15)=-1/6*6^(1/2)*i*(-1)*(-f1);
C(8,19)=1/18*18^(1/2)*(-1)*(f1);
C(8,20)=1/18*18^(1/2)*(-1)*(-2*f1);
C(8,21)=1/18*18^(1/2)*(-1)*(-f1);
C(8,22)=-1/6*6^(1/2)*i*(-1)*(-f2);
C(8,23)=-1/6*6^(1/2)*i*(-1)*(-f2-f1);
C(8,24)=1/6*6^(1/2)*i*(-1)*(-f2+f1);
C(8,25)=-1/6*3^(1/2)*i*(-1)*(-f2);
C(8,26)=-1/6*3^(1/2)*i*(-1)*(-f2-f1);
C(8,27)=1/6*3^(1/2)*i*(-1)*(-f2+f1);
C(8,31)=-1/6*3^(1/2)*i*(-1)*(-f2);
C(8,32)=-1/6*3^(1/2)*i*(-1)*(-f2+f1);
C(8,33)=1/6*3^(1/2)*i*(-1)*(-f2-f1);
C(8,34)=1/6*6^(1/2)*(-1)*(-f2);
C(8,35)=1/6*6^(1/2)*(-1)*(-f2-f1);
C(8,36)=1/6*6^(1/2)*(-1)*(-f2-f1);
C(8,37)=1/6*3^(1/2)*(-1)*(-f2);
C(8,38)=1/6*3^(1/2)*(-1)*(-f2-f1);
C(8,39)=-1/6*3^(1/2)*(-1)*(-f2-f1);
C(8,43)=1/6*3^(1/2)*(-1)*(-f2);
C(8,44)=1/6*3^(1/2)*(-1)*(-f2+f1);
C(8,45)=-1/6*3^(1/2)*(-1)*(-f2+f1);


C(9,9)=1/3*(f1);
C(9,10)=-1/3*3^(1/2)*i*(-1)*(-f1);
C(9,12)=-1/3*3^(1/2)*i*(-1)*(-f1);
C(9,13)=-1/6*6^(1/2)*i*(-1)*(-f1);
C(9,15)=-1/6*6^(1/2)*i*(-1)*(-f1);
C(9,19)=1/18*18^(1/2)*(-1)*(f1);
C(9,20)=1/18*18^(1/2)*(-1)*(-2*f1);
C(9,21)=1/18*18^(1/2)*(-1)*(-f1);
C(9,22)=-1/6*6^(1/2)*i*(-1)*(-f2);
C(9,23)=-1/6*6^(1/2)*i*(-1)*(-f2-f1);
C(9,24)=1/6*6^(1/2)*i*(-1)*(-f2+f1);
C(9,25)=-1/6*3^(1/2)*i*(-1)*(-f2);
C(9,26)=-1/6*3^(1/2)*i*(-1)*(-f2-f1);
C(9,27)=1/6*3^(1/2)*i*(-1)*(-f2+f1);
C(9,31)=-1/6*3^(1/2)*i*(-1)*(-f2);
C(9,32)=-1/6*3^(1/2)*i*(-1)*(-f2+f1);
C(9,33)=1/6*3^(1/2)*i*(-1)*(-f2-f1);
C(9,34)=1/6*6^(1/2)*(-1)*(-f2);
C(9,35)=1/6*6^(1/2)*(-1)*(-f2-f1);
C(9,36)=1/6*6^(1/2)*(-1)*(-f2-f1);
C(9,37)=1/6*3^(1/2)*(-1)*(-f2);
C(9,38)=1/6*3^(1/2)*(-1)*(-f2-f1);
C(9,39)=-1/6*3^(1/2)*(-1)*(-f2-f1);
C(9,43)=1/6*3^(1/2)*(-1)*(-f2);
C(9,44)=1/6*3^(1/2)*(-1)*(-f2+f1);
C(9,45)=-1/6*3^(1/2)*(-1)*(-f2+f1);

C(10,10)=(-f1);
C(10,12)=(-f1);
C(10,13)=-1/2*2^(1/2)*(f1);
C(10,15)=-1/2*2^(1/2)*(f1);
C(10,19)=-1/6*6^(1/2)*(i)*(-f1);
C(10,20)=-1/6*6^(1/2)*(i)*(2*f1);
C(10,21)=-1/6*6^(1/2)*(i)*(f1);
C(10,22)=-1/2*2^(1/2)*(f2);
C(10,23)=-1/2*2^(1/2)*(f2+f1);
C(10,24)=1/2*2^(1/2)*(f2-f1);
C(10,25)=-1/2*(f2);
C(10,26)=-1/2*(f2+f1);
C(10,27)=1/2*(f2-f1);
C(10,31)=-1/2*(f2);
C(10,32)=-1/2*(f2-f1);
C(10,33)=1/2*(f2+f1);
C(10,34)=-1/2*2^(1/2)*(i)*(f2);
C(10,35)=-1/2*2^(1/2)*(i)*(f2+f1);
C(10,36)=1/2*2^(1/2)*(i)*(f2+f1);
C(10,37)=-1/2*(i)*(f2);
C(10,38)=-1/2*(i)*(f2+f1);
C(10,39)=1/2*(i)*(f2+f1);
C(10,43)=-1/2*(i)*(f2);
C(10,44)=-1/2*(i)*(f2+-f1);
C(10,45)=1/2*(i)*(f2-f1);


C(12,12)=(-f1);
C(12,13)=-1/2*2^(1/2)*(f1);
C(12,15)=-1/2*2^(1/2)*(f1);
C(12,19)=-1/6*6^(1/2)*(i)*(-f1);
C(12,20)=-1/6*6^(1/2)*(i)*(2*f1);
C(12,21)=-1/6*6^(1/2)*(i)*(f1);
C(12,22)=-1/2*2^(1/2)*(f2);
C(12,23)=-1/2*2^(1/2)*(f2+f1);
C(12,24)=1/2*2^(1/2)*(f2-f1);
C(12,25)=-1/2*(f2);
C(12,26)=-1/2*(f2+f1);
C(12,27)=1/2*(f2-f1);
C(12,31)=-1/2*(f2);
C(12,32)=-1/2*(f2-f1);
C(12,33)=1/2*(f2+f1);
C(12,34)=-1/2*2^(1/2)*(i)*(f2);
C(12,35)=-1/2*2^(1/2)*(i)*(f2+f1);
C(12,36)=1/2*2^(1/2)*(i)*(f2+f1);
C(12,37)=-1/2*(i)*(f2);
C(12,38)=-1/2*(i)*(f2+f1);
C(12,39)=1/2*(i)*(f2+f1);
C(12,43)=-1/2*(i)*(f2);
C(12,44)=-1/2*(i)*(f2+-f1);
C(12,45)=1/2*(i)*(f2-f1);

C(13,13)=-1/2*(f1);
C(13,15)=-1/2*2^(1/2)*1/2*2^(1/2)*(f1);
C(13,19)=-1/6*6^(1/2)*(i)*1/2*2^(1/2)*(-f1);
C(13,20)=-1/6*6^(1/2)*(i)*1/2*2^(1/2)*(2*f1);
C(13,21)=-1/6*6^(1/2)*(i)*1/2*2^(1/2)*(f1);
C(13,22)=-1/2*2^(1/2)*1/2*2^(1/2)*(f2);
C(13,23)=-1/2*2^(1/2)*1/2*2^(1/2)*(f2+f1);
C(13,24)=1/2*2^(1/2)*1/2*2^(1/2)*(f2-f1);
C(13,25)=-1/2*1/2*2^(1/2)*(f2);
C(13,26)=-1/2*1/2*2^(1/2)*(f2+f1);
C(13,27)=1/2*1/2*2^(1/2)*(f2-f1);
C(13,31)=-1/2*1/2*2^(1/2)*(f2);
C(13,32)=-1/2*1/2*2^(1/2)*(f2-f1);
C(13,33)=1/2*1/2*2^(1/2)*(f2+f1);
C(13,34)=-1/2*2^(1/2)*(i)*1/2*2^(1/2)*(f2);
C(13,35)=-1/2*2^(1/2)*(i)*1/2*2^(1/2)*(f2+f1);
C(13,36)=1/2*2^(1/2)*(i)*1/2*2^(1/2)*(f2+f1);
C(13,37)=-1/2*(i)*1/2*2^(1/2)*(f2);
C(13,38)=-1/2*(i)*1/2*2^(1/2)*(f2+f1);
C(13,39)=1/2*(i)*1/2*2^(1/2)*(f2+f1);
C(13,43)=-1/2*(i)*1/2*2^(1/2)*(f2);
C(13,44)=-1/2*(i)*1/2*2^(1/2)*(f2+-f1);
C(13,45)=1/2*(i)*1/2*2^(1/2)*(f2-f1);

C(15,15)=-1/2*(f1);
C(15,19)=-1/6*6^(1/2)*(i)*1/2*2^(1/2)*(-f1);
C(15,20)=-1/6*6^(1/2)*(i)*1/2*2^(1/2)*(2*f1);
C(15,21)=-1/6*6^(1/2)*(i)*1/2*2^(1/2)*(f1);
C(15,22)=-1/2*2^(1/2)*1/2*2^(1/2)*(f2);
C(15,23)=-1/2*2^(1/2)*1/2*2^(1/2)*(f2+f1);
C(15,24)=1/2*2^(1/2)*1/2*2^(1/2)*(f2-f1);
C(15,25)=-1/2*1/2*2^(1/2)*(f2);
C(15,26)=-1/2*1/2*2^(1/2)*(f2+f1);
C(15,27)=1/2*1/2*2^(1/2)*(f2-f1);
C(15,31)=-1/2*1/2*2^(1/2)*(f2);
C(15,32)=-1/2*1/2*2^(1/2)*(f2-f1);
C(15,33)=1/2*1/2*2^(1/2)*(f2+f1);
C(15,34)=-1/2*2^(1/2)*(i)*1/2*2^(1/2)*(f2);
C(15,35)=-1/2*2^(1/2)*(i)*1/2*2^(1/2)*(f2+f1);
C(15,36)=1/2*2^(1/2)*(i)*1/2*2^(1/2)*(f2+f1);
C(15,37)=-1/2*(i)*1/2*2^(1/2)*(f2);
C(15,38)=-1/2*(i)*1/2*2^(1/2)*(f2+f1);
C(15,39)=1/2*(i)*1/2*2^(1/2)*(f2+f1);
C(15,43)=-1/2*(i)*1/2*2^(1/2)*(f2);
C(15,44)=-1/2*(i)*1/2*2^(1/2)*(f2+-f1);
C(15,45)=1/2*(i)*1/2*2^(1/2)*(f2-f1);


C(19,19)=1/6*(f1);
C(19,20)=1/6*(-2*f1);
C(19,21)=1/6*(-f1);
C(19,22)=-1/12*12^(1/2)*(i)*(-f2);
C(19,23)=-1/12*12^(1/2)*(i)*(-f2-f1);
C(19,24)=1/12*12^(1/2)*(i)*(-f2+f1);
C(19,25)=-1/12*6^(1/2)*(i)*(-f2);
C(19,26)=-1/12*6^(1/2)*(i)*(-f2-f1);
C(19,27)=1/12*6^(1/2)*(i)*(-f2+f1);
C(19,31)=-1/12*6^(1/2)*(i)*(-f2);
C(19,32)=-1/12*6^(1/2)*(i)*(-f2+f1);
C(19,33)=1/12*6^(1/2)*(i)*(-f2-f1);
C(19,34)=1/12*12^(1/2)*(-f2);
C(19,35)=1/12*12^(1/2)*(-f2-f1);
C(19,36)=-1/12*12^(1/2)*(-f2-f1);
C(19,37)=1/12*6^(1/2)*(-f2);
C(19,38)=1/12*6^(1/2)*(-f2-f1);
C(19,39)=-1/12*6^(1/2)*(-f2-f1);
C(19,43)=1/12*6^(1/2)*(-f2);
C(19,44)=1/12*6^(1/2)*(-f2+f1);
C(19,45)=-1/12*6^(1/2)*(-f2+f1);


C(20,20)=1/6*(2*f1);
C(20,21)=1/6*(-2)*(-f1);
C(20,22)=-1/12*12^(1/2)*(i)*(-2)*(-f2);
C(20,23)=-1/12*12^(1/2)*(i)*(-2)*(-f2-f1);
C(20,24)=1/12*12^(1/2)*(i)*(-2)*(-f2+f1);
C(20,25)=-1/12*6^(1/2)*(i)*(-2)*(-f2);
C(20,26)=-1/12*6^(1/2)*(i)*(-2)*(-f2-f1);
C(20,27)=1/12*6^(1/2)*(i)*(-2)*(-f2+f1);
C(20,31)=-1/12*6^(1/2)*(i)*(-2)*(-f2);
C(20,32)=-1/12*6^(1/2)*(i)*(-2)*(-f2+f1);
C(20,33)=1/12*6^(1/2)*(i)*(-2)*(-f2-f1);
C(20,34)=1/12*12^(1/2)*(-2)*(-f2);
C(20,35)=1/12*12^(1/2)*(-2)*(-f2-f1);
C(20,36)=-1/12*12^(1/2)*(-2)*(-f2-f1);
C(20,37)=1/12*6^(1/2)*(-2)*(-f2);
C(20,38)=1/12*6^(1/2)*(-2)*(-f2-f1);
C(20,39)=-1/12*6^(1/2)*(-2)*(-f2-f1);
C(20,43)=1/12*6^(1/2)*(-2)*(-f2);
C(20,44)=1/12*6^(1/2)*(-2)*(-f2+f1);
C(20,45)=-1/12*6^(1/2)*(-2)*(-f2+f1);

C(21,21)=1/6*(f1);
C(21,22)=-1/12*12^(1/2)*(i)*(-1)*(-f2);
C(21,23)=-1/12*12^(1/2)*(i)*(-1)*(-f2-f1);
C(21,24)=1/12*12^(1/2)*(i)*(-1)*(-f2+f1);
C(21,25)=-1/12*6^(1/2)*(i)*(-1)*(-f2);
C(21,26)=-1/12*6^(1/2)*(i)*(-1)*(-f2-f1);
C(21,27)=1/12*6^(1/2)*(i)*(-1)*(-f2+f1);
C(21,31)=-1/12*6^(1/2)*(i)*(-1)*(-f2);
C(21,32)=-1/12*6^(1/2)*(i)*(-1)*(-f2+f1);
C(21,33)=1/12*6^(1/2)*(i)*(-1)*(-f2-f1);
C(21,34)=1/12*12^(1/2)*(-1)*(-f2);
C(21,35)=1/12*12^(1/2)*(-1)*(-f2-f1);
C(21,36)=-1/12*12^(1/2)*(-1)*(-f2-f1);
C(21,37)=1/12*6^(1/2)*(-1)*(-f2);
C(21,38)=1/12*6^(1/2)*(-1)*(-f2-f1);
C(21,39)=-1/12*6^(1/2)*(-1)*(-f2-f1);
C(21,43)=1/12*6^(1/2)*(-1)*(-f2);
C(21,44)=1/12*6^(1/2)*(-1)*(-f2+f1);
C(21,45)=-1/12*6^(1/2)*(-1)*(-f2+f1);

C(22,23)=-1/2*(f2);
C(22,24)=-1/2*(f2);
C(22,26)=-1/2*4^(1/2)*(f2);
C(22,27)=-1/2*4^(1/2)*(f2);
C(22,32)=-1/2*4^(1/2)*(-f2);
C(22,33)=1/2*4^(1/2)*(f2);
C(22,35)=-1/2*i*(f2);
C(22,36)=1/2*i*(f2);
C(22,38)=-1/4*2^(1/2)*i*(f2);
C(22,39)=1/4*2^(1/2)*i*(f2);
C(22,44)=-1/4*2^(1/2)*i*(-f2);
C(22,45)=-1/4*2^(1/2)*i*(f2);

C(23,23)=-1/2*(2*f2+f1);
C(23,24)=-1/2*(f1);
C(23,25)=-1/4*2^(1/2)*(f2);
C(23,26)=-1/4*2^(1/2)*(f1+2*f2);
C(23,27)=1/4*2^(1/2)*(-f1);
C(23,31)=-1/4*2^(1/2)*(f2);
C(23,32)=-1/4*2^(1/2)*(-f1);
C(23,33)=1/4*2^(1/2)*(f1+2*f2);
C(23,34)=-1/2*i*(f2);
C(23,35)=-1/2*i*(f1+2*f2);
C(23,36)=1/2*i*(f1+2*f2);
C(23,37)=-1/4*2^(1/2)*i*(f2);
C(23,38)=-1/4*2^(1/2)*i*(f1+2*f2);
C(23,39)=1/4*2^(1/2)*i*(f1+2*f2);
C(23,43)=-1/4*2^(1/2)*i*(f2);
C(23,44)=-1/4*2^(1/2)*i*(-f1);
C(23,45)=1/4*2^(1/2)*i*(-f1);

C(24,25)=1/4*2^(1/2)*(-f2);
C(24,26)=1/4*2^(1/2)*(-f1);
C(24,27)=-1/4*2^(1/2)*(f1-2*f2);
C(24,31)=1/4*2^(1/2)*(-f2);
C(24,32)=1/4*2^(1/2)*(f1-2*f2);
C(24,33)=-1/4*2^(1/2)*(-f1);
C(24,34)=1/2*i*(-f2);
C(24,35)=1/2*i*(-f1);
C(24,36)=-1/2*i*(-f1);
C(24,37)=1/4*2^(1/2)*i*(-f2);
C(24,38)=1/4*2^(1/2)*i*(-f1);
C(24,39)=-1/4*2^(1/2)*i*(-f1);
C(24,43)=1/4*2^(1/2)*i*(-f2);
C(24,44)=1/4*2^(1/2)*i*(f1-2*f2);
C(24,45)=-1/4*2^(1/2)*i*(-f1-2*f2);

C(25,26)=-1/2*4^(1/2)*1/2*2^(1/2)*(f2);
C(25,27)=-1/2*4^(1/2)*1/2*2^(1/2)*(f2);
C(25,32)=-1/2*4^(1/2)*1/2*2^(1/2)*(-f2);
C(25,33)=1/2*4^(1/2)*1/2*2^(1/2)*(f2);
C(25,35)=-1/2*i*1/2*2^(1/2)*(f2);
C(25,36)=1/2*i*1/2*2^(1/2)*(f2);
C(25,38)=-1/4*2^(1/2)*i*1/2*2^(1/2)*(f2);
C(25,39)=1/4*2^(1/2)*i*1/2*2^(1/2)*(f2);
C(25,44)=-1/4*2^(1/2)*i*1/2*2^(1/2)*(-f2);
C(25,45)=-1/4*2^(1/2)*i*1/2*2^(1/2)*(f2);


C(26,26)=-1/4*2^(1/2)*1/2*2^(1/2)*(f1+2*f2);
C(26,27)=1/4*2^(1/2)*1/2*2^(1/2)*(-f1);
C(26,31)=-1/4*2^(1/2)*1/2*2^(1/2)*(f2);
C(26,32)=-1/4*2^(1/2)*1/2*2^(1/2)*(-f1);
C(26,33)=1/4*2^(1/2)*1/2*2^(1/2)*(f1+2*f2);
C(26,34)=-1/2*i*1/2*2^(1/2)*(f2);
C(26,35)=-1/2*i*1/2*2^(1/2)*(f1+2*f2);
C(26,36)=1/2*i*1/2*2^(1/2)*(f1+2*f2);
C(26,37)=-1/4*2^(1/2)*i*1/2*2^(1/2)*(f2);
C(26,38)=-1/4*2^(1/2)*i*1/2*2^(1/2)*(f1+2*f2);
C(26,39)=1/4*2^(1/2)*i*1/2*2^(1/2)*(f1+2*f2);
C(26,43)=-1/4*2^(1/2)*i*1/2*2^(1/2)*(f2);
C(26,44)=-1/4*2^(1/2)*i*1/2*2^(1/2)*(-f1);
C(26,45)=1/4*2^(1/2)*i*1/2*2^(1/2)*(-f1);


C(27,27)=-1/4*(f1-2*f2);
C(27,31)=1/4*2^(1/2)*1/2*2^(1/2)*(-f2);
C(27,32)=1/4*2^(1/2)*1/2*2^(1/2)*(f1-2*f2);
C(27,33)=-1/4*2^(1/2)*1/2*2^(1/2)*(-f1);
C(27,34)=1/2*i*1/2*2^(1/2)*(-f2);
C(27,35)=1/2*i*1/2*2^(1/2)*(-f1);
C(27,36)=-1/2*i*1/2*2^(1/2)*(-f1);
C(27,37)=1/4*2^(1/2)*i*1/2*2^(1/2)*(-f2);
C(27,38)=1/4*2^(1/2)*i*1/2*2^(1/2)*(-f1);
C(27,39)=-1/4*2^(1/2)*i*1/2*2^(1/2)*(-f1);
C(27,43)=1/4*2^(1/2)*i*1/2*2^(1/2)*(-f2);
C(27,44)=1/4*2^(1/2)*i*1/2*2^(1/2)*(f1-2*f2);
C(27,45)=-1/4*2^(1/2)*i*1/2*2^(1/2)*(-f1-2*f2);

C(31,32)=-1/2*4^(1/2)*1/2*2^(1/2)*(-f2);
C(31,33)=1/2*4^(1/2)*1/2*2^(1/2)*(f2);
C(31,35)=-1/2*i*1/2*2^(1/2)*(f2);
C(31,36)=1/2*i*(1/2*2^(1/2)*f2);
C(31,38)=-1/4*2^(1/2)*i*(f2);
C(31,39)=1/4*2^(1/2)*i*1/2*2^(1/2)*(f2);
C(31,44)=-1/4*2^(1/2)*i*1/2*2^(1/2)*(-f2);
C(31,45)=-1/4*2^(1/2)*i*1/2*2^(1/2)*(f2);

C(32,32)=-1/4*(f1-2*f2);
C(32,33)=-1/4*2^(1/2)*1/2*2^(1/2)*(-1)*(-f1);
C(32,34)=1/2*i*1/2*2^(1/2)*(-1)*(-f2);
C(32,35)=1/2*i*1/2*2^(1/2)*(-1)*(-f1);
C(32,36)=-1/2*i*1/2*2^(1/2)*(-1)*(-f1);
C(32,37)=1/4*2^(1/2)*i*1/2*2^(1/2)*(-1)*(-f2);
C(32,38)=1/4*2^(1/2)*i*1/2*2^(1/2)*(-1)*(-f1);
C(32,39)=-1/4*2^(1/2)*i*1/2*2^(1/2)*(-1)*(-f1);
C(32,43)=1/4*2^(1/2)*i*1/2*2^(1/2)*(-1)*(-f2);
C(32,44)=1/4*2^(1/2)*i*1/2*2^(1/2)*(-1)*(f1-2*f2);
C(32,45)=-1/4*2^(1/2)*i*1/2*2^(1/2)*(-1)*(-f1-2*f2);

C(33,33)=-1/4*(f1+2*f2);
C(33,34)=-1/2*i*1/2*2^(1/2)*(-1)*(f2);
C(33,35)=-1/2*i*1/2*2^(1/2)*(-1)*(f1+2*f2);
C(33,36)=1/2*i*1/2*2^(1/2)*(-1)*(f1+2*f2);
C(33,37)=-1/4*2^(1/2)*i*1/2*2^(1/2)*(-1)*(f2);
C(33,38)=-1/4*2^(1/2)*i*(1/2*2^(1/2)*(-1)*f1+2*f2);
C(33,39)=1/4*2^(1/2)*i*(1/2*2^(1/2)*(-1)*f1+2*f2);
C(33,43)=-1/4*2^(1/2)*i*1/2*2^(1/2)*(-1)*(f2);
C(33,44)=-1/4*2^(1/2)*i*1/2*2^(1/2)*(-1)*(-f1);
C(33,45)=1/4*2^(1/2)*i*(1/2*2^(1/2)*(-1)*-f1);


C(34,35)=1/2*i*(-i)*(f2);
C(34,36)=-1/2*i*(-i)*(f2);
C(34,38)=1/4*2^(1/2)*i*(-i)*(f2);
C(34,39)=-1/4*2^(1/2)*i*(-i)*(f2);
C(34,44)=1/4*2^(1/2)*i*(-i)*(-f2);
C(34,45)=1/4*2^(1/2)*i*(-i)*(f2);

C(35,35)=1/2*(f1+2*f2);
C(35,36)=1/2*i*(i)*(f1+2*f2);
C(35,37)=-1/4*2^(1/2)*i*(i)*(f2);
C(35,38)=-1/4*2^(1/2)*i*(i)*(f1+2*f2);
C(35,39)=1/4*2^(1/2)*i*(i)*(f1+2*f2);
C(35,43)=-1/4*2^(1/2)*i*(i)*(f2);
C(35,44)=-1/4*2^(1/2)*i*(i)*(-f1);
C(35,45)=1/4*2^(1/2)*i*(i)*(-f1);

C(36,36)=1/2*(f1+2*f2);
C(36,37)=-1/4*2^(1/2)*i*(i)*(-1)*(f2);
C(36,38)=-1/4*2^(1/2)*i*(i)*((-1)*f1+2*f2);
C(36,39)=1/4*2^(1/2)*i*(i)*(-1)*(f1+2*f2);
C(36,43)=-1/4*2^(1/2)*i*(i)*(-1)*(f2);
C(36,44)=-1/4*2^(1/2)*i*(i)*(-1)*(-f1);
C(36,45)=1/4*2^(1/2)*i*(i)*(-1)*(-f1);

C(37,38)=-1/4*2^(1/2)*i*1/2*2^(1/2)*i*(f2);
C(37,39)=1/4*2^(1/2)*i*1/2*2^(1/2)*i*(f2);
C(37,44)=-1/4*2^(1/2)*i*1/2*2^(1/2)*i*(-f2);
C(37,45)=-1/4*2^(1/2)*i*1/2*2^(1/2)*i*(f2);


C(38,38)=1/4*(f1+2*f2);
C(38,39)=1/4*2^(1/2)*i*1/2*2^(1/2)*i*(f1+2*f2);
C(38,43)=-1/4*2^(1/2)*i*1/2*2^(1/2)*i*(f2);
C(38,44)=-1/4*2^(1/2)*i*1/2*2^(1/2)*i*(-f1);
C(38,45)=1/4*2^(1/2)*i*1/2*2^(1/2)*i*(-f1);

C(39,39)=1/4*(f1+2*f2);
C(39,43)=-1/4*2^(1/2)*i*1/2*2^(1/2)*i*(-1)*(f2);
C(39,44)=-1/4*2^(1/2)*i*1/2*2^(1/2)*i*(-1)*(-f1);
C(39,45)=1/4*2^(1/2)*i*1/2*2^(1/2)*i*(-1)*(-f1);

C(43,44)=1/4*(-f2);
C(43,45)=1/4*(f2);

C(44,44)=1/4*(f1-2*f2);
C(44,45)=-1/4*(f1-2*f2);

C(45,45)=1/4*(f1-2*f2);


D=A+B+C;

[V1,E1]=eig(D);
E1