home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #16 / NN_1992_16.iso / spool / sci / math / symbolic / 2116 < prev    next >
Encoding:
Internet Message Format  |  1992-07-30  |  1.7 KB
  1. Path: sparky!uunet!mcsun!sunic!lunic!eru.mt.luth.se!wucm
  2. From: wucm@mt.luth.se (Wu Changming)
  3. Newsgroups: sci.math.symbolic
  4. Subject: Help needed with numerical inverse of a matrix in Maple
  5. Message-ID: <1985@eru.mt.luth.se>
  6. Date: 30 Jul 92 14:41:00 GMT
  7. Reply-To: Wu Changming <wucm@mt.luth.se>
  8. Organization: University of Lulea, Sweden
  9. Lines: 50
  10. UUCP-Path: {uunet,mcvax}!sunic!mt.luth.se!wucm
  11.  
  12.  
  13. I am using Maple V, but I do not know much about Maple, so I need some help here.
  14.  
  15. I need to do an inverse of a 8x8 matrix. Since the computer never gives the
  16. result, I plan to do the inverse numerically. However I got the error message:
  17.  
  18. ##################################################
  19. Error, (in normal/doit) floats not handled
  20. --------------------------------------------
  21.  
  22. Then I use the following to do a test.
  23.  
  24. ##################################################
  25.  K := array(1..4, 1..4):
  26.  for ii from 1 by 1 to 4 do
  27.   for jj from 1 by 1 to 4 do
  28.     K[ii, jj] := Pi**(1/ii*jj) + I**(ii+jj) :
  29.    od;
  30.  od;
  31.  
  32.  IK := inverse(K):
  33.  IK :=  map(evalc, map(evalf, map(evalm, IK ) ) ):
  34.  print(IK);
  35. ---------------------------------------------------  
  36.  
  37. which gives the symbolic solution first, and then the numerical solution.
  38.  
  39. The I do the above calculation in this way:
  40.  
  41. ##################################################
  42.  K := array(1..4, 1..4):
  43.  for ii from 1 by 1 to 4 do
  44.   for jj from 1 by 1 to 4 do
  45.     K[ii, jj] := evalc( evalf( Pi**(1/ii*jj) + I**(ii+jj) )):
  46.    od;
  47.  od;
  48.  
  49.  IK := inverse(K) :
  50.  print(IK)
  51. -----------------------------------------------------
  52.  
  53.  
  54. The program gives  `Error, (in normal/doit) floats not handled' when it is doing
  55. the inverse.
  56.  
  57. Could someone help me to do numerical inverse of a large matrix?
  58.  
  59. Thanks
  60.  
  61. Wu Changming
  62.