Duh pre second cup of coffee/. Should be / df1b2vector tmp=solve(varcov,diff,ld,sgn); nLL+=.5*(nobs*log(2*M_PI)+ld+ tmp*diff); -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.admb-project.org/pipermail/users/attachments/20140502/d6706235/attachment.html>