function h=pdesurf(p,t,u)
%PDESURF Surface plot of PDE solution.
%
%       PDESURF(P,T,U) plots a 3-D surface of PDE node or triangle data.
%       If U is a column vector, node data is assumed, and continuous
%       style and interpolated shading are used.  If U is a row vector,
%       triangle data is assumed, and discontinuous style and flat shading
%       are used.
%
%       H=PDESURF(P,T,U) additionally returns a handle to the patch objects.
%
%       See also PDEPLOT

%       A. Nordmark 4-26-94, LL 1-30-95.
%       Copyright 1994-2003 The MathWorks, Inc.
%       $Revision: 1.8.4.1 $  $Date: 2003/11/18 03:12:02 $

if size(t,2)==size(u,2)
  cont='discontinuous';
  shad='flat';
elseif size(p,2)==size(u,1)
  cont='continuous';
  shad='interp';
else
  error('PDE:pdesurf:InvalidSolFormat', 'Illegal solution format.');
end

hh=pdeplot(p,[],t,'xydata',u,'xystyle',shad,...
           'zdata',u,'zstyle',cont,'colorbar','off');

if nargout==1
  h=hh;
end