矩阵lu分解.cpp

这是一种牛顿迭代算法

#include<stdio.h>
#include<math.h>

void main()
{
 int i,j,k,r=0;
 int n=10;
 int b[10],y[10],x[10];
 int a[10][10];
 int l[10][10];
 int u[10][10];
 printf("input a[n][n]:\n");
 for(i=0;i<n;i++)
	 for(j=0;j<n;j++)
		 scanf("%d",&n);
	 for(i=0;i<n;i++)
		 a[0][i]=u[0][i];
	 for(i=1;i<n;i++)
		 l[i][0]=a[i][0]/u[0][0];
	 for(k=0;k<r-1;k++)
	 { j=0;
	   j=j+l[r][k]*u[k][i];
	  for(i=r;i<n;i++)
		u[r][i]=a[r][i]-j;
	  for(k=0;k<r-1;k++)
	  { j=0;
	    j=j+l[i][k]*u[k][r];
		for(i=r+1;i<n;i++)
			l[i][r]=(a[i][r]-j)/u[r][r];
	  }
	 }
 for(k=1;k<i-1;k++)
 {j=0;j=j+l[i][k]*y[k];
 for(i=1;i<n;i++)
	 y[i]=b[i]-j;
 for(k=i+1;k<n;k++)
 {j=0;j=j+u[i][k]*x[k];
 for(i=n-1;i>0;i--)
	 x[i]=(y[i]-j)/u[i][i];
 }
 }
}