
d(Nmeoh)/d(t)=rhac*V+Fmeoho
d(V)/d(t)=vo
d(Nhac)/d(t)=rhac*V
d(Nmeac)/d(t)=-rhac*V
d(Nh2o)/d(t)=-rhac*V
Fmeoho=3
Cmeoh=Nmeoh/V
Chac=Nhac/V
Cmeac=Nmeac/V
Ch2o=Nh2o/V
T=350
Cmeoho=5
k=(8.88*10**8)*exp(-7032.1/T)
K=5.2*exp((-8000/1.978)*((1/298)-(1/T)))
vo=Fmeoho/Cmeoho
rhac=-k*((Cmeoh*Chac)-((Cmeac*Ch2o)/K))
t(0)=0
Nmeoh(0)=0
V(0)=150
Nhac(0)=300
Nmeac(0)=0
Nh2o(0)=0
t(f)=120
