NR_power-Flow
代码说明:
x0 = ones(2,2) Make a starting guess at the solution options = optimset( Display , off ) Turn off Display [x,Fval,exitflag] = fsolve(@myfun,x0,options) The solution is x = -0.1291 0.8602 1.2903 1.1612 Fval = 1.0e-009 * -0.1619 0.0776 0.1161 -0.0469 exitflag = 1 and the residual is close to zero. sum(sum(Fval.*Fval)) ans = 4.7915e-020
下载说明:请别用迅雷下载,失败请重下,重下不扣分!