//奇异值分解法求广义逆

//本函数返回值小于0表示在奇异值分解过程,

//中迭代值超过了60次还未满足精度要求.

//返回值大于0表示正常返回。

//a-长度为m*n的数组，返回时其对角线依次给出奇异值，其余元素为0

//m-矩阵的行数

//n-矩阵的列数

//aa-长度为n*m的数组，返回式存放A的广义逆

//eps-精度要求

//u-长度为m*m的数组，返回时存放奇异值分解的左奇异量U

//v-长度为n*n的数组，返回时存放奇异值分解的左奇异量V

//ka-整型变量，其值为max(n,m)+1

//调用函数：dluav()-//Singular value decomposition method for generalized inverse// this function return value is less than 0 in the singular value decomposition process,// in the iterative value of more than 60 are unable to meet the precision requirement.// Return value greater than 0 indicated a normal return.// a-length m* n array, return its diagonal singular value is given followed by the remaining elements to 0// m-matrix rows// n-matrix the number of rows// aa-length n* m array, retrievable storage A generalized inverse// eps-precision// u-length m* m of the array, return deposited the left singular value decomposition singular quantity U// v-length of n* n of the array, return deposited the left singular value decomposition singular volume V// ka-integer variables, its value is max (n, m)+ 1// call the function: dluav ()