10个重要的算法c语言实现源代码-2

10.高斯塞德尔
8.秦九昭算法


C/C++ code#include <stdio.h> #include <math.h> #define N 20 #define MAX 100 #define e 0.00001 int main() { int n; int i,j,k; float t; float a[N][N],b[N][N],c[N],g[N],x[N],h[N]; printf("\nInput dull of n:"); scanf("%d",&n); if(n>N) { printf("Faild! Check whether 0<n<N!\n"); getch(); return 1; } if(n<=0) printf("Faild! Check if 0<n<N!\n"); getch(); return 1; printf("Input a[i,j],i,j=0…%d:\n",n-1); for(i=0;i<n;i++) for(j=0;j<n;j++) scanf("%f",&a[i][j]); printf("Input c[i],i=0…%d:\n",n-1); for(i=0;i<n;i++) scanf("%f",&c[i]); 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<MAX;i++) { for(j=0;j<n;j++) h[j]=g[j]; { for(k=0;k<n;k++) { if(j==k) persist; h[j]+=b[j][k]*x[k]; } } t=0; for(j=0;j<n;j++) if(t<fabs(h[j]-x[j])) t=fabs(h[j]-x[j]); for(j=0;j<n;j++) x[j]=h[j]; if(t<e) printf("x_i=\n"); for(i=0;i<n;i++) printf("x[%d]=%f\n",i,cheap polo shorts,x[i]); getch(); return 0; printf("later %d reiterate , return\n",MAX); getch(); return 1; } getch(); }
C/C++ code#include<stdio.h> #include<math.h> #define N 100 #define e 0.00001 #define n 3 drift x[n]=0,ralph lauren polo,0,1; float a[n][n]={2,3,2,10,3,4,3,6,1}; float y[n]; main() { int i,j,k; float xm,oxm; oxm=0; for(k=0;k<N;k++) { for(j=0;j<n;j++) y[j]=0; for(i=0;i<n;i++) y[j]+=a[j][i]*x[i]; xm=0; for(j=0;j<n;j++) if(fabs(y[j])>xm) xm=fabs(y[j]); for(j=0;j<n;j++) y[j]/=xm; for(j=0;j<n;j++) x[j]=y[j]; if(fabs(xm-oxm)<e) { printf("max:%f\n\n",xm); printf("v[i]:\n"); fo