Sampling angles on a sphere: From a general covariance matrix to concentration parameter

I am trying to extract points on a sphere in python.

Sampling angles on a sphere: From a genral covariance matrix to concentration parameter

I am trying to extract points on a sphere in python.