% Dette er: nh3eq.m
% Bruk:  x = fzero('nh3eq',x0)  der x0 er startverdien (f.eks. x0=0.1)
%
% Beskrivelse: nh3eq beregner verdien av funksjonen z=f(x), der
%    x - reaksjonsomfang  
%    z = Q-K - avvik fra likevekt
%
function z=f(x)   
p=200;            % trykket er her 200 bar (kan endres)
nn2=1-x;
nh2=3-3*x;
nnh3=2*x;
ntot=nn2+nh2+nnh3;
xn2=nn2/ntot
xh2=nh2/ntot
xnh3=nnh3/ntot
pn2=p*xn2;
ph2=p*xh2;
pnh3=p*xnh3;
Q = (pnh3^2 / (pn2 * ph2^3));
K=6.36e-5;    % K er her ved T= 700 K (kan endres)
z=Q-K;

