A sphere of radius R, centered at (x0, y0, z0) can be represented by (x - x0)^2 + (y - y0)^2 + (z - z0)^2 = R^2. Also, any point on the line pointing in the direction of (a, b, c) can be represented by x = ka, y = kb, z = kc.
Plug those in for x, y, z and solve for "k". When I solved for "k", I got k = (G +- √(G² - F*H) )/F where F = a² + b² + c², G = a*x0 + b*y0 + c*z0, H = x0² + y0² + z0² - R². Plug those two values of "k" in the equations of the line and you'll get the two points (x, y, z) where the line intersects the sphere.