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

Appendix 515

Gn_2 = Gn_1;
Gn_1 = Gn;
end
end

function mu = special_function_mu (x, y)
n = 500;
mu = ((special_function_F0(y, 0)) ^ 2 ...
+ (special_function_F0(y, x)) ^ 2) / 2;
for i=1:n-1
xx = i ∗ x/n;
mu = mu + (special_function_F0(y, xx)) ^ 2;
end
mu = (special_function_F1(y, x) - mu ∗ x/n)/y;
end

function mu1 = special_function_mu1 (x, y, z, phi)
n = 500;
mu1 = (special_function_F0(z, phi ∗ x) ∗ special_function_F0(y, x) ...
+ special_function_F0(z, 0) ∗ special_function_F0(y, 0)) / 2;
for i=1:n-1
xx = i ∗ x/ n;
mu1 = mu1 + special_function_F0(z, phi ∗ (x - xx)) ...
∗ special_function_F0(y, (x - xx));
end
mu1 = mu1 / n;
end

function mu2 = special_function_mu2 (x, y, z, phi)
n = 500;
mu2 = (special_function_F0(z, 0) ∗ special_function_F0(y, x) ...
+ special_function_F0(z, phi ∗ x) ∗ special_function_F0(y, 0)) / 2;
for i=1:n-1
xx = i ∗ x/ n;
mu2 = mu2 + special_function_F0(z, phi ∗ xx) ...
∗ special_function_F0(y, (x - xx));
end
mu2 = mu2 / n;
end