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

512 Appendix

function nu_ = special_function_nu_ (x, y)
nu_ = exp(- x ∗ special_function_K(y));
end
function nu = special_function_nu (x, y)
nu = special_function_F1(x, y) / y;
end

function F1 = special_function_F1 (x, y)
F1 = special_function_G0(x, y) + special_function_G1(x, y);
end

function G1 = special_function_G1 (x, y)
G1 = (y - x - 1) ∗ special_function_G0(x, y) ...
+y ∗ (special_function_F_1(x, y) + special_function_G_1(x, y));
end

function F0 = special_function_F0 (x, y)
^
% F0(x,y)=exp(-x-y)∗sum(m = 0:Inf;k = 0:m)(y m/m!)∗(x k/k!)
^
if (x == 0)
F0 = 1;
elseif (y == 0)
F0 = exp(- x);
else
m=1;
sx = 1;
sy = 1;
ay = 1;
while (abs(sx) > abs(sy) ∗ 1E-8)
ay = ay ∗ y/m; % y m/m!
^
ax = ay;
sx = ay;
for k=1 :m
ax = ax ∗ x/k;%(y m/m!)∗(x k/k!)
^
^
sx = sx + ax;
end
sy = sy + sx;
m=m+1;
end
F0 = exp(- x - y) ∗ sy;
end
end