function f=calcf(x);
beta=x(1); alpha=x(2);
z=[ones(111,1),hsdata(2:112,1)./hsdata(1:111,1),(hsdata(2:112,2)+ ..
   hsdata(2:112,3))./hsdata(1:111,2)];
h=beta*z(:,3).*z(:,2).^alpha-1;
z=z(1:110,:);
f=h(2:111)'*z*inv(z'*z)*z'*h(2:111);
[x(1),x(2),f]
