function [hyperfit]=hypfit(xx,yy)

        x_hy=(1/xx)';

        y_hy=(1/yy)';

        x_sum=sum(x_hy);

        y_sum=sum(y_hy);

        slope=y_sum/x_sum;

        intercept=value??; %%?

        m=1/intercept;

        k=slope*m;

        y_hy_fit=k.*x_hy+m;

        hyperfit=(y_hy_fit);

    end