//Variable declaration
chi=-0.42*10**-3;      //magnetic susceptibility
H=1000;      //magnetic field(A/m)
mew0=4*%pi*10**-7;    

//Calculation
M=chi*H;      //magnetisation(A/m)
B=mew0*(H+M);     //flux density(Wb/m**2)

//Result
printf('magnetisation is %0.3f    A/m \n',M)
printf('flux density is %0.3f   *10**-3 Wb/m**2 \n',(B*10**3))