.title ch3_1.sp v1 1 0 1 pwl(0 0 .01ms 0 .02ms 1) r1 1 2 1 a1 (2 0) (3 0) induct1 .model induct1 lcouple(num_turns=100) a2 (3 0) iron_core .model iron_core core (H_array = [-2000 -1000 0 1000 2000] B_array = [-.1 -.1 0 .1 .1] area = 0.00008 length = 0.04) .save all @r1[i] .control tran .01m 2m plot @r1[i] .endc .end