// chapter 8 , Example 8.4 , pg 239
H=1800  //magnetizing field  (in  A/m)
phi=3*10^-5   //magnetic  flux   (in Wb)
A=0.2 *10^-4  //area  (in m^2)
B=phi/A    //magnetic flux density   (in Wb/m^2)
u=B/H    //permeability    (in H/m)
printf("permeability (in  H/m )=")
disp(u)