北航数值分析计算大作业1

#include  
#include  
#define mat_num 501  //定义阶数
#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）;  //求两值较大的数
double min（double xdouble y）;  //求两值较小的数
double maxt（double xdouble ydouble z）; //求三个数中最大的
void creat_c（）;                          //转换为C阵
void move（double a_mindouble a_maxint k）; //平移函数
double fan_mingfa（double extend_num）;  //反幂法求解

double  a[501]u[501]y[501]C[5][501]detA;//定义全局变量


//主函数
void main（）
{
    int i;
	double num_1num_2num_3num_4cond_Atmpans[39];
	FILE *p;
	p = fopen（“c:\\Resul