image=Compile[{{n,_Integer}},tmax=25;dt=0.1;h=0.25;g=0.2;mu=0.07;zlist={Sqrt[3]+I,–Sqrt[3]+I,–2I};Table[z=x+Iy;v=a=a1=0;Do[z+=vdt+(4a–a1)dt^2/6;vpredict=v+(3a–a1)dt/2;a2=Plus@@((zlist–z)/(h^2+Abs[zlist–z]^2)^1.5)–gz–muvpredict;v+=(2a2+5a–a1)dt/6;a1=a;a=a2,{t,0,tmax,dt}];r=Abs[z–zlist],{y,–5.0,5.0,10.0/n},{x,–5.0,5.0,10.0/n}],CompilationTarget–>"C",RuntimeOptions–>"Speed"];
