%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%                            Michael Pokojovy                             % 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[X, Y] = meshgrid(-2:0.1:2);

figure(1);

set(gcf, 'PaperUnits', 'centimeters');
xSize = 12; ySize = 14;
xLeft = (21 - xSize)/2; yTop = (30 - ySize)/2;
set(gcf,'PaperPosition', [xLeft yTop xSize ySize]);
set(gcf,'Position',[0 0 xSize*50 ySize*50]);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
hold on;
t = title('Der Kegel $\big\{\big(x, h\big(1 - \frac{|x|}{r}\big)\big)^{T} \,\big|\, 0 < |x| < r\big\}$, $r = 3$, $h = 5$', 'interpreter', 'latex');
set(t, 'FontSize', 16);

R = 3;
H = 5;

M = 30;
N = 20;

phi = linspace(0, 2*pi, M);
r   = linspace(0, R, N);

X = zeros(M, N);
Y = zeros(M, N);
Z = zeros(M, N);

for i = 1:M
for j = 1:N
    X(i, j) = cos(phi(i))*r(j);
    Y(i, j) = sin(phi(i))*r(j);
    Z(i, j) = H*(1 - r(j)/R);
end
end

grid on;

colormap gray;
surf(X, Y,  Z);

view([138 26]);

xlabel('x_{1}');
ylabel('x_{2}');
zlabel('f(x_{1}, x_{2})');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%