Gravity Field Vis

n = 6;
ms = RandomReal[{0, 0.1}, n];
ps = RandomReal[{0, 1}, {n, 2}];
resolution = 50;
V[x_, y_] :=
  Sum[(-ms[[i]] ({x, y} - ps[[i]]))/(
   ms[[i]]/100 + Norm[({x, y} - ps[[i]])]^2), {i, n}];
LineIntegralConvolutionPlot[{V[x, y], {"noise", resolution,
   resolution}}, {x, 0, 1}, {y, 0, 1}, ColorFunction -> "ArmyColors",
 LightingAngle -> 0, LineIntegralConvolutionScale -> 3,
 Frame -> False, ImageSize -> resolution]