Numerical method for dc Josephson current between d-wave superconductors