Aside from the differences of BLAS and LAPACK in Julia and Matlab, it maybe meaningful to compare the algorithms between Julia and Matlab. For instance, sqrtm function requires to do a Schur decomposition for a normal square matrix which highly depends on LAPACK. How about benchmark sqrtm for an upper triangular matrix?