Как растянуть ось 3D графика

function [T,yma]=sloypr(x,y,z,q,v,del,vp,scor);
% Расчет температуры плоского слоя в предельном состоянии.
% Соответствует расчетам Кархина В.А.
%________________________________________________________
%% Пример задания исходных данных (в командном окне): 
%>>x=[1:-0.01:-10];  (см)
%>>y=[-3:0.01:3];    (см)
%>>z=0;              (см)        
%>>q=10000;          (мощность источника, Вт=Дж/с)
%>>v=1;              (скорость сварки, см/с)
%>>del=1;            (толщина слоя, см).
%>>vp=[1480 1100 1000 Ас3 Ас1 600]; (контролируемые температуры)
%>>scor (признак плоскости анализа: =1(x,y),=2(y,z),=3(x,z)) 
%>>[T,yma]=sloypr(x,y,z,q,v,del,vp,scor); (Запуск программы расчета).
 
%никоуглеродистая_______________________________________________________
lam=0.38;
a=0.079;
cro=lam/a;
yma=1;
 
nf=40; % Количество фиктивных источников
 
nx=length(x);
ny=length(y);
nz=length(z);
 
if (nx~=1)&(ny~=1)&( nz~=1)
    error('***** Координатное пространство должно быть двумерным !!!!')
else
qr=q/(4*pi*lam);
va=v/(2*a);
 
for m=1:1:nz;
    for n=1:1:ny;
        for k=1:1:nx;  
            
          Tc=0;             % Обнуление суммы
          
            for nf1=1:1:21; %Организация циклического суммирования 
                nf=nf1-10;  % фиктивных источников, ограниченных номером 21 
 
   R=sqrt(x(k)^2+y(n)^2+(z(m)-2*nf*del)^2); % Радиус-вектор источника nf
                
                if R<del/1000
                    R=del/1000;
                end% if
                
                dT=qr/R*exp(-va*(R+x(k))); % Формула 2.26
                Tc=Tc+dT;    % Суммирование температур действительного и
                             %  фиктивных источников 
            end  %nf1
            
            if Tc>1500
                Tc=1500;
            end %if
          
            if nz==1
             T(n,k,m)=Tc;     %  Для x,y   
             coor=y;
            elseif nx==1
             T(m,n,k)=Tc;     %  Для y,z
             coor=z;
            else
             T(m,k,n)=Tc;    %  Для x,z
             coor=z;
      
         end
               
      end %x
    end %y
end %z
%____________________________________________________
% Расчет ширины заданных изотерм на поверхности (z=0)
%____________________________________________________
 if nz==1     
       Te=T(:,:,1);      
            nvp=length(vp);
for iv=1:1:nvp
   
    Tr=Te-vp(iv);
    [miT,nn]=min(abs(Tr));
       
        ym=y(nn);
        [my(iv),nym(iv)]=max(abs(ym));
        
    end
    yma=my
end %if 
end %if пространство
 
if scor==1
    pcolor(x,y,T)
    
elseif scor==2    
pcolor(y,z,T)
%set(gca,'YDir','reverse')
else scor==3
 pcolor(x,z,T)
%set(gca,'YDir','reverse')
end %if
set(gca,'YDir','reverse')
axis equal tight
shading interp
grid

x=[1:-0.01:-27];
y=[-3:0.01:3];
z=0;    
q=15300;
v=1;
del=0.5;
vp=[1480 1100 1000 855 715 600];
scor=1;
[T,yma]=sloypr(x,y,z,q,v,del,vp,scor);
contour(x,y,T,vp);
grid;
surf(x,y,T);
shading interp