# Parametre
coil_group = 7
filename = 'example.fem'


# Inputgreier
strom = -1.0;          % Stromstyrken din

Xmax = 5.0             % Største utsving i positiv retning
Xmin = -5.0            % Største utsving i negativ retning
Xinc = 0.5             % Størrelsen på inkrementet
Xpos = Xmin:Xinc:Xmax  % Liste med alle posisjoner som simuleres

# Tom liste for a putte ting i...
bl = [];               % Tom liste for a lagre Bl verdier i...

% open FEMM and initialize the problem
openfemm;
try
    vv=ver; opendocument([cd,strjoin('/',filename)]);
catch
    opendocument(filename);
end
mi_saveas('temp.fem');

# Set stromstyrken til noe
mi_setcurrent('Prove circuit',strom);

mi_selectgroup(7);
mi_movetranslate(0, Xmin);
mi_clearselected;

for i=Xpos,
  mi_selectgroup(coil_group);
  mi_movetranslate(0, Xinc);
  mi_clearselected;

  # Kjor analysen
  mi_analyze;
  # Last opp resultat
  mi_loadsolution;
  # Ta ut en liste med noe greier som har me spolen a gjore
  mo_groupselectblock(coil_group);
  fz = mo_blockintegral(12);
  bl = [bl; [i, fz]];
  mo_close;
  #mo_clearselected;


end
disp(bl)
# Lag pizza
#closefemm

plot(bl(:,1),bl(:,2));
xlabel('Displacement, mm');
ylabel('BL, N/A');
title('BL vs. Displacement');