//determine the value of relative permeability of iron
D=15e-2
l=%pi*15e-2
N=450
I=2
B=1.2
u=B/(4*%pi*10^-7*N*I*l)
disp('value of relative permeability='+string(u)+'  ')