# p1 = (1, 1)
# p2 = (3, 3)
# p3 = (5, 1)
# p4 = (2, 8)
# p5 = (-5, 1)
# p6 = (-7, 3)
# p7 = (-1, 1)
# p8 = (-17, 1)
# p9 = (1, -4)
# p10 = (3, -12)
# p11 = (12, -1)
# p12 = (1, -5)
# p13 = (-4, -1)
# p14 = (-5, -5)
# p15 = (-12, -3)
# p16 = (-2, -7)
H1 = [1,1,1;
3,3,1;
5,1,1;
2,8,1;
-5,1,1;
-7,3,1;
-1,1,1;
-17,1,1;
1,-4,1;
3,-12,1;
12,-1,1;
1,-5,1;
-4,-1,1;
-5,-5,1;
-12,-3,1;
-2,-7,1;
];
Y = [0,0;
0,0;
0,0;
0,0;
0,1;
0,1;
0,1;
0,1;
1,0;
1,0;
1,0;
1,0;
1,1;
1,1;
1,1;
1,1;
];
H1=H1';
#disp(h1);
W = [1,1;
1,1;
1,1];
W=W';
Z2=W*H1;
str1=Z2(1,:);
str2=Z2(2,:);
#disp(str2);
str1_sum=0;
str2_sum=0;
for i=1:16;
str1_sum+=exp(str1(i));
endfor;
for i=1:16;
str2_sum+=exp(str2(i));
endfor;
for j=1:16;
str1_h(j)=exp(str1(j))/str1_sum;
endfor;
for j=1:16;
str2_h(j)=exp(str2(j))/str2_sum;
endfor;
H2 = [str1_h; str2_h];
Y=Y';
disp(Y);