Page 538 - Design and Operation of Heat Exchangers and their Networks

Appendix 521

ay = 1;
ex =-x-y;
while (abs(sx) > abs(sy) ∗ 1E-8)
j=j+ 1;
ay = ay ∗ y/j;
if (abs(ay) > CAS_EXP_350)
ay = ay / CAS_EXP_350;
sy = sy / CAS_EXP_350;
ex = ex + 350;
end
ax = ay;
sx = ay;
for k=1:1: j-1
ax = ax ∗ x/ k;
if (abs(ax) > CAS_EXP_350)
ax = ax / CAS_EXP_350;
sx = sx / CAS_EXP_350;
ay = ay / CAS_EXP_350;
sy = sy / CAS_EXP_350;
ex = ex + 350;
end
sx = sx + ax;
end
sy = sy + sx;
end
while (ex > 350)
sy = sy ∗ CAS_EXP_350;
ex = ex - 350;
end
while (ex < -350)
sy = sy / CAS_EXP_350;
ex = ex + 350;
end
t_c = exp(ex) ∗ sy;
end
end