BATCH DISTILLATION :part 1( equilibrium eq: y=k*x)
------------------
------------------
part1 : Integral batch distillation.
[email protected] , 18/02/2017
(%i1) | kill(all)$ |
(%i1) | eq1:integrate(x^2, x); |
(%i2) | eq1:integrate(x^2,x,a,b); |
(%i6) |
assume(x >= 0 and x <= 1); assume(y >= 0 and y <= 1); assume( b>a); assume( a >0 and a < 1 and b < 1); |
(%i7) | eq2:y=k*x; |
(%i8) | eq3:op1='integrate(1/(rhs(eq2)-x),x,a,b); |
(%i9) | eq4:combine(ev(eq3,nouns)); |
(%i10) | eq5:radcan((log(b*abs(k-1))-log(a*abs(k-1)))/(k-1)); |
(%i11) | eq6:op1=logcontract(eq5); |
example :
k= 2.3
a= 0.4
b= 0.7
(%i12) | eq7:ev(eq6,k=2.3,a=0.4,b=0.7); |