Сетка-сетка, которую вы создали, основана на тэте и радиусе.Однако радиус является постоянным для внешней части трубы, поэтому вместо этого он должен основываться на тэте и z, поскольку это две независимые переменные, определяющие сетку.Исходя из этих рассуждений, я считаю, что вам нужно следующее:
r = 0:0.1:3;
z = 0:0.1:10;
% set cut planes angles
theta1 = 0;
theta2 = pi*135/180;
nt = 101; % angle resolution
figure(1);
clf;
% create a grid over theta and z
t3 = linspace(theta1, (theta2 - 2*pi), nt);
[tt3, zz3] = meshgrid(t3, z);
% convert from cylindical to Cartesian coordinates
xx5 = r(end) * cos(tt3);
yy5 = r(end) * sin(tt3);
% plot surface
h5 = surface(xx5, yy5, zz3, 'EdgeColor', 'none');
% extra stuff to make plot prettier
axis vis3d
axis equal
view(3)
camzoom(0.7);