// Winding for Ferromagnetic Core
module magnetic_winding(e1,e2,m1,m2);
 inout e1,e2,m1,m2;
 electrical e1,e2;
 magnetic m1, m2;
 parameter real turns=1 from (0:inf);
 parameter real r=0 from [0:inf);

 analog
 begin
   V(e1,e2)   <+ turns * ddt( Phi(m1,m2) ) + r * I(e1,e2);
   MMF(m1,m2) <+ turns * I(e1,e2);
 end
endmodule