资源预览内容
第1页 / 共4页
第2页 / 共4页
第3页 / 共4页
第4页 / 共4页
亲,该文档总共4页全部预览完了,如果喜欢就下载吧!
资源描述
【这是自己编写的程序,请自己阅读着使用。】/*shuzhi1*/#include #define C 3void main()float ACC; float bC; float XC; float t,r,s; int i,j,k,m,n; printf(Please input A:n); for(i=0;i=C-1;i+)for(j=0;j=C-1;j+) scanf(%f,&Aij); printf(The data in the computer is:nA=n); for(i=0;i=C-1;i+)for(j=0;j=C-1;j+)printf(%4.4f ,Aij);printf(n); /*A数据的输入与检验*/printf(Is there any mistake? Press 1 for yes,Press 0 for no! );scanf(%d,&k);for(;k=1;)printf(Please input the adress(i,j) of the error data: );scanf(%d,&i);scanf(%d,&j);printf(Please input the right data: );scanf(%f,&t);Ai-1j-1=t;printf(The data in the computer is:nA=n); for(i=0;i=C-1;i+) for(j=0;j=C-1;j+) printf(%4.4f ,Aij); printf(n);printf(Is there any mistake? Press 1 for yes,Press 0 for no! ); scanf(%d,&k); /*修改出错数据*/printf(Please input b:n);for(i=0;i=C-1;i+)scanf(%f,&bi); printf(The data in the computer is:nb=n);for(i=0;i=C-1;i+)printf(%4.4f ,bi); /*b数据的输入与检验*/printf(nIs there any mistake? Press 1 for yes,Press 0 for no! );scanf(%d,&k);for(;k=1;)printf(Please input the adress(i) of the error data: );scanf(%d,&i);printf(Please input the right data: );scanf(%f,&t);bi-1=t;for(i=0;i=C-1;i+)printf(%4.4f ,bi); printf(nIs there any mistake? Press 1 for yes,Press 0 for no! ); scanf(%d,&k); /*修改出错数据*/ /*选择主元*/for(i=0;i=C-2;i+)for(j=i+1;j=C-1;j+)if(Aii=0)r=-Aii;elser=Aii;if(Aji=0)s=-Aji;elses=Aji;if (rs)for(k=0;k=C-1;k+) t=0; t=Ajk; Ajk=Aik; Aik=t; /*交换Ai与Aj*/ t=0; t=bj;bj=bi;bi=t; /*交换b*/ printf(After %d change,the the A and b are:nn,i+1);for(m=0;m=C-1;m+)for(n=0;n=C-1;n+)printf(%4.4f ,Amn); printf(b%d=%4.4f ,m+1,bm);printf(n);printf(n); /*打印变化后的数组*/ for(m=i+1;m=C-1;m+)t=(-Ami);for(n=0;n=C-1;n+) Amn=Amn+(Ain/Aii)*t;bm=bm+bi/Aii*t;for(m=0;m=C-1;m+)for(n=0;n=C-1;n+)printf(%4.4f ,Amn); printf(b%d=%4.4f ,m+1,bm);printf(n);printf(n); /*主元消零并打印当前数组*/ /*迭代求XC*/for(i=0;i=0;i-)t=0;for(j=i+1;j=C-1;j+) t=t+Aij*Xj;Xi=(bi-t)/Aii;printf(The answer of Ax=b is:n);for(i=0;i=C-1;i+) printf(X%d=%4.4f ,i,Xi);printf(n);/*程序结束*/
收藏 下载该资源
网站客服QQ:2055934822
金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号