clf
axis([-5 5 0 32])
box on
hold on

x = -5:1/64:5;
p = plot(x,[x.^2; 2.^x],'linewidth',2);

syms x real
z = solve(x^2 == 2^x);
plot(z,z.^2,'.','markersize',24,'markeredgecolor','k');
set(gca,'xtick',sort([-5; double(z); 5]))

x = 2:1/64:4;
y = [2.^x fliplr(x.^2)];
x = [x fliplr(x)];
silver = (192/256)*[1 1 1];
fill(x,y,silver)

legend(p,'string',{'x^2','2^x'},'location','north');
xlabel('x')
title('x^2 = 2^x')
hold off
shg