Question: comparing solutions of an ode

I have an ode y''(x) + ...=0 that depends on a parameter a and I would like to see how the solutions vary with a. I created a procedure using odeplot to yield a sequence of plots and display them. Not bad but I really need to look more closely. Plotting y for a=7 and overlaying a plot for y=7.1 is too crude. I want to look at y(x,7.1) - y (x,7) or y(x,7.1) /y (x,7) for x=0..3 say. The obvious solution is to compute y(x,7.1) for x=0..3 and save it as a vector and also y(x, 7). How can I do that so that the pointsx at which y(x, 7.1) are the same as the pointsx for which y(x,7) are evaluated. In other words, how can I specify to dsolve the intermediate points xj for which y(x) is calculated? 

Please Wait...