function ezcontourf(f,varargin)
%EZCONTOURF Easy to use filled contour plotter.
%   EZCONTOURF(f) plots the contour lines of f(x,y) using CONTOURF where
%   f is a string or symbolic expression representing a mathematical
%   function of two variables, say 'x' and 'y'.  The function f is 
%   plotted over the default domain -2*pi < x < 2*pi, -2*pi < y < 2*pi.
%   The computational grid is chosen according to the amount
%   of variation that occurs.
%
%   EZCONTOURF(f,DOMAIN) plots f over the specified DOMAIN instead of the
%   default DOMAIN = [-2*pi,2*pi,-2*pi,2*pi].  The DOMAIN can be the
%   4-by-1 vector [xmin,xmax,ymin,ymax] or the 2-by-1 vector [a,b] (to
%   plot over a < x < b, a < y < b).
%
%   If f is a function of the variables u and v (rather than x and y),
%   then the domain endpoints umin, umax, vmin, and vmax are sorted
%   alphabetically.  Thus, EZCONTOURF(u^2 - v^3,[0,1],[3,6]) plots the
%   contour lines for u^2 - v^3 over 0 < u < 1, 3 < v < 6.
%
%   EZCONTOURF(...,N) plots f over the default domain using an N-by-N grid.
%   The default value for N is 60.
%
%   ROTATE3D always on.  To rotate the graph, click and drag with your mouse.
%
%   Examples:
%    f = ['3*(1-x)^2*exp(-(x^2) - (y+1)^2)' ... 
%       '- 10*(x/5 - x^3 - y^5)*exp(-x^2-y^2)' ... 
%       '- 1/3*exp(-(x+1)^2 - y^2)'];
%    ezcontourf(f,[-pi,pi])
%    ezcontourf('sin(sqrt(x^2+y^2))/sqrt(x^2+y^2)',[-6*pi,6*pi])
%    ezcontourf('x*exp(-x^2 - y^2)')
%    ezcontourf('-3*z/(1 + t^2 - z^2)',[-4,4],120)
%    ezcontourf('sin(u)*sin(v)',[-2*pi,2*pi])
%
%   See also EZPLOT, EZPLOT3, EZPOLAR, EZCONTOUR, EZSURF, EZMESH,
%            EZSURFC, EZMESHC, CONTOURF.

%   Copyright (c) 1984-98 by The MathWorks, Inc.
%   $Revision: 1.3 $  $Date: 1998/04/09 21:44:58 $
%   P. J. Costa and Clay Thompson

ezgraph3('contourf',f,varargin{:});
rotate3d off;
