Evaluate a differential equation at a given point I have a second order differential equation $x^2\cdot y''(x)+x\cdot y'(x)+(x^2-10)=0$ with the initial condition $y(100)=1$ and $y'(100)=0$. I want to evaluate $y(x)$ when $x=103$ with 8 digit accuracy.
The first step is rewrite this equation to a system of first order equation.
$\begin{cases} y'=w\\ w' = -\frac{w}{x}-\frac{10}{x^2}+1\end{cases}$ with the initial condition $\begin{cases} y(100)=1\\ w(100)=0\end{cases}$. Then I can use fourth-order Runge-Kutta methods to solve the system. In order to achieve 8 digit accuracy, I need to decide the step size $h$. Solving $h^4\le 10^{-8}$, I get $h=0.01$.
Is this seems a reasonable approach? Is there anything special with such a differential equation with the initial condition so that I need use other method to evaluate it?
 A: It is reasonable to begin with $h=0.01$. Since the theoretical bounds for RK4 are very conservative and also difficult to use, the step halving method may be preferable. It is based on running the computations with two or several step sizes and  comparing the results. 
For illustration purposes, I used Euler's method (rather than RK4) for this equation. The results were: 


*

*with $h=0.02$, $y(103)\approx 5.4225$

*with $h=0.01$, $y(103)\approx 5.4371$

*with $h=0.005$, $y(103)\approx 5.4443$


Since Euler's method is of 1st order, each result is off by about $Ch$ from the true value. We can estimate $C$ as follows
$$
\frac{5.4225-5.4371}{0.02-0.01} \approx -1.46
$$
$$
\frac{5.4371-5.4443}{0.01-0.005} \approx -1.44
$$
The fact that two values of $C$ are close allows us to conclude confidently that the error is bounded by, say, $2h$ when $h$ is small. 
For RK4 the error behaves as $Ch^4$, so you would have the differences like $0.02^4-0.01^4$ in the denominator. Try it (and report the results, if you don't mind).  
A: I think that the method of Frobenius is applicable to this problem.
