% This is about the sense and sensibility.
% x^2 + (y - x^(2/3))^2 = 1

% By Limin Shao on Tuesday of Jan. 17, 2012.

x = linspace(-1, 1, 400);

temp = x.^2;
y1 = sqrt(1 - temp) + exp(log2(temp)/3);
y2 = -sqrt(1 - temp) + exp(log2(temp)/3);

figure; hold on;
fill([x*1.03 x(end:-1:1)*1.03], [y1 y2(end:-1:1)], 'r', ...
    'EdgeColor', [1 0 0], 'LineWidth', 2);
fill([x x(end:-1:1)], [y1 - abs(x)*0.05 y2(end:-1:1)], 'w', ...
    'EdgeColor', 'r', 'LineWidth', 0.5)

set(gca, 'XLim', [-2 2], 'YLim', [-2 2], ...
    'Box', 'on', 'XTick', [-2 2], 'YTick', [-2 2], ...
    'XTickLabel', [], 'YTickLabel', []);

%plot(x, y1, 'r', x, y2, 'r', 'LineWidth', 2);