《列主元高斯消元法的C语言编程(3页).doc》由会员分享,可在线阅读,更多相关《列主元高斯消元法的C语言编程(3页).doc(3页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、-列主元高斯消元法的C语言编程-第 3 页列主元高斯消元法基本思想:用高斯消元法求解线性方程组时,为避免小的主元,在进行第步消元前,应该在第列元素中找出第一个出现的绝对值最大者,例如,再把第个方程与第个方程进行交换,使成为主元。我们称这个过程为选主元。由于只在第列元素中选主元,通常也称为按列选主元。列主元高斯消元法的C语言编程列主元高斯消元法的C语言程序代码如下:#include#include#include#define N 4void Gause_pivot(int n,double ANN+1,double x)/高斯消元int i,j,k; for(k=1;k=n-1;k+)/选主元
2、int ii,jj,kk,row;double max,temp;ii=1;kk=k;max=Aii-1+(kk-1)kk-1;for(jj=ii+1;jjfabs(max)max=Ajj-1+(kk-1)kk-1;row=jj+(kk-1);for(ii=1;ii=n+1;ii+)temp=Akk-1ii-1;Akk-1ii-1=Arow-1ii-1;Arow-1ii-1=temp; for(i=k+1;i=n;i+) for(j=k+1;j=1;k-)xk-1=0.0;for(j=k+1;j=n;j+) xk-1+=Ak-1j-1*xj-1;xk-1=(Ak-1n-xk-1)/Ak-1k-1; /在屏幕上输出结果for(i=0;i=n-1;i+) printf(%fn,xi);void main()double ANN+1=0,2,-0.1,7,76.2,6,0,3,-5,15,3,1,2,5,86,2,4,1,0,48; double xN=0;Gause_pivot(N,A,x);