function values = tspval(tsp,x,y)
%TSPVAL	Evaluate a tensor product spline.
%
%        values = tspval(tsp,x,y)
%
% returns the array  values(i,j) := f(x(i),y(j)) , with  f  the tensor product 
% spline  f  whose B-form is contained in  tsp .

% C. de Boor: Feb.4, 1991
% Copyright (c) 1990-92 by Carl de Boor and The MathWorks, Inc.

[knotsx,knotsy,coefs]=tspbrk(tsp);
values=fnval(spmak(knotsx,fnval(spmak(knotsy,coefs,1),y)',1),x)';
