showconsole()
mydir="./"
open(mydir .. "mcm_2coil.fem")

mi_saveas(mydir .. "temp.fem")
clearconsole()

mi_analyze()
mi_loadsolution()

for n=0,4,0.1 do

  A, B1, B2 =mo_getpointvalues(n,10.1)
  bbot=sqrt(B1*B1 + B2*B2)
  A, B1, B2 =mo_getpointvalues(n,11.4)
  btop=sqrt(B1*B1 + B2*B2)
  print(n,",",bbot,",",btop)
  
  end

mo_close()
mi_close()
