m = 9.11*10^-31; //mass of electron ml = 0.98*m; ms = 0.19*m; mn = 6^(2/3)*(ml*ms*ms)^(1/3); //density of states effective mass calculation mn0 = mn/m; disp(mn,"density of states effective mass (in kilogram)=") disp(mn0,"density of states effective mass in proportion to mass of electron= ")