#include??
#include??
#define?mat_num?501??//定義階數(shù)
#define?b?0.16?????//定義b
#define?c?-0.064???//定義c

double?mingfa（double?extend_num）;???//冪法求解
void?chushihua（）;?????????????????????//矩陣初始化
void?LU_resolve（double?p[5][501]）;????//LU矩陣分解
void?LU_solve（double?p[5][501]double?B[501]double?d[501]）;//LU求解
double?max（double?xdouble?y）;??//求兩值較大的數(shù)
double?min（double?xdouble?y）;??//求兩值較小的數(shù)
double?maxt（double?xdouble?ydouble?z）;?//求三個(gè)數(shù)中最大的
void?creat_c（）;??????????????????????????//轉(zhuǎn)換為C陣
void?move（double?a_mindouble?a_maxint?k）;?//平移函數(shù)
double?fan_mingfa（double?extend_num）;??//反冪法求解

double??a[501]u[501]y[501]C[5][501]detA;//定義全局變量


//主函數(shù)
void?main（）
{
????int?i;
	double?num_1num_2num_3num_4cond_Atmpans[39];
	FILE?*p;
	p?=?fopen（“c:\\Resul