1 SUBROUTINE solsy(wm, iwm, x, tem)
8 INTEGER,
DIMENSION(*) :: iwm
9 REAL(rprec),
DIMENSION(*) :: wm, x, tem
14 COMMON /ls0001/ rowns(209), ccmax, el0, h, hmin, hmxi, hu, rc, tn
15 1 , uround, iownd(14), iowns(6), icf, ierpj, iersl, jcur, jstart
16 2 , kflag, l, meth, miter, maxord, maxcor, msbp, mxncf, n, nq,
18 INTEGER iownd, iowns, icf, ierpj, iersl, jcur, jstart, kflag, l
19 1 , meth, miter, maxord, maxcor, msbp, mxncf, n, nq, nst, nfe,
21 REAL(rprec) ::rowns, ccmax, el0, h, hmin, hmxi, hu, rc, tn,
26 INTEGER :: i, meband, ml, mu
27 REAL(rprec) :: di, hl0, phl0, r
57 CALL li_gesl (wm(3:3+n*n), n, n, iwm(21:21+n), x, 0)
67 di = 1.0_dp - r*(1.0_dp - 1.0_dp/wm(i+2))
68 IF (abs(di) == 0.0_dp)
go to 390
72 x(:n) = wm(3:n+2)*x(:n)
81 meband = 2*ml + mu + 1
82 CALL li_gbsl (wm(3:3+meband*n), meband, n, ml, mu,