计算数学通讯 2004-1 第3页
