How to get coordinates from an ellipsoid function

I've got the function of an ellipsoid. The function looks like: Ax^2 + By^2 + Cz^2 + 2Dxy + 2Exz + 2Fyz + 2Gx + 2Hy + 2Iz = 1. I want to get coordinates of a large number of points(say 10000) located on the surface of this ellipsoid. These points don't have to be random, as long as they are evenly distributed (for example, if z ranges from 0 to 100, for every z=0:1:100, I can get 100 points). The output format should be 10000 rows, each with 3 columns for x, y, z. And also the numbers should be rounded to integers.

I'd appreciate your kind help. Thanks in advance.