数值分析 各个算法的C语言程序
double A[3][3]={5,2,1,-1,4,2,2,-3,10},x[50][3],b[3]={-12,20,3},f1=0,f2=0,L=1; int i,j,k=1,n=3;
for (i=0;i<3;i++) x[0][i]=0; while(L>0.00001||L<-0.00001) {
for(i=0;i<n;i++)
{ for(j=0;j<n;j++)
{if(j<i) f1+=A[i][j]*x[k][j];
else if(j>i) f2+=A[i][j]*x[k-1][j];} x[k][i]=(b[i]-f1-f2)/A[i][i];}
L=x[k][0]-x[k-1][0];
for(i=1;i<n;i++)
if((x[k][i]-x[k-1][i])/(x[k][i-1]-x[k-1][i-1])>1||(x[k][i]-x[k-1][i])/(x[k][i-1]-x[k-1][i-1])<-1) L=x[k][i]-x[k-1][i]; k++;
}
printf("输出迭代X矩阵:\n"); for(i=1;i<k;i++)
{for(j=0;j<n;j++) printf("%f ",x[i][j]); printf("\n");} printf("\n输出迭代次数:%d\n",k-1);
}
十二、 超松弛迭代法
#include<stdio.h> #include<stdlib.h> #include<math.h> void main() {
百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读,免费范文网,提供经典小说教育文库数值分析 算法C语言程序(10)在线全文阅读。
相关推荐: