I don't have any experience with either python or blender. I've been plotting implicit and parametric functions through a 3D graphing calculator, called calcplot3D.
I also started off with my imagination, too. After a long, continued interest, a lot of math experiments and playing around, along with input from others, I developed a few methods for getting equations for shapes, and using them to render slices and projections.
Don't get me wrong, you can still extrapolate a lot of info without the equation, and find accurate details about the shape. The equation ran through a computer is, however, much more revealing. It will show you things you didn't know, too!
If you want to use python/blender with an implicit equation, you'll likely have to implement a marching cubes algorithm.
-----------------
I just so happen to have rendered the coninder myself, recently. It was one of many experiments to see if my new parametric function algorithm is working. It seems to be doing a great job! Check out this projection of a rotating coninder:
The parametric function I discovered for a unit coninder is :
r(x,y,z,w) = { (v-1)u*cos(t)√3 , (v-1)u*sin(t)√3 , 3v+1 , 2s√3 } | u,v,s ∈ [-1,1] ; t ∈ [0,π]
Now, you may notice this is a 4-manifold embedded in E^4 , so it's a solid 4D shape. I found another algorithm that decomposes this equation into 0,1,2, and 3-manifolds embedded in E^4 . These are the faces and edges that you see in the animation (only the 1D and 2D surfaces).