add: solve linear equations with lu decomposition