function hh = ezsurf(varargin)
%EZSURF   Easy to use 3-D colored surface plotter
%   EZSURF(FUN) plots a graph of the function FUN(X,Y) using SURF. FUN is
%   plotted over the default domain -2*PI < X < 2*PI, -2*PI < Y < 2*PI.
% 
%   EZSURF(FUN,DOMAIN) plots FUN over the specified DOMAIN instead of the
%   default domain.  DOMAIN can be the vector [XMIN,XMAX,YMIN,YMAX] or the
%   vector [A,B] (to plot over A < X < B, A < Y < B).
%
%   EZSURF(FUNX,FUNY,FUNZ) plots the parametric surface FUNX(S,T),
%   FUNY(S,T), and FUNZ(S,T) over the square -2*PI < S < 2*PI and
%   -2*PI < T < 2*PI. 
%
%   EZSURF(FUNX,FUNY,FUNZ,[SMIN,SMAX,TMIN,TMAX]) or
%   EZSURF(FUNX,FUNY,FUNZ,[A,B]) uses the specified domain.
%
%   EZSURF(...,N) plots f over the default domain using an N-by-N grid.
%   The default value for N is 60.
%
%   EZSURF(...,'circ') plots f over a disk centered on the domain.
%
%   EZSURF(AX,...) plots into AX instead of GCA.
%
%   H = EZSURF(...) returns a handle to the surface object in H.
%
%   Examples:
%   The easiest way to express a function is via a string:
%      ezsurf('x*exp(-x^2 - y^2)')
%
%   One programming technique is to vectorize the string expression using
%   the array operators .* (TIMES), ./ (RDIVIDE), .\ (LDIVIDE), .^ (POWER).
%   This makes the algorithm more efficient since it can perform multiple
%   function evaluations at once.
%      ezsurf('x.*exp(-x.^2 - y.^2)')
%
%   You may also use a function handle to an existing function. Function
%   handles are more powerful and efficient than string expressions.
%       ezsurf(@peaks)
%
%   EZSURF plots the variables in string expressions alphabetically.
%      subplot(1,2,1), ezsurf('u.*(v.^2)./(u.^2 + v.^4)')
%   To avoid this ambiguity, specify the order with an anonymous function:
%      subplot(1,2,2), ezsurf(@(v,u)u.*(v.^2)./(u.^2 + v.^4))
%
%   If your function has additional parameters, for example k in myfun:
%      %------------------------------%
%      function z = myfun(x,y,k1,k2,k3)
%      z = x.*(y.^k1)./(x.^k2 + y.^k3);
%      %------------------------------%
%   then you may use an anonymous function to specify that parameter:
%      ezsurf(@(x,y)myfun(x,y,2,2,4))
%
%   See also EZPLOT, EZPLOT3, EZPOLAR, EZCONTOUR, EZCONTOURF, EZMESH, 
%            EZSURFC, EZMESHC, SURF, VECTORIZE, FUNCTION_HANDLE.

%   Copyright 1984-2003 The MathWorks, Inc.
%   $Revision: 1.9.4.4 $  $Date: 2004/11/29 23:33:14 $

% Parse possible Axes input
[cax,args,nargs] = axescheck(varargin{:});

if ~isempty(cax)
    h = ezgraph3(cax,'surf',args{:});
else
    h = ezgraph3('surf',args{:});
end

if nargout > 0
    hh = h;
end