https://github.com/numpy/numpy/issues/15298 [reference](https://math.stackexchange.com/a/1916083/1896) I would like to work on this issue