matrix.h

求矩阵A i,j)=1/ i+j+1)的最大、最小特征值

#ifndef _MATRIX_H
#define _MATRIX_H

enum Fill_Way{crs,row,col,all};

class Matrix{					//仅为方阵!! 
private:
	int dimension;
	int cross; 
	double **elements;
	
public: 
	Matrix(int d=1);								//仅指定维数			
	Matrix(int,double **);						//用一个二维数组来填充矩阵 
	Matrix(int d,double e,int c=0);			//用一个数来填充矩阵,所有元素都=这个数 
	Matrix(int d,double *&es,Fill_Way fw=all);	//对角则用一个数组来填充对角,否则 
	Matrix(int d,double *es,Fill_Way fw=all);
	Matrix(const Matrix &,const Matrix &,const Matrix &,const Matrix &);
	~Matrix(void);
	int getDimension(void);
	double getElement(int,int);
	Matrix operator*(const Matrix &);		//矩阵乘法 
	Matrix operator*(double &);				//数乘 
	Matrix operator-(Matrix &);
	double *operator*(double *);
	Matrix operator/(double);
	int operator==(Matrix &);		//重载== 
	void LU(Matrix *&,Matrix *&);	//LU分解
	void exRow(int,int);			//交换两行 
	Matrix trans(void);
	int getCross(void);
	void setCross(void);
	void setElement(double,int,int);
	double *getRow(int);
	double *getCol(int);
	double *getCrossE(void);
	void setRow(double *,int);
	void setCol(double *,int);
	void setCrossE(double *);
	void print(int);
	void print(void);
};
#endif