%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%  PM610.m  -  Determines Epstein loss given value of B
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function y = PM610(Bx)
% B-Wlb values are valid for M-6, 10 deg, 29 ga,60 Hz
B=[0 38.7 51.6 64.5 96.75 120 140 167.7]*1000;
Wlb=[0 0.137 0.222 0.325 0.725 1.5 2.5 6.7 ];
% Activate to plot B-Wlb curve
% m=8; plot(Wlb(1:m),B(1:m)); grid; % Linear plot
% m=8; semilogx(Wlb(2:m),B(2:m)); grid;  % Semilog plot
n=length(B); k=0;
if Bx==0; k=-1; y=0; end
if Bx<0; k=-1; y=0; disp('WARNING - Bx < 0, Wlb = 0 returned'); end
if Bx>B(n); y=Wlb(n); k=-1; disp('CAUTION - Beyond B-Wlb curve');end
for i=1:n
if k==0 & (Bx-B(i))<=0; k=i; break; end
end
if k>0;
y=Wlb(k-1)+(Bx-B(k-1))/(B(k)-B(k-1))*(Wlb(k)-Wlb(k-1));
else;
end