gcsgau1(3f) - [M_math] solve a system of simultaneous linear equations (LICENSE:PD)
Synopsis
Description
Pedigree
Written By
License
subroutine gcsgau1(n,a,b)
integer,parameter :: dp=kind(0.0d0) integer,intent(in) :: n real(kind=dp) :: a(11,11) real(kind=dp) :: b(*)
Solve a system of simultaneous linear equations of the form
** ** ** ** ** ** * A(1,1) A(1,2) ... A(1,N) * * X(1) * * B(1) * * * * * * * * A(2,1) A(2,2) ... A(2,N) * * X(2) * * B(2) * * * * * * * * . . . * * . * = * . * * . . . * * . * * . * * . . . * * . * * . * * * * * * * * A(N,1) A(N,2) ... A(N,N) * * X(N) * * B(N) * ** ** ** ** ** **where matrices A and B are known and matrix X is the set of unknowns to be determined. N is the number of equations.
Based on
Graphics Compatibility System 3-D Device-Dependent Subroutine GCSGAU1 Level 7
Fred Taylor, Computer Analysis Branch USAEWES, Vicksburg, MS. 39180
Public Domain
Nemo Release 3.1 | gcsgau1 (3) | February 23, 2025 |