clear;
x=[1,2,3,4,5,6,7,8,9,10];
y=[1,4,9,16,25,36,49,64,81,100];
scatter(x,y);
len_x = length(x);
max_x = max(x);
min_x = min(x);
num = 100;
delta = (max_x-min_x)/num;
R = (max_x-min_x)*3/10;
a =[];
b =[];
for i=0:num
x1 = min_x + i*delta;
A = zeros(2,2);
B = [];
