Numerical test for ordinary differential equation

¸ñ·Ï