function f=minRSS (x,y)
minRSS=999999;% 初始RSS值為一極大數
bestlambda=0;% 初始lambda
for lambda=-2:0.005:2;% 比較[-2，2]之間的lambda值
transdat=boxcox (lambda,y);
J=1;
n=length (y);
for i=1:n  % 這道題有53個數據，這個循環是為了將53個數據連乘
    J=J*(y (i)^(lambda-1));% 變換的Jacobi行列式
end
    z=transdat./J^(1/n);
    E=eye (n);  % 單位陣
    RSS=z'*(E-x*inv (x'*x)*x')*z;
    if (RSS<minRSS)
        minRSS=RSS;
        bestlambda=lambda;
    end
end
fprintf ('the best lambda=% g',bestlambda)