clear
//
//
//

//Variable declaration
chi=0.5*10**-5    //magnetic susceptibility
H=10**6     //field strength(ampere/m)
mew0=4*%pi*10**-7

//Calculation
I=chi*H     //intensity of magnetisation(ampere/m)
B=mew0*(I+H)    //flux density in material(weber/m**2)

//Result
printf("\n intensity of magnetisation is %0.3f  ampere/m",I)
printf("\n flux density in material is %0.3f  weber/m**2",B)