追赶法解三对角线方程组（C++附注释）

数值分析/计算方法，欢迎交流

#include <iostream>
using namespace std;
double m[100][100],n[100];//m存放系数矩阵，n存放常数列
double a[100],b[100],c[100],B[100];//分别对应公式里的a,b,c,β;
double y[100],x[100];//公式里的y和x
void chase()
{
	int num;
	cout<<"请输入系数矩阵的阶数：";
	cin>>num;
	cout<<"请输入系数矩阵："<<endl;
	for(int i=0;i<num;i++)
	{
		for(int j=0;j<num;j++)
			cin>>m[i][j];
	}
	cout<<"请输入常数列："<<endl;
	for(int i=0;i<num;i++)
		cin>>n[i];
	for(int i=0;i<num;i++)
		b[i]=m[i][i];
	for(int i=0;i<num-1;i++)
		c[i]=m[i][i+1];
	for(int i=0;i<num-1;i++)
		a[i]=m[i+1][i];
	B[0]=c[0]/b[0];
	for(int i=1;i<num;i++)
		B[i]=c[i]/(b[i]-a[i]*B[i-1]);
	y[0]=n[0]/b[0];
	for(int i=1;i<num;i++)
		y[i]=(n[i]-a[i-1]*y[i-1])/(b[i]-a[i-1]*B[i-1]);
	x[num-1]=y[num-1];
	for(int i=num-2;i>=0;i--)
		x[i]=y[i]-B[i]*x[i+1];
	cout<<"解向量的转置为：";
	for(int i=0;i<num;i++)
		cout<<x[i]<<" ";
}
int main()
{
	chase();
	system("pause");
	return 0;
}