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

142 Programming for Chemical Engineers

x

b

= 0.15

q = 1

R

D

= 3.5

Answer: y-intercept is 0.47, therefore, R

Dmin

= 0.81.

Program Listing 4.5: Distillation program.

1 #include <stdio.h>

2 #include <stdlib.h>

3 #include “benz_tol.h”

4 int main()

5 {

6 double xd=0,xf=0,xb=0,xcoor[100],ycoor[100],q,ycoorf

,xcoorr;

7 double ycoorr,mr,mb,rmin,bmin,bact,dcoorr,ref,delm,

delb;

8 int j,n=0,fplate;

9 printf(“\nThis program will solve problems on

Distillation of Benzene-Toluene\n”);

10 printf(“\nEnter fraction of benzene in distillate : “);

11 scanf(“%lf”,&xd);

12 printf(“Enter fraction of benzene in feed : “);

FIGURE 4.13: McCabe-Thiele method used

to determine ...