Question: singularity value from dsolve warning

so i got this
> low := 170; high := 195; j := 0; byy := 5; plott := vector((high-low)/byy+1);
with(plots);
for i from low by byy 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 = VXX[i-1], vvx))), diff(Y(t), t) = evalf(subs(x = X(t), y = Y(t), subs(vvy = VYY[i-1], vvy)));
ics := X(0) = Cordx[i-1], Y(0) = Cordy[i-1];
ds := dsolve([ode, ics], type = numeric);
j := j+1;
plott[j] := odeplot(ds, [X(t), Y(t)], t = 0 .. 100, numpoints = 10000, 'legend' = i, color = ColorTools:-Color([rand()/10^12, rand()/10^12, rand()/10^12]), 'thickness' = 1) end do;
plots[display](convert(plott, 'list'));

%;Warning, cannot evaluate the solution further right of 9.3249218, probably a singularity
Warning, cannot evaluate the solution further right of 10.312887, probably a singularity
Warning, cannot evaluate the solution further right of 11.668729, probably a singularity
Warning, cannot evaluate the solution further right of 13.881632, probably a singularity
Warning, cannot evaluate the solution further right of 21.248509, probably a singularity 

 

i need somehow get singularity value from warnings, or mb from dsolve directly. i just need to get X(signularity-eps) and Y(signularity-eps)

Please Wait...