Utility

The script on this page requires the utility module tompy.py

* * *

Python Script

The following script solves a system of equations of the form Ax=b

linear_system.py

The key lines are

*from scipy import linalg*

*x = linalg.solve(A, b)*

* * *

Fortran

A Fortran program for this problem is: SLE.F

The program is compiled via

gfortran -o SLE SLE.F -llapack

It uses the Lapack subroutine DGESV.

* * *

Tom Irvine