Solving matrix equations in Julia