计算方法实验数据.docxVIP

下载本文档

1
0
约1.02万字
约 23页
2019-02-21 发布于广东
举报
版权申诉

计算方法实验数据.docx

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

计算方法实验数据 Exp 1 #include stdio.h #define M 4 #define N 3 #define V 3 #define Y 2 int main() { FILE *fp1; float A[M][N],B[N][V],C[Y][M],D[M][N],E[M][N],F[Y][M]; int i,j,k; if((fp1=fopen(f1.txt,r))==NULL) { printf(error\n); return 0; } for(i=0; iM; i++) for(j=0; jV; j++) fscanf(fp1,%f,A[i][j]); FILE *fp2; if((fp2=fopen(f2.txt,r))==NULL) { printf(error\n); return 0; } for(i=0; iN; i++) for(j=0; jV; j++) fscanf(fp2,%f,B[i][j]); FILE *r1; r1=fopen(result1.txt,w); for(i=0; iM; i++) { for(j=0; jN; j++) { C[i][j]=0; for(k=0; kV; k++) C[i][j]=C[i][j]+A[i][k]*B[k][j]; printf(%10.6f,C[i][j]); fprintf(r1,%10.2f ,C[i][j]); } printf(\n); fprintf(r1,\n); } Exp 2 #include stdio.h #include math.h float D(double x) { return(sqrt(3+(1/x))); //return((3*x+1)/(x*x)); //return((x*x*x-1)/3); //return(1/(x*x-3)); } int main() { double x1,x2,x3; float eps=0.000000005; scanf(%f,x2); do{ x1=x2; x2=D(x1); printf(\n); printf(x=%10.6f\n,x2); }while(fabs(x1-x2)=eps); } Exp 3.1 #include iostream #include stdio.h int main() { FILE *fp1; fp1=fopen(f1.txt,r); FILE *fp2; fp2=fopen(f2.txt,r); double A[15][15],M[15][15]; double B[15],s,X[15]; int i,j,k=0,n=10; for(i=0; i=n-1; i++) { for(j=0; j=n-1; j++) fscanf(fp1,%lf,A[i][j]); fscanf(fp2,%lf,B[i]); } fclose(fp1); fclose(fp2); for(k=0; k=n-2; k++) { for(i=k+1; i=n-1; i++) { M[i][k]=A[i][k]/A[k][k]; for(j=k+1; j=n-1; j++) A[i][j]-=M[i][k]*A[k][j]; B[i]-=M[i][k]*B[k]; } } X[n-1]=B[n-1]/A[n-1][n-1]; for(i=n-2; i=0; i--) { s=0; for(j=i+1; j=n-1; j++) s