Energy industry

Development Platform:

C/C++

array_solution.cpp：Code Content

							// array_solution.cpp : Defines the entry point for the console application.
//
#include "stdafx.h"
#define n 3
void array_solution(float array[][n+1],float x[n])
{
	int i,j,k;
	for(i=0;i<n-1;i++)
	for(j=i+1;j<n;j++)
		for(k=j;k<n+1;k++)
		array[j][k]=array[j][k]-array[i][k]*(array[j][i]/array[i][i]);
	for(i=0;i<n;i++)
	{
		for(j=0;j<n+1;j++)
		{
			printf("%f",array[i][j]);
			printf("  ");
		}
		printf("n");
	}
	x[n-1]=array[n-1][n]/array[n-1][n-1];
	for(j=n-2;j>=0;j--)
	{
		float sum=0.0;
		for(i=n-1;i>j;i--)
			sum=sum+array[j][i]*x[i];
		x[j]=(array[j][n]-sum)/array[j][j];
	}
	for(i=0;i<n;i++)
		printf("x[%d]=%f  ",i,x[i]);
}
int main(int argc, char* argv[])
{
	float a[3][4]={{1,1,1,6},{1,1,-1,0},{-1,1,-1,-2}},x[n];
	array_solution(a,x);
	return 0;
}

						