I have an ODE equation which is sort of 

    y''[x] + 2 y'[x]/x + .0001 (y[x])^3 ==0 

subject to the boundary conditions 

    y'[0]==0  and y[Infinity]==0

Can anyone please suggest what would be a reliable process for solving this numerically in Mathematica?

`y[Infinity]` can obviously be truncated down to , say, `y[20]`.