%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%    H1010.m  -  B-H interpolation routine
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function y = H1010(Bx)
% B-H values that follow are valid for 1010 steel plate
B=[0 2.6 5.2 7.7 10.3 12.9 18.1 24.5 31 38.7 51.6 64.5 ...
      71 77 83.8 90.3 97 103 110 116 122 129 135 142 ...
      150 260]*1000;     % kilolines/sq_in
H=[0 0.6 1.3 2 2.5 2.8 3.4 4 4.7 5.5 6.9 8.5 9.9 11.9 ...
      14.3 19.2 28.3 46.5 86.9 155.6 242.5 444.6 647 ...
      950 2021 4e4];   % A_t/in
% Activate to plot B-H curve
% m=26; plot(H(1:m),B(1:m)); grid; % Linear plot
% m=26; semilogx(H(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, Hx = 0 returned'); end
if Bx>B(n); y=H(n); k=-1; disp('CAUTION - Beyond B-H curve'); end
for i=1:n
if k==0 & (Bx-B(i))<=0; k=i; break; end
end
if k>0;
y=H(k-1)+(Bx-B(k-1))/(B(k)-B(k-1))*(H(k)-H(k-1));
else;
end