x = [0; 1; 2; 3.5; 4];  y = [-0.5; 1; 2.4; 2.0; 3.1];
F = [ones(length(x),1) x];
[p,y_v,r,p_v] = LinearRegression(F,y)
sqrt(p_v)
