// Example 19.11, page no-548
clear
clc
awt=157.25//atomic weight
an=64//atomic number
d=7860//density
k=9.27*10^-24
avg=6.023*10^26
N=d*8*k*avg/awt
printf("The saturation magnetisation of gadolinium is %.2f*10^6 A/m",N*10^-6)