Electricity, Complex Systems, Numerical Analysis, Algorithms, Network.
var=c(U_new,phi_new) } return(var) } The iterative function for the second method is: #Jacobian, condition 2 Jacobian<-function(Y_eq,U_old,phi_old){ J=matrix(nrow=2*length(U_old),ncol=2*length(U_old)) for(i in 1:length(U_old)){ #partial derivatives with respect to phi #P