home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 7 / 07.iso / c / c129 / 1.ddi / GJDEMO.C < prev    next >
Encoding:
C/C++ Source or Header  |  1988-01-19  |  1.7 KB  |  80 lines
  1. # include <stdio.h>
  2. # include <math.h>
  3. # include "gj.h"
  4.  
  5.  
  6. float x[100];
  7. float inv[100];
  8. float y[10];
  9. float solution[10];
  10. int numvar;
  11. int i;
  12. int j;
  13. float matdet;
  14. char sign;
  15. char answer;
  16.  
  17.  
  18.  
  19. void GetData()
  20. {
  21. float temp;
  22.    ClrScr();
  23.    answer = 'N';
  24.    while ( answer == 'N' ) {
  25.       printf("enter number of equations - must equal number of variables: ");
  26.       scanf("%d", &numvar);
  27.       printf("\n");
  28.       for ( i = 0; i <= numvar-1; ++i ) {
  29.          printf("equation # %2d\n", i);
  30.          for ( j = 0; j <= numvar-1; ++j ) {
  31.             printf("coefficient # %2d = ",  j);
  32.             scanf("%f", &x[i*numvar+j]);
  33.          }
  34.          printf("constant = ");
  35.          scanf("%f", &y[i]);
  36.          printf("\n");
  37.       }
  38.       ClrScr();
  39.       printf("the following system of equations was entered:\n");
  40.       printf("\n");
  41.       for ( i = 0; i <= numvar-1; ++i ) {
  42.          for ( j = 0; j <= numvar-1; ++j ) {
  43.             if ( x[i*numvar+j] < 0 ) {
  44.                sign = '-';
  45.             }
  46.             else {
  47.                sign = '+';
  48.             }
  49.             temp = fabs(x[i*numvar+j]);
  50.             printf("  %c %8.2f X%2d", sign, temp, j);
  51.          }
  52.          printf(" = %8.2f", y[i]);
  53.          printf("\n");
  54.       }
  55.       printf( "is this right (y OR n) ?  ");
  56.       answer = getche();
  57.       printf("\n");
  58.       answer = toupper(answer);
  59.  
  60.    }
  61. }
  62.  
  63.  
  64. void main()
  65. {
  66.    GetData();
  67.    printf("\n");
  68.    printf("\n");
  69.    GaussJordan(x,y,numvar,solution,inv,&matdet);
  70.    printf( "the solution is: \n");
  71.    printf("\n");
  72.    for ( i = 0; i <= numvar-1; ++i ) {
  73.       printf("x %1d = %8.3f\n", i,  solution[i]);
  74.    }
  75.    printf("\n");
  76.    printf("matrix determinant = %f\n",  matdet);
  77.  
  78. }
  79.  
  80.