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

516 Appendix

function A_inv = special_function_A_inv (x, y, N)
N1 =N+1;
A_inv = zeros(N1, N1);
Fn = special_function_Fn(x, y, 2 ∗ N + 1);
am = -1;
for m=0: N
am =-am ∗ (m + 1) / x;
for n=0:N
s=0;
for k=0:m
if (k == 0)
ak = 1;
else
ak =-ak ∗ x/ k;
end
ss = 0;
for j=0:m-k
ss = ss + nchoosek(m - k, j) ∗ y ^ (n + j + 1) ...
/ factorial(n +j+1);
end
s=s+ak ∗ ss - nchoosek(m, k) ∗ Fn(n +k+1);
end
A_inv(m + 1, n + 1) = am ∗ s;
end
end
end

function A_id = special_function_A_id (x, y, N)
N1 =N+1;
A_id = zeros(N1, N1);
Fn = special_function_Fn(x, y, 2 ∗ N + 1);
am = 1;
for m=0: N
am = am ∗ (m + 1) / x;
for n=0:N
s=0;
for k=0:m
ss = 0;
for j=0:k
ss = ss + nchoosek(k, j) ∗ Fn(n +j+1);