Best numerical method to integrate a differential equation like $\frac{\partial y}{\partial t}=f\left(t,\frac{\partial y}{\partial r}\right)$? What is the best numerical method to solve an equation like this:
$$\frac{\partial y}{\partial t}=f\left(t,\frac{\partial y}{\partial r}\right)\quad?$$
Can somebody give, at least, a reference?

This is quite difficult to answer without further details. Assuming that you have no clue about what to do, I would recommend you to look into finite difference methods. These are probably the easiest methods to understand. Other methods are finite elements (bit more complicated and in 1d often similar to finite differences, though boundary conditions are handled differently) and (pseudo)spectral methods (which may give you better accuracy if everything is smooth).
A good reference is Morton & Mayers, Numerical solutions of partial difference methods (though as Ross says, any numerical analysis book will have something about this, often in the last chapter). There is also a list of links at http://math.fullerton.edu/mathews/n2003/finitediffpde/FiniteDifferencePDEBib/Links/FiniteDifferencePDEBib_lnk_1.html .
