function[L]=likeli(B,D,Omega,T,k1,x);
n=size(Omega,2);
L=-(-(T*n/2)*log(2*pi)+(T/2)*log(det(B(x(1:k1))).^2)-(T/2)*log(det(D(x(k1+1:size(x,2))))) -(T/2)*trace(B(x(1:k1))'*inv(D(x(k1+1:size(x,2))))*B(x(1:k1))*Omega));