华北电力大学科技学院数值计算方法上机报告

《华北电力大学科技学院数值计算方法上机报告》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、华北电力大学科技学院农电12K1G&P华北电力大学科技学院数值计算方法上机报告姓名：G&P班级：农电12K1学号：2014年11月28日华北电力大学科技学院农电12K1G&P华北电力大学科技学院农电12K1G&P华北电力大学科技学院农电12K1G&P华北电力大学科技学院农电12K1G&P华北电力大学科技学院农电12K1G&P华北电力大学科技学院农电12K1G&P华北电力大学科技学院农电12K1G&P华北电力大学科技学院农电12K1G&P华北电力大学科技学院农电12K1G&P华北电力大学科技学院农电12K1G&P华北电力大学科技学院农电12K1G&P华北电力大学科技学院农电12K1G&P华

2、北电力大学科技学院农电12K1G&P华北电力大学科技学院农电12K1G&P华北电力大学科技学院农电12K1G&P华北电力大学科技学院农电12K1G&P华北电力大学科技学院农电12K1G&P华北电力大学科技学院农电12K1G&P实验程序为实验一#include"stdio.h"#include"math.h"#defineN30voidmain(){inti;floatx[N],c;scanf("%f",&c);scanf("%f",&x[0]);//初值printf("****倒数表****");for(i=0;i

3、2*c*x[i]-1);printf("k=%dtx(%d)=%.5f",i,i,x[i]);if(fabs(x[i+1]-x[i])<=0.0005)break;else华北电力大学科技学院农电12K1G&Pcontinue;printf("k=%dtx(%d)=%.5f",i+1,i+1,x[i]);printf("计算结果：1/%f=%.3f",c,x[i+1]);}}实验二#include"stdio.h"#include"math.h"#defineN10voidmain(){inti,j,k,l,n;floatb[N],a[N][N],t,d,de

4、t=1.0;//***数据输入*/printf("输入n");scanf("%d",&n);printf("输入Ai");for(i=0;i

5、k++)//从第一次消元到第N-1次消元{d=a[k][k];l=k;for(i=k+1;ifabs(d)){d=a[i][k];i=l;}}if(i==n)//判断是否奇异，不奇异进行行交换{if(d==0)printf("奇异");//如果所有行的"首列"都为0，为奇异else{if(l!=k)//如果第k行的"首列"并不是最大{det=det*(-1);华北电力大学科技学院农电12K1G&Pfor(j=k;j<=n;j++)//交换系数矩阵中的两行{t=a[l][j];a[l][j]=a[k]

6、[j];a[k][j]=t;}t=b[l];b[l]=b[k];b[k]=t;//交换右端常向量中的两行}}}//************列选主元函数*/for(i=k+1;i

7、b[n-1]/a[n-1][n-1];//计算x(N)的解for(i=n-2;i>=0;i--)//从倒数第二项开始依次回代N-1次{t=0;for(j=i+1;j