//Example 18.5//

a=18.4;// measured value of nickel ferrite
ub=9.274*10^-24;//A m^2// ampere-meters square //Moment
v=(0.833*10^-9);//m //meter // volume of unit cell
Ms=(a*ub)/v^3
mprintf("Ms = %e A/m",Ms)