home *** CD-ROM | disk | FTP | other *** search
/ Power-Programmierung / CD2.mdf / c / library / dos / diverses / leda / prog / basic / matrix.c < prev    next >
Encoding:
C/C++ Source or Header  |  1991-11-15  |  997 b   |  64 lines
  1. #include <LEDA/vector.h>
  2. #include <LEDA/matrix.h>
  3. #include <LEDA/list.h>
  4.  
  5. declare(list,vector)
  6.  
  7.  
  8. main()
  9.  
  10. { 
  11.   int d = read_int("dimension = ");
  12.  
  13.   matrix       A(d,d);
  14.   vector       b(d);
  15.   vector       x(d);
  16.  
  17.   list(vector) B;  // List of b-vectors
  18.   list(vector) X;  // List of x-vectors
  19.  
  20.   printf("give (%d x %d) - matrix A:\n",d,d);
  21.  
  22.   cin >> A;
  23.  
  24.   cout << "A = \n";
  25.   cout << A << "\n";
  26.   newline;
  27.  
  28.   cout << "A.trans() = \n";
  29.   cout << A.trans() << "\n";
  30.   newline;
  31.  
  32.  
  33.   cout << "A.inv() = \n";
  34.   cout << A.inv() << "\n";
  35.   newline;
  36.  
  37.   cout << "A*A.inv() = \n";
  38.   cout << A*A.inv() << "\n";
  39.   newline;
  40.   
  41.  
  42.   cout << "We solve the system of linear equations A*x = b.\n";
  43.   cout << "Give a list of vectors b (terminated by ctrl-d):\n";
  44.  
  45.   while (cin >> b) B.append(b);
  46.   newline;
  47.  
  48.   B.print("b's = ",'\n');
  49.   newline;
  50.   newline;
  51.  
  52.   forall(b,B)  X.append(A.solve(b));
  53.   
  54.   X.print("x's = ",'\n');
  55.   newline;
  56.   newline;
  57.  
  58.  
  59.   cout << "A*x's = \n";
  60.   forall(x,X) cout << A*x << "\n";
  61.   newline;
  62.  
  63. }
  64.