r.h

数值计算利用LU分解法来求解线性方程组,简明方便

void R_XY(float *X,int M,float *H,int N,float *R)
{   
	int k,i;
	for(k=0;k<M;k++)
	{	*(R+k)=0.0;
	    for(i=0;i<N;i++)
		{
			if((k+i)<M)
				*(R+k)=*(R+k)+(*(X+i))*(*(H+k+i));
		}
	}
}
void R_XX(float *X ,int M,float *R)
{   
	int k,i;
	for(k=0;k<M;k++)
	{	*(R+k)=0.0;
	    for(i=0;i<M;i++)
		{
			if((k+i)<M)
				*(R+k)=*(R+k)+(*(X+i))*(*(X+k+i));
		}
	}
}