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

02