#include<stdio.h> #include<stdlib.h> #include<time.h> #include<math.h> #define n 4 INT main() { float r[n][n+1]={ {8,1,3,-2,-3}, {-2,6,1,1,15}, {-1,5,9,2,-13}, {2,3,-3,9,70} }; float x[n] = {}; float c, max,d; float eps = 1e-5; INT iteracja = 1; INT i,k,j; srand(time(NULL)); printf("DANE TESTOWE:\n"); FOR(i = 0; i < n; i++) { FOR(j = 0; j <= n; j++) { printf("%f, ", r[i][j]); } printf("\n"); } printf("\n"); FOR(i=0;i<n;i++) { x[i]=0; c=r[i][i]; FOR(k=0;k<=n;k++) { r[i][k]=r[i][k]/c; } } do { max=0; FOR(i=0;i<n;i++) { c=r[i][n]; FOR(j=0;j<i;j++) { c=c-r[i][j]*x[j]; } FOR(j=i+1;j<n;j++) { c=c-r[i][j]*x[j]; } d = abs(x[i] - c); IF(d>max) { max=d; } x[i] = c; } printf("Iteracja %d\n",iteracja); FOR(i = 0; i < n; i++) { printf("X[%d] = %f ",i+1,x[i]); } printf("\n"); iteracja++; }WHILE(max>eps); printf("\nWynik:\n\n"); FOR(i = 0; i < n; i++) { printf("X[%d] = %f ",i+1,x[i]); printf("\n"); } RETURN 0; }