z_solve_LU_sp-mz 3.31; sp-mz.D.16