A Fortran 90 subroutine to calculate an inverse of a matrix using a subroutine for Gaussian elimination with backwards substitution:
Fortran 90 Uniform Random Number Generators
NOTE: The Numerical Recipes demo and tutorial are now new and improved!!
There were several glaring errors before the improvements.
Here is the silly little program I gave in class
to demonstrate passing a function name into a subroutine.
Here is that same program with all external
subroutines and functions.
In addition, you should check this web page for when the INTERFACE statement is
required. Specifically, if you want to use assumed-shape arrays, then
you have to use an INTERFACE!
And finally, Here is that same program with all
of the sub-programs contained in a MODULE, just like our random number
generator, except that this MODULE lives in the same file as the program
which USEs it.