矩阵的三角分解c程序

#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<<“请输入矩阵维数：“<	cin>>n;

    p=new float*[n];                      //分配行指针的内存
	for（i=0;i		p[i]=new float[n];                //分配每行元素的内存

	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]