d(Nmeoh)/d(t)=rhac*V+3-Fmeoh
d(V)/d(t)=vo-((Fmeac*74.08)+(Fmeoh*32.04)+(Fhac*60.05)+(Fh2o*18))/943.4
d(Nhac)/d(t)=rhac*V-Fhac
d(Nmeac)/d(t)=-rhac*V-Fmeac
d(Nh2o)/d(t)=-rhac*V-Fh2o
d(I)/d(t)=Fhac
vo=3/5
X=(300-Nhac-I)/300
Cmeoh=Nmeoh/V
Chac=Nhac/V
Cmeac=Nmeac/V
Ch2o=Nh2o/V
T=325
Ptotal=101.3
Ftotal=100
Nt=Nmeoh+Nhac+Nmeac+Nh2o
Xmeac=Nmeac/Nt
k=(8.88*10**8)*exp(-7032.1/T)
K=5.2*exp((-8000/1.978)*((1/298)-(1/T)))
Xmeoh=Nmeoh/(Nt)
Xhac=Nhac/Nt
Xh2o=Nh2o/Nt
Pvmeoh=8092*exp(14.413-(14.8248/(T/512.6))-8.8032*(ln(T/512.6))+0.4118*(T/512.6)**6)
Pvhac=5786*exp(12.446-(12.8016/(T/592.7))-7.1135*(ln(T/592.7))+0.3556*(T/592.7)**6)
Pvh2o=22060*exp(11.06-(11.376/(T/647.1))-5.9233*(ln(T/647.1))+0.316*(T/647.1)**6)
Pvmeac=4750*exp(10.703-(11.0088/(T/506.5))-5.4361*(ln(T/506.5))+0.3058*(T/506.5)**6)
Fmeoh=Xmeoh*(Pvmeoh/Ptotal)*Ftotal
Fmeac=Xmeac*(Pvmeac/Ptotal)*Ftotal
Fhac=Xhac*(Pvhac/Ptotal)*Ftotal
Fh2o=Xh2o*(Pvh2o/Ptotal)*Ftotal
rhac=-k*((Cmeoh*Chac)-((Cmeac*Ch2o)/K))
t(0)=1e-08
Nmeoh(0)=0
V(0)=150
Nhac(0)=300
Nmeac(0)=0
Nh2o(0)=0
I(0)=0
t(f)=38