function R = R_matrix_IM(alpha)

R(1,51) = 1; R(1,52) = 1;
R(2,53) = 1; R(2,58) = 1;
R(3,53) = 0; R(3,58) = 0;
R(4,54) = 1; R(4,60) = 1;
R(5,79) = 1; R(5,80) = 1;
R(6,55) = 1; R(6,56) = 1;
R(7,57) = 1; R(7,76) = 1;
R(8,57) = 0; R(8,76) = 0;
R(9,77) = 1; R(9,59) = 1; 
R(10,77) = 0; R(10,59) = 0; 
R(11,61) = 1; R(11,81) = 1;
R(12,65) = 1; R(12,82) = 1;
R(13,78) = 1; R(13,83) = 1;
R(14,78) = 0; R(14,62) = 0;
R(15,63) = 0; R(15,67) = 0;
R(16,63) = 1; R(16,84) = 1;
R(17,64) = 1; R(17,69) = 1;
R(18,66) = 0; R(18,73) = 0;
R(19,66) = 1; R(19,73) = 1;
R(20,68) = 0; R(20,74) = 0;
R(21,68) = 1; R(21,86) = 1;
R(22,70) = 1; R(22,72) = 1;
R(23,71) = 0; R(23,75) = 0;
R(24,85) = 1; R(24,75) = 1;
R(25,1) = 1; R(25,2) = 1;
R(26,22) = 1; R(26,24) = 1; 
R(27,23) = 1; R(27,21) = 1; R(27,12) = 1;
R(28,3) = 1; R(28,4) = 1;  R(28,7) = 1; 
R(29,6) = 1; R(29,8) = 1;  R(29,9) = 1; 
R(30,10) = 1; R(30,18) = 1;  
R(31,14) = 1;   
R(32,15) = 1;   
R(33,17) = 1; R(33,19) = 1;  R(33,20) = 1; 
R(34,5) = 1; R(34,11) = 1;  R(34,13) = 1;  R(34,16) = 1;
R(35,25) = 1; 
R(36,26) = 1; 
R(37,27) = 1/2;
R(38,27) = 1/2;
R(39,28) = 1/2;
R(40,28) = 1/2;
R(41,34) = 1/2;
R(42,34) = 1/2;
R(43,29) = 1/2;
R(44,29) = 1/2;
R(45,30) = 1/2;
R(46,30) = 1/2;
R(47,31) = 1;
R(48,32) = 1;
R(49,33) = 1/2;
R(50,33) = 1/2;
R(51,35) = 1-alpha;
R(52,39) = 1-alpha;
R(53,40) = 1-alpha;
R(54,42) = 1-alpha;
R(55,41) = 1-alpha;
R(56,43) = 1-alpha;
R(57,44) = 1-alpha;
R(58,48) = 1-alpha;
R(59,50) = 1-alpha;
R(60,49) = 1-alpha; 
R(61,38) = 1-alpha; 
R(62,37) = 0;
R(63,40) = alpha;
R(64,42) = alpha;
R(65,41) = alpha; 
R(66,44) = alpha;
R(67,48) = 0;
R(68,50) = alpha;
R(69,49) = alpha;
R(70,38) = alpha;
R(71,37) = 0;
R(72,36) = alpha;
R(73,45) = alpha;
R(74,46) = 0;
R(75,47) = alpha;
R(76,45) = 1-alpha;
R(77,46) = 1-alpha;
R(78,47) = 1-alpha;
R(79,35) = alpha;
R(80,39) = alpha;
R(81,36) = 1-alpha;
R(82,43) = alpha;
R(83,37) = 1-alpha;
R(84,48) = alpha;
R(85,37) = alpha;
R(86,46) = alpha;






























