cpp2.h

高斯消去法解线性方程组

#include<stdlib.h>
#include<math.h>
#include<stdio.h>

int treede(int n,double a[],double b[])
{
	int k;

	for(k=0;k<=n-2;k++)
	{
		//归一化
		a[3*k+1]=a[3*k+1]/a[3*k]; 
		b[k]=b[k]/a[3*k];
		//消元法
		a[(k+1)*3]=a[(k+1)*3]-a[3*k+2]*a[3*k+1];
		b[k+1]=b[k+1]-b[k]*a[3*k+2];
	}

	b[n-1]=b[n-1]/a[(n-1)*3];
	for(k=n-2;k>=0;k--)
		b[k]=b[k]-a[3*k+1]*b[k+1];
	return 1;
}