Y=X * beta
(y, n * 1 vector, X for n * p matrix, beta vector of p * 1)

Beta estimate=argmin | y - X * beta | solving

Note: see quantreg function in python, but I this is containing matrix, so don't know how to write parameters to coefficient of beta,