a) In Python, define a set of 5 4x4 matrices containing random floating-point numbers between 0 and 1. Use the appropriate routines in numpy.linalg and scipy.linalg to find the determinants of each matrix, in the case of numpy without using a loop. Make a quantitative comparison of the results between the two methods.

b)

Use the function from part a and the linalg functions to solve for the equilibrium positions of the masses (where F_i=0) for the cases

(i) k = [1 2 3 4], L = [1 1 1 1], L_w=10

(ii) k = [0 1 1 0], L = [2 2 1 1], L_w=4