function cyl_surf()

clf

N = 20;
u = linspace(-0.1,1,N);
v = linspace(-0.1,1,N);

[um,vm] = meshgrid(u,v);

% Plot parabolic cylinder x = y^2
xm = um.^2;
ym = um;
zm = vm;

p = patch(surf2patch(xm,ym,zm));
set(p,'edgecolor','none');
set(p,'facecolor','r');
set(p,'cdata',1);
grid on

% Plot the plane y = pi/4
xm = um;
ym = 0*um + pi/4;
zm = vm;
p = patch(surf2patch(xm,ym,zm));
set(p,'edgecolor','none');
set(p,'facecolor','g','facealpha',0.5);
set(p,'cdata',1);
grid on

% Plot the plane y = 0
xm = um;
ym = 0*um;
zm = vm;
p = patch(surf2patch(xm,ym,zm));
set(p,'edgecolor','none');
set(p,'facecolor','g','facealpha',0.5);
set(p,'cdata',1);
grid on


% Plot the plane z = y
hold on;
xm = um;
ym = vm;
zm = ym;
p = patch(surf2patch(xm,ym,zm));
set(p,'edgecolor','none');
set(p,'facecolor','b');
set(p,'facealpha',0.5)
set(p,'cdata',1);

% Plot the plane z = 0
hold on;
xm = um;
ym = vm;
zm = 0*um;
p = patch(surf2patch(xm,ym,zm));
set(p,'edgecolor','none');
set(p,'facecolor','b');
set(p,'facealpha',0.5)
set(p,'cdata',1);

% Add labels, etc
xlabel('x','fontsize',16);
ylabel('y','fontsize',16);
zlabel('z','fontsize',16);
title('Problem 7','fontsize',18);

set(gca,'zlim',[-0.1,1]);



daspect([1,1,1]);
view(3)


shg






end