void displaymatrix(float* ptr,int order);
	void LU_decompose(float matrix_a[][500],float*matrix_L,float*matrix_U,float*vector_b,int order);
	void solve_L(float*matrix,float*vector,float*ans,int order);
	void solve_U(float*matrix,float*vector,float*ans,int order);
	void gauss(float*matrix,float*vector,int n,int order);
	float absolute(float num);