function [err,n5ntvalue,n1,n2,n3,n4,n5,n6] = nEquations(n5ntvalue,...
    pressure,fuelequiv,carbon,oxygen,hydrogen,...
    molar_air_stoich,kco2,kh2o)

% This function calculates the values for n, a, and error.


   a = sqrt(n5ntvalue * pressure);
   n3 = (carbon * fuelequiv) / (kco2 * a + 1);
   n4 = ((hydrogen / 2) * fuelequiv) / (kh2o * a + 1);
   n1 = carbon * fuelequiv - n3;
   n2 = (hydrogen / 2) * fuelequiv - n4;
   n5 = (molar_air_stoich * 2 + oxygen * fuelequiv - 2 * n1 - n2 - n3) / 2;
   n6 = molar_air_stoich * 3.76;
   ntsum = n1 + n2 + n3 + n4 + n5 + n6;
   err = n5ntvalue - n5 / ntsum;

end