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

514 Appendix

F_1 = exp(- x - y) ∗ s;
end
end
function G_1 = special_function_G_1 (x, y)
if (x == 0)
G_1 = exp(- y);
elseif (y == 0)
G_1 = exp(- x);
else
m=1;
s=1;
a=1;
while (abs(a) > abs(s) ∗ 1E-8)
a=a ∗ x ∗ y/ m ^ 2;
s=s+a;
m=m+1;
end
G_1 = exp(- x - y) ∗ s;
end
end

function h = special_function_h (x, y, z)
h= 0;
a= 1;
Gn_3 = special_function_F_1(x, y) - special_function_G_1(x, y); % G_2
Gn_2 = special_function_G_1(x, y);
Gn_1 = special_function_G0(x, y);
for n=1: 10000
a=a ∗ z;
Gn =((y -x-2 ∗ n+1) ∗ Gn_1 + (2 ∗ y-n +1) ∗ Gn_2 ...
+y ∗ Gn_3) / n;
Fn = Gn + Gn_1;
as = a ∗ Fn;
h=h +as;
if (abs(as) < abs(h) ∗ 1E-8 && n >= 10)
break;
end
Gn_3 = Gn_2;