function y = table1(tab,x0)
%TABLE1	Table look-up.
%	Y = TABLE1(TAB,X0) returns a table of linearly interpolated rows from
%	table TAB, looking up X0 in the first column of TAB.
%
%	See also INTERP2, TABLE2.

%	NOTE:  TAB's 1st column is checked for monotonicity.
%	Different algorithms are used depending on whether X0 is scalar
%	so that vectorization can be exploited in the case of vector X0.
%	It is an error to request a value outside the range of the first
%	column of TAB for X0.

%	Tomas Schoenthal 5-1-85
%	Egbert Kankeleit 1-15-87
%	Revised by L. Shure 2-3-87
%	Copyright (c) 1984-93 by The MathWorks, Inc.

if (nargin ~= 2), error('Wrong number of input arguments.'), end

[m,n]=size(tab)  ;  k0=max(size(x0))  ;

if k0==1
   dx=tab(2:m,1)-tab(1:m-1,1);
   sig=sign(dx(1,1));
   if any(sign(dx(:))-sig);
   error('First column of the table must be monotonic.')
   end
   if sig > 0
        ii = find(tab(:,1) >= x0);
        if size(ii)==0; error('x0 larger than all values in first column'),end
        if x0 < tab(1,1)
           error('x0 smaller than all values in first column') 
        end
   else
   ii = find(tab(:,1) <= x0);
        if size(ii)==0;
           error('x0 smaller than all values in first column')
        end
        if x0 > tab(1,1)
           error('x0 larger than all values in first column')
        end
   end
   kk=ii(1);
   xhigh = tab(kk,1);
   if kk > 1 
      xlow = tab(kk-1,1);
   else
      xlow = xhigh;
   end
   if xlow == xhigh
      y = tab(kk,:);
   else
      alpha = (xhigh-x0)/(xhigh-xlow);
      beta = 1. - alpha;
      y = alpha*tab(kk-1,:) + beta*tab(kk,:);
   end
   y = y(2:n);
   return
else
   dx=tab(2:m,:)-tab(1:m-1,:);
   dx=[dx;dx(m-1,:)];sig=sign(dx(1,1));
   if any(sign(dx(:,1))-sig);
   error('First column of the table must be monotonic.')
   end
   y=zeros(k0,n-1);
   for k=1:k0
     if sig>0; ii=max(find(tab(:,1) <= x0(k) ));
        if size(ii)==0; error('x0 smaller than all values in first column'),end
        if x0(k) > tab(m,1)
           error('x0 larger than all values in first column') 
        end
     else ; ii=max(find(tab(:,1) >= x0(k) ));
        if size(ii)==0;
           error('x0 larger than all values in first column')
        end
        if x0(k) < tab(m,1)
           error('x0 smaller than all values in first column')
        end
     end
     y(k,:)=tab(ii,2:n)+dx(ii,2:n)*(x0(k)-tab(ii,1))/dx(ii,1);
   end
end
