雅克比迭代法 数值分析

#include <stdio.h> #include <math.h> #define MAX_N 20 #define MAXREPT 20 #define epsilon 0.00001 int main() { int n; int i,j,k,h; double err; printf(" 雅克比迭代法\n"); look: static double a[MAX_N][MAX_N],b[MAX_N][MAX_N],c[MAX_N],g[MAX_N]; static double x[MAX_N],nx[MAX_N]; printf("输入维数n:"); scanf("%d",&n); if(n>MAX_N) { printf("输入的维数n大于MAX_N,请重新定义MAX_N!!\n"); goto look; } if(n<=0) { printf("请输入一个数在1和%d 之间,否则不合题意!\n",MAX_N); goto look; } printf("请输入系数矩阵a(i,j):\n"); for(i=0;i<n;i++) for(j=0;j<n;j++) scanf("%lf",&a[i][j]); printf("请输入b[i]:\n"); for(i=0;i<n;i++) scanf("%lf",&c[i]); printf("\n"); printf(" k X(k) ||X(k)-X(k-1)||\n"); for(i=0;i<n;i++) for(j=0;j<n;j++) { b[i][j]=-a[i][j]/a[i][i]; g[i]=c[i]/a[i][i]; } for(i=0;i<MAXREPT;i++) { printf("%d ",i+1); for(j=0;j<n;j++) nx[j]=g[j]; for(j=0;j<n;j++) { for(k=0;k<n;k++) { if(j==k) continue; nx[j]+=b[j][k]*x[k]; } } printf("( "); for(j=0;j<n;j++) printf(" %lf",nx[j]); printf(" )"); err=0; for(j=0;j<n;j++) { if(err<fabs(nx[j]-x[j])) err=fabs(nx[j]-x[j]); } printf(" %lf\n",err); for(j=0;j<n;j++) x[j]=nx[j]; if(err<epsilon) { printf("方程的解是:\n"); for(i=0;i<n;i++) printf("X%d=%lf\n",i+1,x[i]); printf("是否继续?是(1),否(0):"); scanf("%d",&h); for(i=0;i<n;i++) x[i]=0; switch(h) { case 1:goto look; case 0:break; } return 0; } } printf("在MAX_N的定义下,得不出结果,请重新定义!!!\n"); return 1; }