//牛拉法解潮流程序

#include
#include
#define?N?5??????????????//節點數
#define?n_PQ?3?????????//PQ節點數
#define?n_PV?1?????????//PV節點數
#define?n_br?5?????????//串聯支路數
void?main（）
{
void?disp_matrix（float?*disp_pint?disp_mint?disp_n）;???????//矩陣顯示函數
float?Us[2*N]={1.001.001.001.0501.050};?????????????????????//電壓初值
float?Ps[N]={-1.6-2-3.75};??????????????????????????????????????//有功初值
float?Qs[N]={-0.8-0-1-1.3};?????????????????????????????????????????????????//無功初值????????
float?G[N][N]B[N][N];?????????????????????????????????????????????????//各幾點電導電納
struct??????????????????????????????????????????????????????????????????????//阻抗參數
???{
????int?nl;?????????????//左節點
????int?nr;????????????//右節點
????float?R;??????????//串聯電阻值
????float?X;???????