The problem can be resolved by using the information T.
Let's say we have [T,Y] = ode45(myode,tspan,y0) and tspan has only two elements, e.g. tspan = [a,b]. Then, to get a 2-dimensional solution for ODE23/45, we write
Y = Y(abs(T-a)<0.0001 | abs(T-b)<0.0001)
which would result in a 2-dimensional solution as desired.