Пояснение к скрипту

syms x ;%определение символьной переменной x%
y=-1.1*x^5+4.5*x^4-1.6*x^3-2.5*x^2+0.5*x;%определение символьной переменной y через x%
k=sym2poly(y);%преобразовать в вектор-строку полиномиальных коэффициентов%
x=roots(k);%нахождение корней уравнения %
plot(x,0,'ok');%строит на одном графике функции%
hold all;%обьединение в одной окне Figure(1) бвух трехмерных рисунков%
n=length(x);%Нахождение количества корней уравнения%
x=sort([-3;x;5]);%отсортировать значения вектора-столбца по возрастанию %
x((x<-3)|(x>5))=[];%удаление из вектора x значений меньше -3 или больше 5%
lx=x;
lx(length(x))=[];
rx=x;
rx(1)=[];
mx=(lx+rx)/2;%Вектор mx содержит экстремумы функции%
my=subs(y,'x',mx);%Вычисления значений функции в экстремумах%
by=(my>0);%нахождение максимума и минимума%
for i=1:length(lx);
    xx=lx(i):0.01:rx(i);
    yy=subs(y,'x',xx);
    if(by(i))
      plot(xx,yy,'-r');
    else
      plot(xx,yy,'-b'); 
    end;
end;    
title({'Анализ функции f=-1.1*x^5+4.5*x^4-1.6*x^3-2.5*x^2+0.5*x'});
dy=diff(y, 'x');%Нахождение производной функции в символическом виде%
kd=sym2poly(dy);%преобразовать в вектор-строку полиномиальных коэффициентов%
xe=roots(kd);%нахождение корней уравнения %
ye=subs(y, 'x', xe);%отсортировать значения вектора-столбца по возрастанию %
te=subs(diff(y,2),'x',xe)<0;%нахождение максимума и минимума%
for i=1:length(xe);
    if  te(i)
        str={'max';num2str(xe(i),'%2.2f');num2str(ye(i),'%2.2f');'\downarrow'};
    else
        str={'min';num2str(xe(i),'%2.2f');num2str(ye(i),'%2.2f');'\downarrow'};
    end;
    text(xe(i),ye(i),str,'HorizontalAlignment','center','VerticalAlignment','bottom');
end;
p = sym2poly(y);
g = @(x) abs(polyval(p,x));
zzz=quadl(g, -3, 5);
s=num2str(zzz);
xlabel(strcat('Площадь под кривой равна:',s));

lx=x; % что тут комментировать? просто присвоение
lx(length(x))=[]; % удаляется последний элемент
rx=x; % тоже
rx(1)=[]; % удаляется первый элемент