Page 485 - Advanced Mine Ventilation

454 Appendix C: Ventilation Network Analyzer in Cþþ With Input and Output

for(int i¼1; i<¼n1; iþþ)
{
k ¼ iþ1;
0;
z[i] ¼
for(int j¼1; j<¼np; jþþ)
{
z[i] ¼ z[i] þ fx[j] * p[k][j];
}
s[i] ¼ z[i]/w[i];
dsq ¼ dsq - s[i] * s[i] * w[i];
sigma2[i] ¼ dsq/(float)(np-i);

if(sigma2[i] - 1e-6 > 0)
{
if (i-n1 < 0)
{
alpha[i] ¼ 0;

for(int j¼1; j<¼np; jþþ)
{
alpha[i] þ x[j] * p[k][j] * p
alpha[i] ¼
[k][j];
}
alpha[i] ¼ alpha[i]/w[i];
w[k] ¼ 0;

for(int j¼1; j<¼np; jþþ)
{
p[kþ1][j] ¼ (x[j] - alpha[i]) * p[k][j] -
beta[i] * p[i][j];
w[k] ¼ w[k] þ p[kþ1][j] * p[kþ1][j];
}
beta[k] ¼ w[k]/w[i];
}
else
{
break;
}
}
else
{
goto label78;
}
}