function [r,J] = loadFlowOne(x,V0,Z,S)
%S = 
r = zeros(2,1);
J = zeros(2);

V1 = x(1) + 1i*x(2);
I = (V0 - V1)/Z;
S1 = conj(I)*V1;

r(1) = real(S1 - S);
r(2) = imag(S1 - S);

dS1dV1 = conj(I)*[1 1i] - V1/conj(Z)*[1 -1i];
J(1,:) = real(dS1dV1);
J(2,:) = imag(dS1dV1);





