clear;
clc;
a=5*10^-10 //lattice constatnt in m

//Calculation
n111=1/(a^2*sqrt(3))

mprintf("n(111)= %.1e atoms/m^2",n111)
//2.3e+18 is 2.3*10^18