
xcoorr= compute(xf,1);
do
{
xcoorr=xcoorr+0.00000005;
ycoorf=(((-1*q)/(1-q))*xcoorr)+(xf/(1-q));
ycoorr= compute(xcoorr,2);
dcoorr=ycoorf-ycoorr;
}while(dcoorr > 0.000001 || dcoorr< -0.000001);
}
mr=(xd-ycoorr)/(xd-xcoorr);
bmin=ycoorr-(mr*xcoorr);
rmin=(xd/bmin)-1;
cout<<»Minimum reflux is = «<<rmin<<endl;
}
void DISTILL::actualPlates(double ref){
double xcoorr,delm,delb,ycoorr,mb,bact,ycoor[100],xcoor
[100];
int j,fplate,n;
if(q==1)
(xcoorr=xf);
else
{
delm=(xd/(ref+1))-(xf/(1-q));
delb=((-1*q)/(1-q))-(ref/(ref+1));
xcoorr=delm/delb;
}
ycoorr=((ref/(ref+1))*xcoorr)+(xd/(ref+1));
mb=(ycoorr-xb)/(xcoorr-xb);
bact=ycoorr-(mb*xcoorr);
j=0;ycoor[0]=xd; ...