該程序實現簡單的牛拉法極坐標下的潮流計算。首先將支路導納矩陣轉化為節點導納矩陣，初始化功率初值，電壓初值，接著形成Jacobian矩陣，對矩陣進行LDU分解，前代回代求解修正量，最后得出各個節點的電壓幅值與角度。