function ezcontourfill(f,range,varargin)
% f is symbolic expression in x,y
%   ezcontourfill(f,[xmin xmax ymin ymax])            plot 19 contours
%   ezcontourfill(f,[xmin xmax ymin ymax],30)         plot 30 contours
%   ezcontourfill(f,[xmin xmax ymin ymax],[-1 0 1])   plot contours for f=-1,f=0,f=1
%
% Use "colorbar" to add colorbar to plot
%
% This is actually a surface viewed from above, so you can rotate this to see the surface.

ezsurf(f,range);                                      
surfcontour(varargin{:})   
material([1 0 0])
view(0,90); 
camlight(0,0,'infinite');                           
lighting phong
 


