So i got rhis
with(plots); low := 1015; high := 1019; j := 0; for i from low to high do proc (X) options operator, arrow; X(t) end proc; proc (Y) options operator, arrow; Y(t) end proc; ode := diff(X(t), t) = evalf(subs(x = X(t), y = Y(t), subs(vvx = VX[i], vvx))), diff(Y(t), t) = evalf(subs(x = X(t), y = Y(t), subs(vvy = VY[i], vvy))); ics := X(0) = Cordx[i+1], Y(0) = Cordy[i+1]; j := j+1; plott[j] := odeplot(ds, [X(t), Y(t)], t = 0 .. 100, numpoints = 1000) end do; display(plott, jj)