Please, I would like to know how to use a matrix function in maximum likelihood estimation of the below specifications ( asterik variables are function of parameters) 

series e1 = y - y*

series e2 = x - x*

series v1 = var(e1)

series v2 = var(e2)

series v12 = cov(e1,e2)

matrix U = {e1;e2}

matrix V = {v1,v12;v12,v2}

mle ll = -N/2*log(2π) - 1/2*log|V| - 1/2*U'*inv(V)*U

Is the above specifications possible? Please kindly help. 

Regards. 

Timmy