The program of Dr. Wen-Hsiung Li provides such an estimation. His address: li at gsbs18.gs.uth.tmc.edu His program is in FORTRAN. I have unix-C version of it. Andrey
