% Bruk: z = fmins('syntesegass1',[70 20])
% Dette er modellen for likevekt for syntesegass for eksemplet i Kap.3
%
function z = f(x)
x1 = x(1); x2 = x(2);
% føde (etter at de tunge hydrokarbonene har reagert)
   n0ch4 = 95.5; n0co = 9.1; n0h2 =  22.1; n0h2o = 240.9; n0n2 = 0.6; n0co2=0;
% massebalanser
   nch4 = n0ch4 - x1 
   nh2o = n0h2o - x1 - x2
   nco  = n0co  + x1 - x2
   nh2  = n0h2  + 3*x1 + x2
   nco2 = n0co2 + x2
   nn2  = n0n2;
   ntot = nch4+nh2o+nco+nh2+nco2+nn2
% Likevekt
   p=21;  K1=710; K2=0.81;
   Q1a= (nco * nh2^3) / (nch4*nh2o);
   Q1 = Q1a * p^2/ntot^2
   Q2 = (nco2* nh2) / (nco * nh2o)
% 
z = (Q1/K1-1)^2 + (Q2/K2-1)^2


