数值分析-c语言实现高斯赛德尔解法(共6页).doc

《数值分析-c语言实现高斯赛德尔解法(共6页).doc》由会员分享，可在线阅读，更多相关《数值分析-c语言实现高斯赛德尔解法(共6页).doc（6页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上 数值分析 程 序 设 计 学院：计算机学院 姓名：袁薪洋 1 实验目的： 1熟练掌握C语言程序设计，编程求解问题。 2.运用高斯-赛德尔迭代公式求解线性方程组。2 实验内容： 用高斯-赛德尔迭代公式求解方程组。 10x1-x2-2x3=7.2-x1+10x2-2x3=8.3 -x1-x2+5x3=4.2 程序的核心代码：#includemath.h#include#define NUMBER 20float ANUMBERNUMBER+1 ;float ark;int flag,n;void exchange(int r,int k);float max(int k

2、);void main() float xNUMBER; /*此数组用于存放方程解*/ int r,k,i,j; printf(*); printf(nn用高斯-赛德尔迭代法解线性方程组nn); printf(*); printf(nn 请输入方程组的维数：n=); scanf(%d,&n); printf( nn请输入系数矩阵A和向量b:); for(i=1;i=n;i+) printf(nn请输入a%d1-a%d%d系数和向量b%d(数之间用空格格开):,i,i,n,i); /实现将每一行中的系数和向量一次性输入，数之间用空格格开，输完后回车确定 for(j=1;j=n+1;j+) /将刚

3、才输入的数存入数组 scanf(%f,&Aij); for(k=1;k=n-1;k+) ark=max(k); if(ark=0) /判断方程是否为线性方程 printf(nn此方程组不合法!); else if(flag!=k) exchange(flag,k); for(i=k+1;i=n;i+) for(j=k+1;j=1;k-) float me=0; for(j=k+1;j=n;j+) me=me+Akj*xj; xk=(Akn+1-me)/Akk; for(i=1;i=n;i+) printf( nnx%d=%f,i,xi); void exchange(int r,int k) /交换行的矩函数 int i; for(i=1;i=n+1;i+) A0i=Ari; for(i=1;i=n+1;i+) Ari=Aki; for(i=1;i=n+1;i+) Aki=A0i;float max(int k) /比校系数大小的函数 int i; float temp=0; for(i=k;itemp) temp=fabs(Aik); flag=i; return temp;3. 运行结果：截图如下： 专心-专注-专业