Problem D11a3 - Approach 2 - Differential Equations for Binary Batch Distillation  |POLVER05_0
d(L)/d(x2) =  L/(k2*x2-x2)  #
d(T)/d(x2) =  (k2-k1)/(ln(10)*(x1*k1*(-1211.033)/(220.79+T)^2+x2*k2*(-1344.8)/(219.482+T)^2))  #
k2 =  10^(6.95464-1344.8/(T+219.482))/(760*1.2)  #
k1 =  10^(6.90565-1211.033/(T+220.79))/(760*1.2)  #
x1 =  1-x2  #
x2(0)=0.4
L(0)=100
T(0)=95.5851
x2(f)=0.8