home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 7 / 07.iso / c / c129 / 1.ddi / CGJDEMO.C < prev    next >
Encoding:
C/C++ Source or Header  |  1988-02-01  |  2.8 KB  |  114 lines
  1. # include <math.h>
  2. # include "complex.h"
  3. # include "cgj.h"
  4. # include "miscio.h"
  5. # include <ctype.h>
  6. # include <stdio.h>
  7.  
  8. struct complext xary[100];
  9. struct complext inv[100];
  10. struct complext yary[100];
  11. struct complext solution[100];
  12. int numvar;
  13. int i;
  14. int j;
  15. float matdet;
  16. char signr;
  17. char signi;
  18. char answer;
  19. struct complext com;
  20. float angle;
  21.  
  22.  
  23.  
  24.  
  25.  
  26. void GetData()
  27.  
  28. {
  29.   int j,i;
  30.  
  31.    ClrScr();
  32.    answer = 'n';
  33.    while ( answer == 'n' ) {
  34.       printf( "Enter number of equations - must equal number of variables: ");
  35.       scanf("%d", &numvar);
  36.       printf("\n");
  37.       printf("You will be asked to enter the complex coefficients.\n");
  38.       printf("Enter the real part first, followed by a space, then the\n");
  39.       printf( "coefficient of the imaginary part.\n");
  40.       for ( i = 0; i <= numvar-1; ++i ) {
  41.          printf("Equation # %  d\n",  i);
  42.          for ( j = 0; j <= numvar-1; ++j ) {
  43.             printf("Complex coefficient # %   d = ",  j );
  44.             scanf("%f %f", &xary[i*numvar+j].x,&xary[i*numvar+j].y);
  45.          }
  46.          printf("    constant = ");
  47.          scanf("%f %f", &yary[i].x, &yary[i].y);
  48.          printf("\n");
  49.       }
  50.       ClrScr();
  51.       printf( "The following system of equations was entered:\n");
  52.       printf("\n");
  53.       for ( i = 0; i <= numvar-1; ++i ) {
  54.          for ( j = 0; j <= numvar-1; ++j ) {
  55.             printf("(%5.2f %5.2fi X%1d)",  xary[i*numvar+j].x,
  56.                      xary[i*numvar+j].y, j );
  57.             if ( j != numvar ) {
  58.                printf( " + ");
  59.             }
  60.          }
  61.          printf(" = (%5.2f %5.2fi)",  yary[i].x, yary[i].y);
  62.          printf("\n");
  63.       }
  64.       printf( "is this right (y OR n) ?  ");
  65.       answer = getch();
  66.       answer = toupper(answer);
  67.    }
  68. }
  69.  
  70. void CheckAnswer()
  71. {
  72.    struct complext sum[100];
  73.    struct complext temp;
  74.    int i;
  75.    int j;
  76.  
  77.    printf( "Checking the answer by calculating y values using calculated\n");
  78.    printf( "x values and original coefficients.\n");
  79.    for ( i = 0; i <= numvar-1; ++i ) {
  80.       sum[i].x = 0;
  81.       sum[i].y = 0;
  82.       for ( j = 0; j <= numvar-1; ++j ) {
  83.          ComplexMath(xary[i*numvar+j],'*',solution[j],&temp);
  84.          ComplexMath(sum[i],'+',temp,&sum[i]);
  85.       }
  86.       printf("y%2d = ",  i );
  87.       printf("(%6.2f  %6.2fj )", sum[i].x, sum[i].y);
  88.       printf("\n");
  89.    }
  90. }
  91.  
  92.  
  93. void main()
  94.  
  95. {
  96.    ClrScr();
  97.    GetData();
  98.    printf("\n");
  99.    printf("\n");
  100.    ComplexGaussJordan(xary,yary,numvar,solution,inv,&matdet);
  101.   
  102.    printf( "The solution is: \n");
  103.    printf("\n");
  104.    for ( i = 0; i <= numvar-1; ++i ) {
  105.       printf("x %1d = %6.3f %6.3fj\n", i, solution[i].x, solution[i].y);
  106.    }
  107.    printf("\n");
  108.    printf("matrix determinant = %6.3f\n",  matdet);
  109.    CheckAnswer();
  110.  
  111.  
  112.  
  113. }
  114.