Построить усечённый конус

r1 = str2double(get(handles.edit5,'string')); %Нижнее основание
r2 = str2double(get(handles.edit6,'string')); %Верхнее основание
x0 = str2double(get(handles.edit7,'string')); % Координата точки
y0 = str2double(get(handles.edit8,'string')); % Координата точки
z0 = str2double(get(handles.edit14,'string')); % Координата точки
h = str2double(get(handles.edit12,'string')); % Высота
 
h1=r1/h;
h2=r2/h;
z2=h+z0;
x1 = -r1:h1:r1;
x2 = -r2:h2:r2;
x = [-r1:h1:r1;-r2:h2:r2]-x0;
y1=[sqrt(r1^2-x1.^2); sqrt(r2^2-x2.^2)]-y0;
y2=[-sqrt(r1^2-x1.^2); -sqrt(r2^2-x2.^2)]-y0;
z1(1:length(x1)) = z0;
z2(1:length(x1)) = z2;
z = [z1;z2];
 
surf(x,y1,z,x,'EdgeColor','none')
surf(x,y2,z,x,'EdgeColor','none')
set(gca,'PlotBoxAspectRatio',[1 1 1]);
meshc(x,y,z);

%Конус
clc
clear all
r1=5;
r2=2;
h=5;
z0=0;
y0=0;
x0=0;
n=20;
 
h1=r1/n;
h2=r2/n;
z2=h+z0;
x10 = -r1:h1:r1;
x11 = r1:-h1:-r1;
x20 = -r2:h2:r2;
x21 = r2:-h2:-r2;
xx1 = [x10 x11];
xx2 = [x20 x21];
x=[xx1;xx2];
y10=sqrt(r1^2-x10.^2); y11=-sqrt(r1^2-x11.^2);
y20=sqrt(r2^2-x20.^2); y21=-sqrt(r2^2-x21.^2);
y=[[y10 y11];[y20 y21]];
z1(1:length(xx1)) = 0;
z2(1:length(xx1)) = h;
z = [z1;z2];
x=x+x0; y=y+y0; z=z+z0;
 
surf(x,y,z,x,'EdgeColor','none')
%surf(x,y2,z,x,'EdgeColor','none')
set(gca,'PlotBoxAspectRatio',[1 1 1]);
%meshc(x,y,z);