123456789101112131415161718192021222324252627282930 |
- *DECK PSGF
- FUNCTION PSGF (X, IZ, C, A, BH)
- C***BEGIN PROLOGUE PSGF
- C***SUBSIDIARY
- C***PURPOSE Subsidiary to BLKTRI
- C***LIBRARY SLATEC
- C***TYPE SINGLE PRECISION (PSGF-S)
- C***AUTHOR (UNKNOWN)
- C***SEE ALSO BLKTRI
- C***ROUTINES CALLED (NONE)
- C***REVISION HISTORY (YYMMDD)
- C 801001 DATE WRITTEN
- C 891214 Prologue converted to Version 4.0 format. (BAB)
- C 900402 Added TYPE section. (WRB)
- C***END PROLOGUE PSGF
- DIMENSION A(*) ,C(*) ,BH(*)
- C***FIRST EXECUTABLE STATEMENT PSGF
- FSG = 1.
- HSG = 1.
- DO 101 J=1,IZ
- DD = 1./(X-BH(J))
- FSG = FSG*A(J)*DD
- HSG = HSG*C(J)*DD
- 101 CONTINUE
- IF (MOD(IZ,2)) 103,102,103
- 102 PSGF = 1.-FSG-HSG
- RETURN
- 103 PSGF = 1.+FSG+HSG
- RETURN
- END
|