Changed boundary condition at r=0 from zero-gradient to reflecting.