Solving System of Linear Equations

[ahad29]

Is there a simple command in Fortran for solving system of linear equations. Like matlab has x = A\B to solve Ax=B. Thanks

Ahad

 

[IRstuff]

no

TTFN

 

[johnhors]

Ahad

As IRstuff said the answer is "NO" - you have to write the source yourself.

But just about every introductory book on Fortran uses the Gaussian elemination to solve linear equations in their text as an example. Otherwise a Google search throws up sites like

http://library.lanl.gov/numerical/bookfpdf.html

which gives complete fortran listings of various methods to solve this rather trivial problem.

 

[nbucska]

Digital Equipment Corp. had a Scientific Subroutine Package
(SSP) written in F. and I am sure it has what you need.


<nbucska@pcperipherals DOT com> subj: eng-tips
read FAQ240-1032

 

[xwb]

If you want something faster than Gaussian, look up Cholesky/Choleski LU Decomposition.

 

[IRstuff]

If you want to do it yourself:

http://www.library.cornell.edu/nr/

TTFN