Using mathematica, make a animated plot of the vector field

associated with a binary star.

The gravitational field of a single star can be taken as

P = (r-r0)/(abs(r-r0)^2+.1)^1.5

The position of the first star is r0= [0.5 cos(t),0.5 sin(t),0], the positionof the second star is r0 = [0.5 cos(t),0.5 sin(t),0]. Make the plot for x ranging from -1 to 1, y ranging from -1 to 1 and t ranging from 0 to 2. Use the VectorPlot[] and Animate[] commands.