Reproduces the results of CPL 269, 161 (1997), and are written in FORTRAN 90