Numerical approximation of axisymmetric positive solutions of semilinear elliptic equations in axisymmetric domains of R 3