clf;
x1 = -2:0.01:2;
y1 = x1.^2;
y2 = ones(size(x1))*4;
x = [x1; x1];
y = [y1; y2];
patch(x, y, 'y')

