%This function models the first Eigen vector of Covariance matrices.
%Input:     b = 4x1 vector with the four parameters
%     :     x = point at which to calculate FitU1_new2
%Output:    U1(x)
function U1 = FitU1_new2(b,x)

alpha = b(1);
beta = b(2);
gamma = b(3);
delta = b(4);

U1 = alpha./((beta./x) + gamma ).^(delta);
