数值计算方法实验报告(太原理工大学).docVIP

PAGE PAGE 16 《数值计算方法》 实验报告 学 院：太原理工大学现代科技学院 专业班级： 姓 名： 学 号： 指导教师： 时 间： 学院 现代科技学院 专 业 日期 2010年04月 学号 姓 名 成绩 实验题目 线性方程组高斯解法 指导老师 一．实验目的 熟练运用已学计算方法求解方程组。 加深对计算方法技巧，选择正确的计算方法来求解各种方程组。 培养使用电子计算机进行科学计算和解决问题的能力。 二．实验环境 DEV-C++ 实验语言：c语言 三．实验内容 1．合理利用Gauss消元法求解下列方程组： ① 四．实验公式 高斯分解法： lik=aik/akk aij= aij- lik* akj k=1,2,…,n-1 i=k+1,k+2, …,n j=k+1,k+2, …,n+1 xn= ann+1/ ann xk=( akn+1-∑akj xj)/ akk (k=n-1,n-2, …,2,1) 五．源程序 #include Stdio.h #include Conio.h #include stdlib.h #define MAX 20 float fArrx[MAX],fArra[MAX][MAX],fArrl[MAX][MAX],fArrt[MAX][MAX]; void welcome(); void input(int num); void output(int num,int iFlag); void ArratoArrt(int num); void printArrt(int num); void Gauss(int num); void printresult(int num,int iFlag); int main(void) { int num,iFlag; char cFlag; welcome(); printf(\n 请输入方程个数值2~10： ); iFlag=scanf(%d,num); fflush(stdin); while(num2||num10 || iFlag==0) { iFlag=1; printf(\n**** 输入有误，请重新输入 ****); printf(\n 请输入方程个数值2~10： ); iFlag=scanf(%d,num); fflush(stdin); } do { input(num); printf(\n_______________________________________________________________________________\n); output(num,iFlag); printf(\n_______________________________________________________________________________\n); printf(\n 确认请输入 Y 重新输入请输入 N : ); cFlag=getchar(); fflush(stdin); while(cFlag!=Y cFlag!=y cFlag!=N cFlag!=n) { printf(\n**** 输入有误，请重新输入 ****); printf(\n 确认请输入 Y 重新输入请输入 N : ); cFlag=getchar(); fflush(stdin); } }while(cFlag==N || cFlag==n); ArratoArrt(num); iFlag=100; Gauss(num); printresult(num,iFlag); system(pause); return 0; } void welcome() { system(cls); printf(数值计算方法：高斯消元法); printf(\n\nWELCOME); printf(\n); } void input(int num) { int i,j; printf(\n 请按行顺序输入线性方程组各数值：\n); for(i=1;i=num;i++) { for(j=1;j=num+1;j++) { scanf(%f,fArra[i][j]); } } fflush(stdin); } void output(int num,int iFlag) {