function y = f1(x)
L1 = 3;
L2 = 3*sqrt(2);
L3 = 3;
gamma = pi/4;
p1 = 5;
p2 = 5;
p3 = 3;
x1 = 5;
x2 = 0;
y2 = 6;
A2 = (L3*cos(x))-x1;
A3 = L2*((cos(x)*cos(gamma))-(sin(x)*sin(gamma)))-x2;
B2 = L3*sin(x);
B3 = L2*((cos(x)*sin(gamma))+(sin(x)*cos(gamma)))-y2;
N1 = B3.*(p2^2-p1^2-A2.^2-B2.^2)-B2.*(p3^2-p1^2-A3.^2-B3.^2);
N2 = -A3.*(p2^2-p1^2-A2.^2-B2.^2)+A2.*(p3^2-p1^2-A3.^2-B3.^2);
D  = 2*(A2.*B3-B2.*A3);
y = N1.^2+N2.^2-(p1^2*D.^2);