p1 = [0,0];
p2 = [1,0];
p3 = [5,5];
p4 = [10,0];
y1=0;
y2=0;
y3=1;
y4=1;
function sigma(x)
return (1/(1+exp(-x)));
endfunction
y = [y1,y2,y3,y4];
disp(y')
w11 = 1;
w12 = 1;
w13 = 1;
w = [w11, w12, w13];
h11 = [p1(1),p1(2),1];
h12 = [p2(1),p2(2),1];
h13 = [p3(1),p3(2),1];
h14 = [p4(1),p4(2),1];
h1 = [h11',h12',h13',h14']
s = sigma(0)