clear(),clc()
xe1=[0;0;0;0;0;;]
xe2=[361;625;961;1444;1936;;]
xe3=[19;25;31;38;44;;]
Y=[19.0;32.3;49.0;73.3;97.8]
X=[ones(length(xe1),1) xe1 xe2]
B,resnorm,residual, exitflag, output=lsqlin(X, Y)
B0=B[1];B1=B[2]
er= Y.- mean(Y)
er=(norm(er,2))^2
R2=1-resnorm/er
println("B0=", B[1],'\n'
,"B1=", B[2],'\n'
,"B2=", B[3])
println("resnorm=",resnorm)
println("residual=",residual)
println("R2=",R2)