trid.h

利用gauss分解求解线性方程组

//解三对角方程组的追赶法
int trid(float a[],float b[],float c[],float f[])
{
	float s[n],t[n],temp;
	int i;
	s[0]=0;t[0]=1;
	for(i=0;i<=n-2;i++)
	{   
		temp=b[i+1]+a[i+1]*s[i];
		if(fabs(temp)<eps)
		{ 
			printf("The equation has no solution!");
			return(-1);
		}
		else
		{   
			s[i+1]=-c[i+1]/temp;
			t[i+1]=(f[i+1]-a[i+1]*t[i])/temp;
		}
	}
	f[n-1]=t[n-1];
	for(i=n-2;i>=0;i--)
       f[i]=s[i]*f[i+1]+t[i];
	return(0);
}