function [lambda,lik, x_box]=Box_Cox(a,b,data)
 
   n=length(data);
   lambda_series=a:0.01:b;
   
   Lik=zeros(1,length(lambda_series));
   for i=1:length(lambda_series)
       l=lambda_series(i);
       [Lik(i),x_box]=Box_Cox_lik(l, data);
   end
   [Lik_max, index]=max(Lik);
   lambda=lambda_series(index);
   [lik,x_box]=Box_Cox_lik(lambda, data);