Need opinion on why a simple Newton-Raphson algorithm is not converging correctly

I'm trying to solve a simple 4 bus voltage equation, where the powers for the buses and admittances between the buses are known. I think this code is close to solving, but something just isn't quite right as it gives the wrong solution. The values for the voltage I'm expecting are % What the voltages should be: % The voltage for bus 1 is 1.0000 V real and 0.0000 V reactive, with angle 0.00 degrees % The voltage for bus 2 is 0.9085 V real and -0.0476 V reactive, with angle -3.00 degrees % The voltage for bus 3 is 0.8925 V real and -0.0369 V reactive, with angle -2.37 degrees % The voltage for bus 4 is 0.9208 V real and -0.0377 V reactive, with angle -2.34 degrees % % The powers should be: % Power of bus 1 is -0.4363 W and -0.5159 VAR. % Power of bus 2 is 0.2000 W and 0.1500 VAR. % Power of bus 3 is 0.0200 W and 0.1500 VAR. % Power of bus 4 is 0.2000 W and 0.1500 VAR.

Please see the attached file.

Thank you for any input.

[Walter Roberson]

No files attached