function lsr=lsfct(xx)
x = [10 14 18 22 26 30 34 38 42 46 50 54 58 62 66];
temp = [ 115.95 101.94 90.18 80.01 71.37 63.75 56.94 51.97 48.21 44.48 40.66 38.51 37.81 36.11 35.18];
lsr = 0;
for i=1:15
   lsr = lsr + (temp(i)-sol(xx(1),xx(2),x(i)))^2;
end