gauss.cpp

电路分析小程序

float* *GaussFE(float* *Matrixg,int nViria)
 {
	float b;
	for(int k = 1; k <=nViria; ++k)
	{
		b = Matrixg[k][k];
		for(int j = k; j <=nViria  + 1; ++j)
		{
			Matrixg[k][j] = Matrixg[k][j] / b;
		}

		for(int i = k + 1; i <= nViria; ++i)
		{
			b = Matrixg[i][k];
			for(j = k; j <=nViria +1; ++j)
			{
				Matrixg[i][j] =  Matrixg[i][j] - b * Matrixg[k][j];
			}
		}
	}
	return Matrixg;
 }


float* GaussBS(float** array,int number)
{
	float*p=new float[number+1];
	for(int i=number;i>0;i--)
	{
		p[i]=array[i][number+1];
		for(int j=i+1;j<=number;j++)
		{
			p[i]=p[i]-array[i][j]*p[j];
		}
	}
	return p;
}