取自《C常用算法程序集-徐士良》
/************ 常用矩阵函数****************
//动态分配连续存储的二维数组
double ** ArrayAlloc(int m,int n)
//释放连续存储的二维数组
void ArrayFree(double **y)
//屏幕输出
void Arrayprint(double **array,int m,int n)
//形成单位矩阵
void eye(double *q,int n)
//转置矩阵
void trans(double *A,double *B,int m,int n)
//矩阵相加
void plusAB(double *A,double *B,double *C,int am,int an,int bm,int bn)
//矩阵相减
void subAB(double *A,double *B,double *C,int am,int an,int bm,int bn)
//矩阵相乘
void mulAB(double *A,double *B,double *C,int am,int an,int bm,int bn)
//矩阵数乘
void mulkA(double *A,double *B,int m,int n,int k)
//矩阵求逆(全选主元&高斯约当消去法)
double *inv(double *A,double *Ainv,int n)
//矩阵求行列式值
double det(double a[],int n)
//矩阵求秩
int rank(double a[],int m,int n)
//矩阵作初等变换变成行阶梯矩阵
void chhbg(double A[],double a[],int n)
//LU三角分解
int LU(double a[],int n,double l[],double u[])
//求解 AX=b 线性方程组,A是n*n矩阵,b是n*1列向量
int agaus(double A[],double X[],double b[],int n)
//求解 AX=B 线性方程组 A是n*n矩阵,B是n*m矩阵
int agjdn(double A[],double B[],int n,int m)
********************************************************************/
/************************* 函数实现 ***************************************************/
#ifndef MATRIXLIB_H
#define MATRIXLIB_H
#include <math.h>
#include <stdlib.h>
#include <stdio.h>
