module magnetic_leakage(t1,t2);
inout t1,t2;
magnetic t1,t2;
 parameter real area=1 from (0:inf);
 parameter real length=1 from (0:inf);
 parameter real leakage=0 from [0:1];
 parameter real hc=1 from (0:inf);
 parameter real bsat=1 from (0:inf);

 real permeance;

 analog
 begin
   @(initial_step)
     permeance = leakage*area*bsat/(length*hc);

   Phi(t1,t2) <+ -permeance * MMF(t1,t2);
 end
endmodule