Laplace's equation is nice to work with because it is homogeneous and linear and hence always separable. It's solutions are many, but they are always of a polynomial or exponential nature. Physically, Laplace's equation usually implies some superposition of wave-like behavior in a scalar field.
The equation for a solenoidal vector field (that the vector field has zero divergence) is also very easy to solve and is always separable, but it also doesn't reveal the entire picture. It cannot tell you how the x component of the vector field depends on an additive term containing the y and z variables, for example. In order to figure out the full behavior of the vector field, you must know its curl field as well.
EDIT: I should correct myself. The solutions to Laplace's equation are always of a polynomial or exponential nature for coordinates in a holonomic basis.






