write the function my lu solve 2 that uses lu decomposition with permutation and triangular substitutions to solve the linear system of equations, this time using multiple right-hand sides. the function should return a 1d numpy array containing the value of displacement at the free end of the slinky for each given right-hand side.
