反幂法C语言

//反幂法第一个分数矩阵，先输入系数矩阵的分子,再输入初始向量
#include<stdio.h>
#include<math.h>
#define n 5                     //定义系数矩阵行数                    
int main(){

	double a[n][n]={0};                     //初始化矩阵数组
	int i,j,k;
	int f=0;
	double L[n][n]={0};
	double U[n][n]={0};
	double X[n]={0};
	double Y[n]={0};
	double y[n]={0};
	double Z[n]={0};
	double c;
	double b1=1.000000;
	double b2=0.000000;
	//输入系数矩阵的分母
    printf("Please enter the denominator of the matrix:\n");
	printf("(The denominator of 1 is 1)\n");
	for (i=0; i<n;i++){
		for(j=0;j<n;j++){
			scanf("%lf",&c);
			a[i][j]=1.0/c;
		}     	
	}
	//同时寻找LU中的待定元素
	for(k=0;k<n;k++){	
		for(i=k;i<n;i++){
			for(j=0;j<k;j++){
					U[k][i]=U[k][i]-L[k][j]*U[j][i];
	    	        L[i][k]=L[i][k]-L[i][j]*U[j][k];	
			}
		U[k][i]=U[k][i]+a[k][i];
		L[i][k]=(L[i][k