Programming for Chemical Engineers Using C, C++, and MATLAB®

60 hd(j+1)=enthalpy_l(temp_t);

61 hd(j+1)=hd(j+1)/(1055.056*2.2046226);

62 end

63 xmatrx = zeros(k);

64 constant(1,1)=f*(hl(1)-hl(2));

65 for j=2:k-1

66 constant(j,1)=f*(hl(j)-hl(j+1));

67 end

68 constant(k,1)=f*(hl(k)-hv(k+1)+(conc(1)*(hv(k+1)-

hl(k+1))/conc(k+1)));

69 xmatrx(k,1)=0;

70 for j=2:k-1

71 xmatrx(k,j)=hl(k)-hv(k+1);

72 end

73 xmatrx(k,k)=hd(k+1)+hl(k)-hv(k)-hv(k+1);

74 xmatrx(1,1)=hd(2)-hv(1);

75 for j=1:k-1

76 xmatrx(j,j+1)=hv(j+1)-hl(j+1);

77 end

78 for j=2:k-1

79 xmatrx(j,j)=hd(j+1)+hl(j)-hv(j)-hl(j+1);

80 end

81 for j=3:k-1

82 for w=2:k-2

83 xmatrx(j,w)=hl(j)+hl(j+1);

84 end

85 end

86 vapor=xmatrx\constant;

87 totv=0;

88 for ...