.. .
z
for (j = 1; j <= n; j ++)
if (nobt[j] == -1)
b[j] = 0.0;
i1 = i - 1;
for (im = 1; im <= i1; im ++) {
jj1 = nfd[im]; jj2 = nfd[im + 1] - 1; for (jj = jj1; jj <= jj2; jj ++) { }
if(iu[jj] == i) { }
ai = u[jj] / di[im]; for(k = jj; k <= jj2; k ++) { } break;
j = iu[k];
b[j] = b[j] - ai * u[k];
}
.. .
z
}
x = nfd[i]; di[i] = 1.0 / b[i]; ai = di[i]; k = 0; i1 = i + 1;
for (j = i1; j <= n; j ++) { }
if (fabs(b[j]) > 1.0e-15) { }
u[x] = b[j] * ai; iu[x] = j; k++; x++;
fd[i] = k;
else { }
fd[i] = 0; di[i] = 0.0;
.. .
}
}
nfd[i+1] = nfd[i] + fd[i];
fprintf(fp2, \for (i = 1; i <= x; i ++) { }
fprintf(fp2, \printf1(di, n);
if(i % 3 == 1)
fprintf(fp2, \
fprintf(fp2, \
void printf1(double* aa, int n) {
//**** 本函数输出aa[i],i=1-n ****//
int i;
for (i = 1; i <= n; i ++)
{
z
if(i % 5 == 1)
fprintf(fp2, \
.. .
fprintf(fp2, \
} }
void calc(int* iu, double* u, double* di, int* nfd, double* b) {
//**** 本函数利用因子表解线形程组。(P417图F1-9) ****//
double bi; int i, j, k, i1, i2;
for (i = 1; i <= n; i ++) //前代过程。 fprintf(fp2, \
{
z
bi = b[i]; i1 = nfd[i]; i2 = nfd[i + 1]; for (j = i1; j < i2; j ++) { }
k = iu[j];
b[k] = b[k] - bi * u[j];
.. .
}
b[i] = bi * di[i];
for (i = n; i >= 1; i --) // 回代过程。
{
bi = b[i]; i1 = nfd[i]; i2 = nfd[i + 1] - 1; for (j = i2; j >= i1; j --) { } b[i] = bi;
k = iu[j];
bi = bi - b[k] * u[j];
} }
void zlsort(int* nnew) ////。 {
//**** 本函数进行支路数据排序
****//
//小节点号放左边，大节点号放右边// //左右皆按从小到大顺序排列//
z
电气工程及其自动化专业综合实训一[新]
