电力系统潮流计算 迭代算法

#include <iostream.h>
#include <fstream.h>
#include <math.h>
#include <ctype.h>
#include <string.h>
#include <stdlib.h>
void input();
double num[3];       //结点个数，线路个数，变压器个数
double jd1[5];       //节点1数据
double jd2[5];       //节点2数据
double jd3[5];       //节点3数据
double xl1[5];       //线路1数据
double byq[7];      //变压器数据
double jd ;

void input()        //迭代精度
{double a[30];
ifstream infile("data.txt",ios::in);
for(int i=0;i<31;i++)
infile>>a[i];

for(i=0;i<3;i++)
                   num[i]=a[i];
for(i=0;i<5;i++)
                   jd1[i]=a[i+3];
for(i=0;i<5;i++)
                   jd2[i]=a[i+8];
for(i=0;i<5;i++)
                   jd3[i]=a[i+13];
for(i=0;i<5;i++)
                   xl1[i]=a[i+18];
for(i=0;i<7;i++)
                   byq[i]=a[i+23];
                         jd=a[30];
}
void main()
{
input();
double str,stm,star,stam,ut1,ut2,ul1,ul2,a0=0,a1=0,a2=0,a3=0,QB1,QB2,slr,slm,ut12,ut22,U[3]