Electricity, Complex Systems, Numerical Analysis, Algorithms, Network.
J=Jacobian(Y_eq[1:(nrow(Y_eq)-1),1:(ncol(Y_eq)-1)],U_old[1:(nrow(Y_eq)-1)],phi_old[1:(nrow(Y_eq)-1)]) delta_var=GS(-J,c(Re(delta_S)[1:(nrow(Y_eq)-1)],Im(delta_S)[1:(nrow(Y_eq)-1)])) phi_new=phi_old+c(delta_var[1:(nrow(Y_eq)-1)],0) U_new=U_old+c(delta_var[nrow(Y_eq):length(delta_var)],0)