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

Appendix 517

end
s=s+ nchoosek(m, k) ∗ y ^ (n +k+1) ...
/ factorial(n +k+1) ...
-x ^ (m - k) / factorial(m - k) ∗ ss;
end
A_id(m + 1, n + 1) = am ∗ s;
end
end
end
function B = special_function_B (x, N)
N1 =N+1;
B = zeros(N1, N1);
for n=0:N
if (n == 0)
an = 1;
else
an = an ∗ x/ n;
end
for m=1:N1
B(m, n + 1) = an ∗ m/(n+m);
end
end
end

function Fn = special_function_Fn (x, y, N) % Fn(x,y), n = 1:N
Fn = zeros(N);
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:N
Gn =((y -x-2 ∗ n+1) ∗ Gn_1 + (2 ∗ y-n+1) ∗ Gn_2 ...
+y ∗ Gn_3) / n;
Fn(n) = Gn + Gn_1;
Gn_3 = Gn_2;
Gn_2 = Gn_1;
Gn_1 = Gn;
end
end