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