home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 7 / 07.iso / c / c129 / 1.ddi / CMMATHDE.C < prev    next >
Encoding:
C/C++ Source or Header  |  1988-02-18  |  2.5 KB  |  97 lines
  1. # include <math.h>
  2. # include <stdio.h>
  3. # include "complex.h"
  4. # include "matmath.h"
  5. # include "miscio.h"
  6.  
  7. struct complext mat1[3][3];
  8. struct complext mat2[3][3];
  9. struct complext mat3[3][3];
  10. int order;
  11. float det;
  12. int i;
  13. int j;
  14. int k;
  15. char message[80];
  16. struct complext vm;
  17.  
  18.  
  19. void main()
  20. {
  21.    mat1[0][0].x = 13;
  22.    mat1[0][1].x = -8;
  23.    mat1[0][2].x = -3;
  24.    mat1[1][0].x = -8;
  25.    mat1[1][1].x = 10;
  26.    mat1[1][2].x = -1;
  27.    mat1[2][0].x = -3;
  28.    mat1[2][1].x = -1;
  29.    mat1[2][2].x = 11;
  30.  
  31.    mat1[0][0].y = 1;
  32.    mat1[0][1].y = 2;
  33.    mat1[0][2].y = 3;
  34.    mat1[1][0].y = 2;
  35.    mat1[1][1].y = 2;
  36.    mat1[1][2].y = 2;
  37.    mat1[2][0].y = 3;
  38.    mat1[2][1].y = 2;
  39.    mat1[2][2].y = 3;
  40.  
  41.    mat2[0][0].x = 1.0;
  42.    mat2[0][1].x = 5.0;
  43.    mat2[0][2].x = 1.0;
  44.    mat2[1][0].x = 3.0;
  45.    mat2[1][1].x = 6.0;
  46.    mat2[1][2].x = 3.0;
  47.    mat2[2][0].x = 0.0;
  48.    mat2[2][1].x = 1.0;
  49.    mat2[2][2].x = -9.0;
  50.  
  51.    mat2[0][0].y = 4.0;
  52.    mat2[0][1].y = 2.0;
  53.    mat2[0][2].y = 3.0;
  54.    mat2[1][0].y = 4.0;
  55.    mat2[1][1].y = 3.0;
  56.    mat2[1][2].y = 1.0;
  57.    mat2[2][0].y = 0.0;
  58.    mat2[2][1].y = 3.0;
  59.    mat2[2][2].y = -6.0;
  60.  
  61.    strcpy(message, "Original mat1 = ");
  62.    CMatPrint(&mat1[0][0],3,3,message);
  63.    printf("\n");
  64.    strcpy(message, "Original mat2 = ");
  65.    CMatPrint(&mat2[0][0],3,3,message);
  66.    printf("\n");
  67.    CMatProd(&mat1[0][0],&mat2[0][0],3,3,3,&mat3[0][0]);
  68.    strcpy(message, "matrix product mat1 x mat2 = ");
  69.    CMatPrint(&mat3[0][0],3,3,message);
  70.    printf( "press carriage return to advance to next procedure \n");
  71.    getch();
  72.    vm.x = 11;
  73.    vm.y = 0;
  74.    CMatScalarProd(&mat1[0][0],vm,3,3,&mat3[0][0]);
  75.    strcpy(message, "matrix scalar product mat1 x 11+ j0 = " );
  76.    CMatPrint(&mat3[0][0],3,3,message);
  77.    printf( "press carriage return to advance to next procedure \n");
  78.    getch();
  79.    CMatAdd(&mat1[0][0],&mat2[0][0],3,3,&mat3[0][0]);
  80.    strcpy(message, "matrix add mat1 + mat2 = ");
  81.    CMatPrint(&mat3[0][0],3,3,message);
  82.    printf("press carriage return to advance to next procedure \n");
  83.    getch();
  84.    CMatTranspose(&mat1[0][0],3,3,&mat3[0][0]);
  85.    strcpy(message, "matrix transpose of mat1  = ");
  86.    CMatPrint(&mat3[0][0],3,3,message);
  87.    printf( "press carriage return to advance to next procedure \n");
  88.    getch();
  89.    CMatInvert(&mat1[0][0],3,&det,&mat3[0][0]);
  90.    strcpy(message, "inverse of matrix mat1  = ");
  91.    CMatPrint(&mat3[0][0],3,3,message);
  92.    CMatProd(&mat1[0][0],&mat3[0][0],3,3,3,&mat2[0][0]);
  93.    strcpy(message, "matrix product mat1 x inverse(mat1) = ");
  94.    CMatPrint(&mat2[0][0],3,3,message);
  95. }
  96.  
  97.