电力系统潮流计算注：这是一个基于 N－R 法的潮流计算通用程序，仅提供了子程序，需要做些处理才能成为一个可运行的计算程序
此程序非我原创，仅与大家共享
/***＊*＊＊＊＊*＊**＊＊＊**＊*＊＊＊**＊*＊＊＊＊*＊***＊**＊＊＊****＊＊＊*＊*****＊＊***＊＊＊*＊** 这里提供的是电力系统潮流计算机解法的五个子程序,采纳的方法是 *＊ Newton_Raphson 法
** 程序中所用的变量说明如下： ＊＊ N：网络节点总数
M：网络的 PQ 节点数
*＊ L：网络的支路总数
N0：雅可比矩阵的行数
＊* N1：N0+1 K：打印开关
K=1，则打印；否则,不打印
** K1：子程序 PLSC 中推断输入电压的形式
K1=1,则为极座标形式
否则** 为直角坐标形式
＊* D：有功及无功功率误差的最大值
＊＊ G(I，J):Ybus 的电导元素（实部）
*＊ B(I，J)：Ybus 的电纳元素（虚部）
*＊ G1(I） :第 I 支路的串联电导
B1(I）:第 I 支路的串联电纳
*＊ C1（I) :第 I 支路的 pie 型对称接地电纳
＊＊ C(I,J)：第 I 节点 J 支路不对称接地电纳
＊* CO（I) :第 I 节点的接地电纳
＊＊ S1(I) ：第 I 节点的起始节点号
E1(I)：第 I 节点的终止节点号
*＊ P（I) ：第 I 节点的注入有功功率
Q(I）:第 I 节点的注入无功功率
** P0（I） ：第 I 节点有功功率误差
Q0（I):第 I 节点无功功率误差
** V0（I) :第 I 节点(PV 节点)的电压误差(平方误差）
＊＊ V（I） ：第 I 节点的电压误差幅值
＊* E（I) ：第 I 节点的电压的实部
F(I):第 I 节点的电压的虚部
＊＊ JM（I,J）：J
