BATCH DISTILLATION :part 1( equilibrium eq: y=k*x)
------------------
------------------

part1 : Integral batch distillation.  
peter.vlasschaert@gmail.com , 18/02/2017

kill(all)$

eq1:integrate(x^2, x);

eq1:integrate(x^2,x,a,b);

assume(x >= 0 and x <= 1);
assume(y >= 0 and y <= 1);
assume( b>a);
assume( a >0 and a < 1 and b < 1);

eq2:y=k*x;

eq3:op1='integrate(1/(rhs(eq2)-x),x,a,b);

eq4:combine(ev(eq3,nouns));

eq5:radcan((log(b*abs(k-1))-log(a*abs(k-1)))/(k-1));

eq6:op1=logcontract(eq5);

example :

k= 2.3
a= 0.4
b= 0.7

eq7:ev(eq6,k=2.3,a=0.4,b=0.7);