該程序?qū)崿F(xiàn)簡(jiǎn)單的牛拉法極坐標(biāo)下的潮流計(jì)算。首先將支路導(dǎo)納矩陣轉(zhuǎn)化為節(jié)點(diǎn)導(dǎo)納矩陣，初始化功率初值，電壓初值，接著形成Jacobian矩陣，對(duì)矩陣進(jìn)行LDU分解，前代回代求解修正量，最后得出各個(gè)節(jié)點(diǎn)的電壓幅值與角度。