//chapter6,Example6_2,pg 121 Jc=1 sig=5.8*10^7 E=(Jc)/sig printf("electric field established\n") disp(E)