function n = N_A()

% N_A returns the Avogadro constant 6.022 141 79(30) x 10^23 mol^(-1) 

n = 6.02214179e23;
return