home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1993 #1 / NN_1993_1.iso / spool / comp / lang / fortran / 4960 < prev    next >
Encoding:
Text File  |  1993-01-10  |  1.8 KB  |  55 lines
  1. Newsgroups: comp.lang.fortran
  2. Path: sparky!uunet!zaphod.mps.ohio-state.edu!caen!destroyer!cs.ubc.ca!unixg.ubc.ca!kakwa.ucs.ualberta.ca!muskwa.ucs.ualberta.ca!dhalliwe
  3. From: userDHAL@mts.ucs.UAlberta.CA (David Halliwell)
  4. Subject: Re: simultaneous equations...solutions
  5. Message-ID: <dhalliwe.726710946@muskwa.ucs.ualberta.ca>
  6. Sender: news@kakwa.ucs.ualberta.ca
  7. Nntp-Posting-Host: muskwa.ucs.ualberta.ca
  8. Organization: University Of Alberta, Edmonton Canada
  9. References: <1993Jan10.191952.16384@nuscc.nus.sg>
  10. Date: Mon, 11 Jan 1993 00:09:06 GMT
  11. Lines: 42
  12.  
  13. eng10524@nusunix1.nus.sg (RAMASWAMI VISHWANATH) writes:
  14. >Dear nettors , 
  15. >    writing a program for  a fortran77 project , i have run into a
  16. >bad spot . the prg requires a matrix subroutine to calculate values for
  17. >a set of equations . but the equations are not in a square matrix form .
  18. >most books i have referred to give solutions to n x n matrices only . 
  19. >this prg of mine ,( to calcualte deflections of beams under different
  20. >loading conditions) uses finite differentiaiton to get a set of
  21. >equations . thus the equations are of the form 
  22. >        k = x + y + z
  23.  
  24. >and so on giving me a n x 3 matrix . i am yet to find any concrete
  25. >methods to solve such a system of equations . could anybody suggest
  26. >something ?
  27.  
  28.     Do you mean that you end up with an nx3 matrix that looks like:
  29.  
  30.     a1   b1   c1
  31.     a2   b2   c2
  32. ...
  33.     aN   bN   cN
  34.  
  35. or a matrix that looks like:
  36.  
  37.     b1   c1    0    0   ....
  38.     a2   b2   c2    0   ....
  39.      0   a3   b3   c3   ....
  40. ...
  41.      0    0  ...        0   aN   bN
  42.  
  43.  
  44. The second is typical of finite difference and finite element methods,
  45. and is referred to as a tri-diagonal matrix. Many numerical methods
  46. texts should give a selection of rapid solution techniques for tri-
  47. diagonal systems.
  48.  
  49.  
  50.  
  51. Dave Halliwell
  52. Department of Geography
  53. University of Alberta
  54. Edmonton, Alberta
  55.