function lsr1=lsfct1(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];
lsr1 = 0;
for i=1:15
   lsr1 = lsr1 + (temp(i)-sol1(xx(1),xx(2),xx(3),x(i)))^2;
end