function [W,resN]=femFreqSweep(fem,lambdaRange,intSDNr)

sfem=fem;

N=size(lambdaRange,2);

% for n=1:N
%     disp(num2str(lambdaRange(n),'%14.13e'));
% end

sfem=SolveFemStaticParaSweep(sfem,'lambda0',lambdaRange);
              
W=postint(sfem,'Weav_rfwh','unit','N','dl',intSDNr,'solnum',1:N);

flclear sfem;
resN=find(W==max(W));