added Ming Gong's fortran code printout and the published paper from Japanese group.