Need help with cylinders

Hi everyone. I want to plot and create the intersection of 2 perpendicular cylinders.

It should look like this.

Any help would be greatly appreciated!

[Walter Roberson]

That does not appear to be a simple plot of perpendicular cylinders to me. Perpendicular cylinders would have a hard boundary not a smoothly curving boundary.

[John D'Errico]

Walter - I think you are mistaken. That is what the intersection would look like, of two cylinders of the same radius, each aligned with one of the axes.

Visualize one cylinder along the x axis, the other along the y axis. Your point of view is above the cylinders, and off to the side.

As you can see from the cylinders, where they intersect, the surface would look as Mario wants it.

Answer 1

[Teja Muppirala]

% Make a figure (I like my figures in black)

colordef(figure,'black');

% Draw a parameterized surface

[THETA,Z] = meshgrid(linspace(-pi,pi,9),-1:0.01:1);
L = sqrt(1-Z.^2);
RHO = min( abs(L./cos(THETA)) , abs(L./sin(THETA)) ) ;

[X,Y,Z] = pol2cart(THETA,RHO,Z);
h = surf(X,Y,Z,'EdgeAlpha',0.1,'FaceColor',[0.9 0.9 0.9]);

axis equal
axis vis3d

% Add some lights

rng(0);
for n = 1:4
    light('Position',randn(1,3),'Color',rand(1,3));
end