用C语言编写程序：N元一次方程组的解.docx_第1页
用C语言编写程序：N元一次方程组的解.docx_第2页
用C语言编写程序：N元一次方程组的解.docx_第3页
用C语言编写程序：N元一次方程组的解.docx_第4页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

/*2014年5月20日 12:23:20功能：已知有一三元一次方程组如下,方程系数由键盘输入，求其3个根。*/# include int main(void)int i, j, m, n, r, k=0, t; /i,j,m,k,t 作为循环计数变量；n表示有n个未知数，即n元；r表示有r个方程组成方程组。/n = r = 3;printf(请输入未知数个数n和方程的个数:n);scanf(%d %d, &n, &r); /n表示有n个未知数，即n元；r表示有r个方程组成方程组。float a1120; /用二维数组存储方程组参数。11和100课改动。这里的程序限 10 元 20 个方程的方程组。float b200; /原来存放之间变量。printf(请输入方程组参数n);for(j=0; jn; j+)for(i=0; ir+1; i+)scanf(%f, &aji);/*for(j=0; jn; j+) /此段程序作为检测该程序正确性的试数作用。for(i=0; in+1; i+)printf(%f , aji);printf(n);float a1020=2,2,-1,6,1,-2,4,3,5,7,1,28;*/if(rn)printf(方程有无穷解n);if(r=n)printf(方程有唯一解n);for(t=0; tn-1; t+) /循环次数：当t=0时，仅将a10,a20,a21的值转换为0；t=1时，才将其转换成最简型阶梯矩阵。 for(m=0; mn; m+) /m表示第m列for(j=0; jn; j+) /j表示第j行if(j!=m)bk = ajm/amm;for(i=0; in+1; i+) /因为一行的每个元素均要进行运算等式两边才会相等。aji -= ami*bk;bk=0; /此处是必不可少的。因为内存中会有垃圾值，必须要处理。k+;for(j=0; jn; j+) /为了将其转换成单位矩阵，见线性代数第98页。ajn /= ajj;ajj /=ajj;for(j=0; jr; j+)printf(x%d = %fn, j, ajn); /输出最终结果，即未知数x,y,z.的值。for(j=0; jn; j+)for(i=0; in+1; i+)printf(%f , aji);printf(n);return 0;/*输入为：2 2 -1 6 1 -2 4 3 5 7 1 28即方程组为：2x+2y- z=6 x-2y+4z=35x+7y+ z=28其正确结果是：x=1,y=3,z=2在VC+6.0中的输出结果是：请输入未知数个数n和方程的个数:3 3请输入方程组参数2 2 -1 6 1 -2 4 3 5 7 1 28方程有唯一解x0 = 1.000000x1 = 3.000000x2 = 2.0000001.000000 -0.000000 0.000000 1.0000000.00000
人人文库> 全部分类> 教育资料 > 中学教育

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论