3557/CH18/EX18.1/Ex18_1.sce


1
2
3
4
5
6
7
8
9

//Example 18.1//
ur=1.01;
u0=4*%pi*10^-7;//henry/m
H=2*10^5;//amperes/m
B=ur*u0*H
mprintf("B = %f weber/m^2",B)
//Using second equality, we obtain
M=(ur-1)*(H)
mprintf("\nM = %e amperes/m",M)