#ifndef MAXTRIX_H  //如果被定義了，則跳到下一個#else或下一個#endif
                   //避免文件中的內容被定義兩次
#define MAXTRIX_H

//兩個矩陣相加(重載)
void matrix_add(double *,double *,double *,int,int);
void matrix_add(double *,double *,double *,int);
    //前面兩個函數的計算結果保存在第三個矩陣
void matrix_add(double *,double *,int,int);
void matrix_add(double *,double *,int);
    //這兩個函數的計算結果保存在第一個矩陣中


//兩個矩陣相減(重載4個函數)
void matrix_minus(double *,double *,double *,int,int);
void matrix_minus(double *,double *,double *,int);
   //前面兩個函數的計算結果保存在第三個矩陣
void matrix_minus(double *,double *,int,int);
void matrix_minus(double *,double *,int);
   //這兩個函數的計算結果保存在第一個矩陣中


//兩個矩陣相乘(重載2個函數）
void matrix_mul(double *,double *,double *,int,int,int);
void matrix_mul(double *,double *,double *,int);
    //第三個矩陣返回結果
void matrix_mul(double *,double *,int);//無法應用


//矩陣復制(重載2個函數)
void matrix_copy(double *,double *,int m,int n);
void matrix_copy(double *,double *,int n);
    //第二個矩陣返回結果


//矩陣求逆的全選主員高斯-約當法
int matrix_qiuni(double *,int n);
//對稱正定矩陣的求逆
int dzmatrix_qiuni(double *,int); 
//這兩個求逆函數還算比較準確，有一定的價值。我覺得。


//矩陣的三角分解(LU分解)
void matrix_fenjie(double *,double *,double *,int);
//一般實矩陣的QR分解
void matrix_qrfenjie(double *,double *,int,int);


//求Hessenberg矩陣全部特征根的QR法
int matrix_hqr(double *a,double *u,double *v,double eps,int n,int jt);
//該式子可以用來求高階方程的根

//約化一般實矩陣為Hessenberg矩陣的初等相似變換
void matrix_hessen(double *a,int n);












#endif