
396 Programming for Chemical Engineers
s = 0.0;
for (i = 0; i < 7; i++)
{
s = s + p[i].ni * p[i].di * p[i].ti * powi (delta,
p[i].di - 1)
* pow (tau, p[i].ti - 1);
}
for (i = 0; i < 44; i++)
{
dc = powi (delta, q[i].ci);
s = s + q[i].ni * q[i].ti * powi (delta, q[i].di - 1)
* powi (tau, q[i].ti - 1) * (q[i].di - q[i].ci *
dc)* exp (-dc);
}
for (i = 0; i < 3; i++)
{
de = delta - z[i].epsiloni;
te = tau - z[i].gammai;
ex = -z[i].alphai * de * de - z[i].betai * te * te;
dd = powi (delta, z[i].di);
s = s + z[i].ni * dd * powi (tau, z[i].ti) *
exp(ex)
* ((z[i].di / delta) - 2.0 * z[i].alphai * de)
*