p1 = [0,0];
p2 = [-1,-2];
p3 = [3,3];
p4 = [5,4];
y1=0;
y2=0;
y3=1;
y4=1;
function s = sigma(x)
s = (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(-5);
z21 = (w*h1)';
disp(z21);
h21 = [sigma(z21(1)), sigma(z21(2)), sigma(z21(3)), sigma(z21(4))]';
disp(h21);
function g = gradient(h1, y)
g =