We give numerical algorithms
for second order differential equations. More specifically we consider the problem of
numerically determining oscillation points and vectors for numerical solutions of the
equation (r(t)x′(t))′ + p(t)x(t) = 0 and focal points and vectors for the quadratic
form J(x) =∫ab(rx′2− px2)dt.