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)