Непонятная ошибка

mstn=round(((sqrt(((xd(pp)-xsh(scnx))^2)+((yd(pp)-ysh(scny))^2)+((zd(pp)-zsh(scnz))^2)))/v)/diskr);

>> a = [1 2 3];
 
>> disp(a(1:3));
 
     1     2     3
 
>> disp(a(1:2.5));
 
Warning: Integer operands are required for colon operator
when used as index 
 
     1     2

function maxpravdt
tic
kolt=30;                             %количество точек kolt x kolt x kolt
grx=3200;                            %граница по оси x от нуля
gry=3200;                            %граница по оси y от нуля
grz=-3200;                           %граница по оси z от нуля
rl=800;                              %размер луча
%координаты расположения датчиков
 
kluc=8;                              %количество лучей
kd=5;                                %количество датчиков на луче
kol=(kd*kluc)+1;                     %суммарное количество датчиков + центральный
clux=1600;                            %центр расположения "солнца" по оси x
cluy=1600;                            %центр расположения "солнца" по оси y
sch=0;
for kold=1:kd
    rad=rl-((rl/kd)*(kold-1));
   for rd=1:kluc
     sch=sch+1;  
   cd(1,sch)=round((rad*cos(((2*pi)/kluc)*(rd-1))))+clux;
   cd(2,sch)=round((rad*sin(((2*pi)/kluc)*(rd-1))))+cluy;
     %save('ydat.mat', 'yd');
   end
end
pm=sch+1;
cd(1,pm)=clux;
cd(2,pm)=cluy;
save('cdat.mat', 'cd');
 
% coord=[400 600 800 1000 1200 1400 800 800 800 800 ; 800 800 800 800 800 800 400 600 800 1200]; 
% disp(coord);
% save('coord.mat', 'coord');
 
% xd = coord(1,:);
% yd=coord(2,:);
xd = cd(1,:);
yd=cd(2,:);
zd=zeros(1,kol);
% plot3(X,Y,Z);
 
v=3000;                             %скорость распространения волны в среде
f=20;                               %частота синуса
koleb=1/f;                          %длительность одного колебания
 diskr=1/160;                        %время дискретизации
% sds=90;                           %сдвиг в массиве vsts
% sds=30;
 
xcr=1600;                             %координата точки удара по оси x
ycr=1600;                              %координата точки удара по оси y
zcr=-1600;                            %координата точки удара по оси z
 
xsh=[linspace(0,grx,kolt)];         %разбитие промежутка по оси x
ysh=[linspace(0,gry,kolt)];         %разбитие промежутка по оси y
zsh=[linspace(0,grz,kolt)];         %разбитие промежутка по оси z
 
% disp(xsh);
%расчет времени отклика для произвольно взятой точки и амплитудного
%коэфф-та
 
for p=1:kol
    tcr=round(((sqrt(((xd(p)-xcr)^2)+((yd(p)-ycr)^2)+((zd(p)-zcr)^2)))/v)/diskr);
    masstcr(p,:) = tcr;
    amp(p,:)=(round(100*(round(sqrt(((xd(p)-xcr)^2)+((yd(p)-ycr)^2)+((zd(p)-zcr)^2)))/1000))/100)-1;
    
end
save('masstcr.mat', 'masstcr');
save('amplify.mat', 'amp');
 
 
% дискретизация функции синус
 
% per=3;
% ckl=((koleb/diskr)*per)+1;            %количество итераций цикла
% for u=1:ckl
%     funct=round((sin(f*2*pi*diskr*(u-1)))*100)/100;
%     massin(u,:) = funct;
% end
 
% % %   disp(massin);
% save('sinus.mat', 'massin');
 
%дискретизация вейвлета
tau=2;
diskre=0.2;
ckl=(((6*tau)/diskre)+1);
 
%(6*(tau)/(diskr))+1;
% for u=1:ckl
%     funct=round((cos(f*2*pi*((u*1)-(3*tau))*diskr)*exp((-(((u*1)-(3*tau))^2))/(tau*tau)))*10000)/10000;
% massin(u,:) = funct;
% end
for u=1:ckl
% funct=round((cos(f*2*pi*((u)-(3*tau))*diskr)*exp((-((((u*1)-(3*tau)))^2))/(tau*tau)))*10000)/10000;
cosin=cos(f*2*pi*(((u*diskre)-diskre)-(3*tau))*0.01);
 expon=exp((-((((u*diskre)-diskre)-(3*tau)).^2))/(tau*tau));
 GRW(u,1)=cosin.*expon;
 end
 save ('GRW.mat','GRW');
 
 %блок добавления шума
%   rnd=(randn(ckl,1))/4;
%   save('rnd.mat','rnd');
%   GRWR=GRW+rnd;
%   save('GRWR.mat','GRWR');
%  
% massin(u,:) = funct;
%  end
%X=0:0.1:(tau*6);
% cosin=cos(f*2*pi*(X-(3*tau))*0.01);
% expon=exp((-((X-(3*tau)).^2))/(tau*tau));
% Y=cosin.*expon;
% 
% plot(X,Y,'r');
%   disp(massin);
% save('sinus.mat', 'massin');
 
%получение массива синус-сигнала плюс задержка от произвольно взятой точки
%со сдвигом sds
 
 
for l = 1:kol
    msvs = masstcr(l,:);
    %disp(msvs);
    SP(l,(msvs+1):(msvs+ckl)) = GRW(1:ckl)*amp(l);
end
save('theory.mat', 'SP');
 
%блок вариации амплитуды
kraz=40;
for vara=1:(kraz)
rando1=rand(1,1);
a1=-1.5; b1=1.5;
varam(1,vara)=(a1+(b1-a1)*rando1);
end
varamp=sort(varam);
save('varamp.mat', 'varamp');
 
 
 
%блок изменения времени
stept=0.01;
tend=5;
for chant=1:tend
timec(chant)=(stept*chant)-stept;
end
 save('timestep.mat','timec');
 
%    msatlabpool
%массивы
stch=0;
for time=1:tend
     swimt=(timec(time))/diskr;
    
  for scnx=1:kolt
    for scny=1:kolt
        for scnz=1:kolt
            stch=stch+1;
            MT1=[];
            
            %блок шума
            rnd=(randn(ckl,1))/4;
            GRWR=GRW+rnd;
            %конец блока шума
            
            for pp=1:kol
                mstn(1,pp)=round(((sqrt(((xd(pp)-xsh(scnx))^2)+((yd(pp)-ysh(scny))^2)+((zd(pp)-zsh(scnz))^2)))/v)/diskr);
            MT1(pp,(mstn(1,pp)+1+swimt):(mstn(1,pp)+ckl+swimt)) = GRWR(1:ckl)*amp(pp);
            end
%           mstn=[];
            for va=1:kraz
                
             
               
%                 masstn(pp,:)=mstn;
%                 msobm=masstn(pp,:);
                MT=MT1*varamp(va);                           
              
            
            
            
          if size(MT, 2) < size(SP, 2)
           big = SP;
           small = MT;
          else
           big = MT;
           small = SP;
          end
   
    
          [rows cols] = size(small);
          MV = big;
          MV(1:rows,1:cols) = big(1:rows,1:cols) - small;
          ME=1;
          
          for i1=1:kol
              for j1=1:length(MV)
                  ME=ME*exp(-(MV(i1,j1)^2));
              end
          end
                  MEA(1,va)=log(ME);
                  
            end
           
%             save(['MEA', num2str(stch), '.mat'], 'MEA');
%             disp(MEA);
%             for vmpa=1:(kraz-1)
%                 if MEA(1,vmpa+1)>MEA(1,vmpa)
%                 peramp=vmpa+1;
%                 else
                    peram=find(MEA==max(MEA));%vmpa;
%                 end
%             end
%             disp(peram);
                  MLN(scnx,scny,scnz)=((varamp(peram)).^2);
                  
%             save(['mst', num2str(stch), '.mat'], 'MT');
              
            
        end
    end
  end
% matlabpool close
   save(['mln', num2str(time), '.mat'], 'MLN');
  
end
savevtk(MLN, 'magic1_n__.vtk');
 
figure;
a=squeeze(MLN(15,:,:));
surf(a,'DisplayName','MLN');figure(gcf)
toc
end

cd = zeros(2, kd*kluc+1);