matrix.h

数值计算中追赶法的C++实现算法

#define deltaT 100      //选取时间间隔deltaT=100秒
#define N 100          //将棒等分成N=100份
#define rou 7300		//常数ρ=7300	
#define K 30			//常数K=30
#define Crou 800		//常数Cρ=800
#define h 100			//常数h=100
#define l 0.1			//棒长l=0.1米
#define T0 1000			//棒的初始温度T0=1000度
#define T_infinite 20	//棒的最终温度T_infinite=20度

class A
{
protected:
	double *a;			//存放追赶法中的a
	double *b;			//存放追赶法中的b
	double *c;			//存放追赶法中的c
	double *f;			//存放追赶法中的f
	double *U;			//存放追赶法中计算得到的U
	double *L;			//存放追赶法中计算得到的L
	double *Y;			//存放追赶法中计算得到的Y
	double *X;			//存放追赶法中计算得到的X
	int n;				//存放追赶法中X的维数
public:
	
	A(int temp_n=2);    //构造函数，给成员变量分配空间
	int setA(double *temp_a,double *temp_b,double *temp_c,double *temp_f); //给a,b,c,f附值
	int getLU();		//计算出U和L
	int getY();			//计算出Y
	double getX();		//计算出X
	int changef();		//用新计算出的X替代上一次的f
	void show();		//将X写入文件
	
};