function [a,b,c,e]=FemLorentzFit(x,y,ResLambda,lambdaBW,femHandle,sweepHandle,args)

e=0;

fem1=femHandle(x,y,args{:},0);

F(:,1)=[ResLambda-lambdaBW/2,ResLambda,ResLambda+lambdaBW/2];

[F(:,2) N]=sweepHandle(fem1,F(:,1)',5); %subdomain 5 in this case only!!!!

[a,b,c]=LorentzFit(F);

if (b>max(F(:,1)) || b<min(F(:,1)))
    disp('error: Fittet Resonance out of bounds')
    disp([num2str(F(1,2),'%10.9e\n') ', ' num2str(F(2,2),'%10.9e\n') ', ' num2str(F(3,2),'%10.9e\n')]);
    e=1;
    return;
end