As /: Table[As[iter___], ℰ__] := Table[ℰ, {iter}];
Rasterize[GraphicsGrid[Table[As[i,1,5], As[j,1,5],
q={a (x^2)+(Sqrt[2])b x y+c (y^2),a (x^2)+(Sqrt[2])b x y+c (y^2)} /. a|b|c :>
RandomVariate[NormalDistribution[]];;
Show[Join[
{VectorPlot[(q/(((q.q))^(3/8))), {x,-1,1}, {y,-1,1}, Frame->None]},
Table[As[k,1,100],
eqns=Join[
Thread[{x'[t], y'[t]} == (q /. {x->x[t], y->y[t]})],
{x[0]==RandomReal[{-1,1}], y[0]==RandomReal[{-1,1}]}];
sol=NDSolve[eqns, {x,y}, {t,-1,1}];
ParametricPlot[
Evaluate[{x[t], y[t]} /. sol], {t,-1,1},
PlotRange -> {{-1,1},{-1,1}}, ColorFunction -> (Red&)]]
]]
]], ImageSize -> 960]