Инициализация массива

mass1(i)=z;

 mass2=x(m);

x=[];

if(i>m) 
z=array_z(i-m);

m=1000

m=1000;
 
if(i>m) 
z=array_z(i-m);

if(i>1000) 
z=array_z(i-m);

clear all;
clc;
 
%Исходные данные
 
m = 1000; %число точек
n = 10; %число оборотов
syms x z xnew znew i;  %объявление символьных переменных
p = 1000;  %частота (Гц)
w = 0.1*p;
R = 20;    %радиус образца(мм)
Vokr = w/R; %линейная скорость в окружном направлении
t1 = 2*pi*R/Vokr; %время полного поворота (на 2Pi)
f = 0.1; %подача (мм/об)
Vos = f/t1; %скорость в осевом направлении
Krez = 1000; %коэффициент жесткости  (МПа)
K = 500; % задаемся
b = 0.1*R; %ширина резания
delt = t1/m; %шаг времени на первом обороте (на всех оборотах)
m = K/p^2; %масса резца
 
%Создаем массивы
 
t=0:delt:t1*n; %массив времени
array_z = [m*n+1]; %Почему подчеркивает скобку(как задавать массивы)?
array_x = [m*n+1];
array_h = [m*n+1]; %не нужен
 
%Цикл вычисления смещений
 
x = 0; z = 0;  % в начальный момент времени
 
for i=1:1:(m*n+1)
    if (-x>z) %Если инструмент вышел
        h=0;
        z=z;
    elseif (-x<z)
        h=z+x;
        array_h(i)=h; %-
        z=-x;
    else h=0;
    end;
    if(i>m) %Если оборот не первый,
        z=array_z(i-m); %тогда ищем значение z в этой точке на предыдущем обороте
    else z=z;
    end;
    znew = z+Vos*delt;
    z=znew;
    array_z(i)=z;
    Frez = -Krez*b*h;
    if (i>m)
        x=array_x(i-m);
    else
        x=x;
    end;
    V_mgn = x/delt;
    xnew = x*cos(p*delt)+V_mgn/p*sin(p*delt)+Frez/(m*p)*2*(sin(p*delt/2))^2/p;
    x=xnew;
    array_x(i)=x;
end;
figure
plot(t,array_x);