d(Nmeoh)/d(t)=rhac*V+Fmeoho
d(V)/d(t)=vo-Fmeac*74/943.4
d(Nhac)/d(t)=rhac*V
d(Nmeac)/d(t)=-rhac*V-Fmeac
d(Nh2o)/d(t)=-rhac*V
Fmeoho=3
Cmeoh=Nmeoh/V
Chac=Nhac/V
Cmeac=Nmeac/V
Ch2o=Nh2o/V
T=350
Cmeoho=5
Xmeac=Nmeac/(Nmeoh+Nhac+Nmeac+Nh2o)
Ptotal=101.3
Ftotal=100
k=(8.88*10**8)*exp(-7032.1/T)
Tr=T/506.5
K=5.2*exp((-8000/1.978)*((1/298)-(1/T)))
vo=Fmeoho/Cmeoho
Pvmeac=4750*exp(10.703-(11.0088/Tr)-5.4361*(ln(Tr))+0.3058*Tr**6)
Fmeac=Xmeac*(Pvmeac/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
t(f)=120