用C语言实现雅可比迭代法_第1页
用C语言实现雅可比迭代法_第2页
用C语言实现雅可比迭代法_第3页
用C语言实现雅可比迭代法_第4页
全文预览已结束

下载本文档

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

文档简介

1、用 C 语言实现雅可比迭代法用 C 语言实现雅可比迭代法, 源程序代码如下 : #include "stdlib.h"#include "stdio.h"#include "conio.h"#include "string.h"#include "math.h"#define N 100float Table(int n,float aNN,float bN) int i,j;float cNN;printf("Please input the matrix A by row!n&qu

2、ot;); for(i=0;i<n;i+) printf("Row %d:",i+1);for(j=0;j<n;j+) scanf("%f",&aij);printf("Please input the vector b:"); for(i=0;i<n;i+) scanf("%f",&bi);for(i=0;i<n;i+) for(j=0;j<n;j+)if(i=j)cij=0; continue;cij=-aij/aii;printf("nThe matri

3、x A and vector b:n"); for(i=0;i<n;i+)for(j=0;j<n;j+) printf("%10.5f",aij);printf("%10.5f",bi); printf("n");printf("nThe Jacobi iterative scheme:n"); for(i=0;i<n;i+)for(j=0;j<n;j+) printf("%10.5f",cij);printf("%10.5f",bi/aii

4、); printf("n");float init_vec(int n,float xN)int i;printf("Please input the initial iteration vector x:"); for(i=0;i<n;i+)scanf("%f",&xi);printf("nThe initial iteration vector x:n"); for(i=0;i<n;i+)printf("%10.5f",xi);printf("n");

5、float jacobi(int n,float aNN,float bN,float xN)int i,j,k;float tmp,x2N;for(k=0;k+)for(i=0;i<n;i+)x2i=xi;for(i=0;i<n;i+)tmp=0.0; for(j=0;j<n;j+) if(j=i) continue; tmp+=aij*x2j; xi=(bi-tmp)/aii; for(i=0,j=0;i<n;i+) if(fabs(x2i-xi)<0.00001) j+;if(j=n) printf("nThis Jacobi iterative

6、scheme is convergent!n"); printf("Number of iterations: %d",k+1);break;if(k=499)printf("nThis Jacobi iterative scheme may be not convergent!"); break; printf("nThe results:n");for(i=0;i<n;i+) printf("%12.7f",xi); int main()int n;float xN,aNN,bN; printf("Input

温馨提示

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

评论

0/150

提交评论