#include?
using?namespace?std;

void?input（float**float*int）;
void?output（float*int）;
void?output1（float**int）;

int?main（）
{
	float?**p;
	float?*x*b*y;
????int??nijks;
	float?gh;

	cout<<“請(qǐng)輸入矩陣維數(shù)：“<	cin>>n;

????p=new?float*[n];??????????????????????//分配行指針的內(nèi)存
	for（i=0;i		p[i]=new?float[n];????????????????//分配每行元素的內(nèi)存

	x=new?float[n];
	b=new?float[n];
	y=new?float[n];

	input（pbn）;
????
	//?Doolittle三角分解
????for（k=0;k	{
	??for（j=k;j	??{
	????g=0;
		for（s=0;s		{
		??g=g+p[k][s]*p[s][j];
		}
		p[k][j]=p[k][j]-g;
	??}
????
	??for（i=k+1;i	??{
	????h=0;
		for（s=0;s		{
		??h=h+p[i][s]*p[s][k];
		}
		p[i][k]=（p[i][k]-h）/p[k][k];
	??}
	}

????cout<<“分解后矩陣為：“<	output1（pn）;

	//求解y
????y[0]=b[0];
	for（k=1;k	{
	??g=0;
	??for（j=0;j		??g=g+p[k][j]*y[j];
	??y[k]