C----------------------------------------------------------------------- C IMSL Name: LSLTQ/DLSLTQ (Single/Double precision version) C C Purpose: Solve a complex tridiagonal system of linear equations. C C Usage: CALL LSLTQ (N, C, D, E, B) C C Example 1: C Declaration of variables INTEGER N PARAMETER (N=4) C COMPLEX B(N), C(N), D(N), E(N) CHARACTER CLABEL(1)*6, FMT*8, RLABEL(1)*4 EXTERNAL LSLTQ, WRCRL C DATA FMT/'(E13.6)'/ DATA CLABEL/'NUMBER'/ DATA RLABEL/'NONE'/ C C(*), D(*), E(*) and B(*) C contain the subdiagonal, C diagonal, superdiagonal and C right hand side. DATA C/(0.0,0.0), (-9.0,3.0), (2.0,7.0), (7.0,-4.0)/ DATA D/(3.0,-5.0), (4.0,-9.0), (-5.0,-7.0), (-2.0,-3.0)/ DATA E/(-9.0,8.0), (1.0,8.0), (8.0,3.0), (0.0,0.0)/ DATA B/(-16.0,-93.0), (128.0,179.0), (-60.0,-12.0), (9.0,-108.0)/ C C CALL LSLTQ (N, C, D, E, B) C Output the solution. CALL WRCRL ('Solution:', 1, N, B, 1, 0, FMT, RLABEL, CLABEL) END C Solution: C 1 2 C (-0.400000E+01,-0.700000E+01) (-0.700000E+01, 0.400000E+01) C C 3 4 C ( 0.700000E+01,-0.700000E+01) ( 0.900000E+01, 0.200000E+01)