МНК. Где ошибка?

clear; 
n=50; %Количество переменных x 
mistake = 1; %ошибка 
temp_x=1+(3-1)*rand(1,n); %Генерирование чисел 
x=sort(temp_x); %упорядочивание элементов массива по возрастанию 
LastX=x.^(-0.7); 
func=3*(x.^(-0.7)); %Функция 
while mistake < 0 || mistake > 0.01 %Ошибка измерения сводиться к 0 
little=-1+(1+1)*rand(1,n); 
mistake=sum(little)/n; 
end 
 
y=func+mistake; %Значения функции с учетом ошибки измерения 
disp(mistake); 
Mx=mean(LastX); 
My=mean(y); 
Mxx=mean(LastX.^2); 
Mxy=mean(LastX.*y); 
a=(Mxy-My*Mx)/(Mxx-Mx*Mx); 
b=(Mxx*My-Mx*Mxy)/(Mxx-Mx*Mx); 
NewFunc=a.*LastX+b; % линейный вид 
fprintf('Ошибка: %f\n', mistake); 
fprintf('Функция: %.4f*x^(-0.7)+%.4f\n', a, b); 
plot(x, NewFunc); 
hold on; % зафиксировали окно, т.е. пока hold не отпустим, все графики рисуются в этом окне. 
plot(temp_x, y, '+'); 
hold on; 
plot(x, y, '*r'); 
grid on; 
legend('Исходная функция','Исходная функция + шум','Метод наименьших квадратов')

plot(temp_x, y, '+');

plot(x, y, '+', x, NewFunc);

clear;
n=50; %Количество переменных x
mistake = 1; %ошибка
temp_x=1+2*rand(1,n); %Генерирование чисел
x=sort(temp_x); %упорядочивание элементов массива по возрастанию
LastX=x.^(-0.7);
func=3*(x.^(-0.7)); %Функция
while mistake < 0 || mistake > 0.01 %Ошибка измерения сводиться к 0
little=-1+2*rand(1,n);
mistake=mean(little);
end
y=func+little; %Значения функции с учетом ошибки измерения
disp(mistake);
Mx=mean(LastX);
My=mean(y);
Mxx=mean(LastX.^2);
Mxy=mean(LastX.*y);
a=(Mxy-My*Mx)/(Mxx-Mx*Mx);
b=(Mxx*My-Mx*Mxy)/(Mxx-Mx*Mx);
NewFunc=a.*LastX+b; % линейный вид
fprintf('Ошибка: %f\n', mistake);
fprintf('Функция: %.4f*x^(-0.7)+%.4f\n', a, b);
plot(x, func,x,y,'+',x,NewFunc,'*r'),grid on,...
legend('Исходная функция','Исходная функция + шум','Метод наименьших квадратов')