%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%% in_loop.m %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function y = in_loop (x)
% a が最小にすべき関数
global C D Lambda Sigma
y=[0.5*(C-D*x)'*(C-D*x)+0.5*Lambda*sum(log(1.+x.*x/(Sigma*Sigma)))];
%%%%%%%%%%%%%%%%%%