%% Fill Area Between Contours
% Fill the area between contours by setting the Fill input of fcontour to 'on'. If you want interpolated shading instead, use the fsurf function with its option 'EdgeColor' set to 'none' followed by the command view(0,90).
% Create a plot that looks like a sunset by filling the contours of
% {\rm erf}((y+2)^3) - e^{(-0.65((x-2)^2+(y-2)^2)}.

syms x y
f = erf((y+2)^3) - exp(-0.65*((x-2)^2+(y-2)^2));
fcontour(f,'Fill','on')

% Copyright 2012 The MathWorks, Inc.