---

clc
close all

% Добавляем знак ' в конце для транспонирования в столбцы
x = [1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7]';
y = [-2.97 -6.18 -11.25 -12.63 -13.51 -14.18 -18.37 -20 -24.44 -27.85 -28.95 -32 -35.5]'; 

plot(x, y, 'o')
xlabel('X')
ylabel('Y')
hold on

% Теперь x.^0 и x — это столбцы, и [x.^0 x] создаст матрицу N x 2
XR = [x.^0, x]; 

[b, bint, r, rint] = regress(y, XR);
Ymodel = b(1) + b(2)*x;
plot(x, Ymodel)

% Вывод уравнения на график
text(3, -5, ['y = ', num2str(b(2)), '*x + ', num2str(b(1))])

figure
rcoplot(r, rint)

ys = mean(y);
SSob = sum((y - ys).^2);
SSos = sum((y - Ymodel).^2);
SSF = sum((Ymodel - ys).^2);
R = 1 - SSos / SSob;
R1 = SSF / SSob;