Use scipy.optimize.minimize_scalar

Replace brute force method with one that uses scipy.optimize.minimize_scalar to find L1 norm by minimizing sum of Euclidean distances