简单极坐标下牛顿拉夫逊法潮流计算

该程序实现简单的牛拉法极坐标下的潮流计算。首先将支路导纳矩阵转化为节点导纳矩阵，初始化功率初值，电压初值，接着形成Jacobian矩阵，对矩阵进行LDU分解，前代回代求解修正量，最后得出各个节点的电压幅值与角度。