计算方法程序题计算方法程序题

西安交通大学计算方法B程序题。希望能够帮助你学习。

#include“stdio.h“
#include“math.h“
void Gauss（double D[]double x[]int n）//D[]盛放增广矩阵，n表示未知数的个数
{
	int ijkm;
	double t;
	for（i=0;i	{
		for（j=0;j		{
		    if（fabs（D[i*（n+1）+i]）			for（m=i;m<=n;m++）
				{
			    t=D[i*（n+1）+m];
		       	D[i*（n+1）+m]=D[（i+j+1）*（n+1）+m];
			    D[（i+j+1）*（n+1）+m]=t;
				}
		}		
		for（m=i;m		{
			t=D[（m+1）*（n+1）+i];
		    for（k=i;k<=n;k++）
			{    
			    
		      D[（m+1）*（n+1）+k]=D[（m+1）*（n+1）+k]+（D[i*（n+1）+k]/D[i*（n+1）+i]）*t*（-1）;
			}
		}
			
		
	}
	x[n-1]=D[（n-1）*（n+1）+n]/D[（n-1）*（n+1）+n-1];
	for（i=n-2;i>=0;i--）
	{	x[i]=D[i*（n+1）+n];
		for（j=i+1;j<=n-1;j++）
		x[i]=x[i]-D[i*（n+1）+j]*x[j];
	    x[i]=x[i]/D[i*（n+1）+i];
		
	}
	for（i=0;i<=n-1;i++）
	printf（“x[%d]=%.8lf\n“ix[i]）;
	
}
void main（）
{
	double x[4]={0};
	int n=4;
    double	D[20]={1.13483.83261.16513.40179.5342
	0.53011.78752.53301.54356.3941
	3.41294.93178.76431.314218.4231
	1.23714.999810.67210.014716.9237};
    Gauss（ Dxn）;
	
}

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1085  2010-12-13 16:15  3110001156\1.c

     文件        785  2010-12-13 21:20  3110001156\2.1.c

     文件        633  2010-12-13 21:21  3110001156\2.2.c

     文件        775  2010-12-13 21:24  3110001156\3.c

     文件        767  2010-12-13 16:56  3110001156\4.1.c

     文件        805  2010-12-13 16:58  3110001156\4.2.c

     文件       3734  2010-12-13 21:32  3110001156\5.c

     文件       3883  2010-12-13 20:50  3110001156\6.c

     文件        767  2010-12-13 21:30  3110001156\7.c

     文件       1900  2010-12-13 22:14  3110001156\8.c

     目录          0  2010-12-13 22:16  3110001156

----------- ---------  ---------- -----  ----

                15134                    11