Plotting a section of a hemisphere using parametric equations

Using:

theta=[0:0.1:2*pi];
phi=[0:0.1:pi/4];
R=1;
[THETA,PHI]=meshgrid(theta,phi);
X1=R.*cos(THETA).*sin(PHI);
Y1=R.*sin(THETA).*sin(PHI);
Z1=R.*cos(PHI);
surf(X1,Y1,Z1)

I get a section of a sphere, but a slit is missing from the plot. How do i fill in the slit?

Answer 1

[David Goodmanson]

Hi Marcus

theta=[0:0.01:1]*2*pi;

ensures that the value of 2*pi is reached.

[Marcus Velling]

Works perfectly, thank you

[darova]

Answer 2

[Carlos Guerrero García]

Another way for changing the theta and/or phi values is

theta=[0:pi/36:2*pi];
phi=[0:pi/36:pi/4];

Anyway, I suggest not to create a vector for theta and another vector for phi only for making the meshgrid. I suggest changing that lines as follows:

[THETA,PHI]=meshgrid(0:pi/36:2*pi,0:pi/36:pi/4);
R=1;
X1=...