%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%  Pc.m  -  Determines Epstein loss given value of B
%           Used by imperf.m
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function y = Pc(Bx)
% B-W/lb values that follow are valid for M-27,24 ga,60 Hz
B=[0 10 19 24 30 35 40 45 50 55 60 64 70 75 80 85 90 ...
95 100 105 110 115 120 125 130 135 138 165 210 ]*1000;
Wlb=[0 0.08 0.15 0.195 0.27 0.35 0.43 0.52 0.625 0.76 ...
0.89 1.0 1.2 1.35 1.55 1.75 2.0 2.25 2.55 2.8 3.15 3.5 ...
3.85 4.2 4.7 5.2 5.6 11.0 30.0 ];
% Activate to plot B-W/lb curve
% m=29; plot(Wlb(1:m),B(1:m)); grid; pause; % Linear plot
% m=29; semilogx(Wlb(2:m),B(2:m)); grid; pause; % 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