I use dsolve on a given system of ODE

sol=dsolve({sys_ode,ics},type=numeric,stiff=true);

and when I ask Maple to compute sol(1), I would like also Maple to display all the time points and the values of the solutions at those points that was computed to get to the endpoint=1 (starting from initial values in ics).

Any suggestions!?

Many thanks!