function y=fc(x)
global inc
%V1= 1.01316;j1=0.18989;p1=3.38;q1=1.34045;V2=1.01299;j2=0.19091;p2=3.4;q2=1.3501;
V1= inc(1);j1=inc(2);p1=inc(3);q1=inc(4);V2=inc(5);j2=inc(6);p2=inc(7);q2=inc(8);
y(1)=(x(1)*V1*cos(x(3)-x(4)-(j2-j1)/6))/x(2)-V1^2*cos(x(3))/x(2)-p1;
y(2)=(x(1)*V1*sin(x(3)-x(4)-(j2-j1)/6))/x(2)-V1^2*sin(x(3))/x(2)-q1;
%y(1)=(x(1)*V1*cos(x(3)-x(4)))/x(2)-V1^2*cos(x(3))/x(2)-p1;
%y(2)=(x(1)*V1*sin(x(3)-x(4)))/x(2)-V1^2*sin(x(3))/x(2)-q1;
y(3)=(x(1)*V2*cos(x(3)-x(4)))/x(2)-V2^2*cos(x(3))/x(2)-p2;
y(4)=(x(1)*V2*sin(x(3)-x(4)))/x(2)-V2^2*sin(x(3))/x(2)-q2;
y=[y(1) y(2) y(3) y(4)];