Subroutine QSeries(s,L,q,qp,q2p) c This subroutine computes the value of q using the power series c solutions for the case omega = 0. Implicit None Integer i,L Real*8 q,qp,q2p,a3,a4,a5,a6,a7,a8 c Note a8 and b8 are not used in this program but are in PSeries Real*8 b2,b3,b4,b5,b6,b7,b8 Real*8 q0,q1,q2,q3,q4,q5,q6,q7 Real*8 alpha,ll1,xl,xi,s Real*8 omi,delom,nom COMMON/PAR/xi,omi,delom,nom Common/coeff/a3,a4,a5,a6,a7,a8,b2,b3,b4,b5,b6,b7,b8 xl = L ll1 = xl*(xl+1.d0) alpha = (-b2 - dsqrt(b2**2 + 4*b2*ll1 + 8*b2*xi - 8*b2**2*xi))/(2.*b2) q0 = 1.d0 q1 = (-4*alpha*b2*q0 - a3*alpha*b2*q0 - 3*alpha*b3*q0 - - 2*alpha**2*b3*q0 - 4*ll1*q0 - 8*q0*xi - 16*b2*q0*xi - - 6*a3*b2*q0*xi - 6*b3*q0*xi)/ - (2.*(2*b2 + 3*alpha*b2 + alpha**2*b2 - ll1 - 2*xi + 2*b2*xi)) q2 = (4*alpha*b2*q0 + a3**2*alpha*b2*q0 - 2*a4*alpha*b2*q0 - - 4*alpha*b3*q0 - a3*alpha*b3*q0 - 4*alpha*b4*q0 - - 2*alpha**2*b4*q0 + 6*ll1*q0 - 4*b2*q1 - a3*b2*q1 - 4*alpha*b2*q1 - - a3*alpha*b2*q1 - 5*b3*q1 - 7*alpha*b3*q1 - 2*alpha**2*b3*q1 - - 4*ll1*q1 + 12*q0*xi + 12*b2*q0*xi - 4*a3*b2*q0*xi + - 7*a3**2*b2*q0*xi - 16*a4*b2*q0*xi - 20*b3*q0*xi - 7*a3*b3*q0*xi - -8*b4*q0*xi - 8*q1*xi - 16*b2*q1*xi - 6*a3*b2*q1*xi - 6*b3*q1*xi)/ - (2.*(6*b2 + 5*alpha*b2 + alpha**2*b2 - ll1 - 2*xi + 2*b2*xi)) q3 = (-4*alpha*b2*q0 - a3**3*alpha*b2*q0 + 3*a3*a4*alpha*b2*q0 - - 3*a5*alpha*b2*q0 + 4*alpha*b3*q0 + a3**2*alpha*b3*q0 - - 2*a4*alpha*b3*q0 - 4*alpha*b4*q0 - a3*alpha*b4*q0 - - 5*alpha*b5*q0 - 2*alpha**2*b5*q0 - 8*ll1*q0 + 4*b2*q1 + - a3**2*b2*q1 - 2*a4*b2*q1 + 4*alpha*b2*q1 + a3**2*alpha*b2*q1 - - 2*a4*alpha*b2*q1 - 4*b3*q1 - a3*b3*q1 - 4*alpha*b3*q1 - - a3*alpha*b3*q1 - 6*b4*q1 - 8*alpha*b4*q1 - 2*alpha**2*b4*q1 + - 6*ll1*q1 - 8*b2*q2 - 2*a3*b2*q2 - 4*alpha*b2*q2 - a3*alpha*b2*q2 - - 14*b3*q2 - 11*alpha*b3*q2 - 2*alpha**2*b3*q2 - 4*ll1*q2 - - 16*q0*xi - 8*b2*q0*xi + 4*a3*b2*q0*xi + 4*a3**2*b2*q0*xi - - 8*a3**3*b2*q0*xi - 8*a4*b2*q0*xi + 26*a3*a4*b2*q0*xi - - 30*a5*b2*q0*xi + 16*b3*q0*xi - 4*a3*b3*q0*xi + 8*a3**2*b3*q0*xi - -18*a4*b3*q0*xi - 24*b4*q0*xi - 8*a3*b4*q0*xi - 10*b5*q0*xi + - 12*q1*xi + 12*b2*q1*xi - 4*a3*b2*q1*xi + 7*a3**2*b2*q1*xi - - 16*a4*b2*q1*xi - 20*b3*q1*xi - 7*a3*b3*q1*xi - 8*b4*q1*xi - - 8*q2*xi - 16*b2*q2*xi - 6*a3*b2*q2*xi - 6*b3*q2*xi)/ - (2.*(12*b2 + 7*alpha*b2 + alpha**2*b2 - ll1 - 2*xi + 2*b2*xi)) q4 = (4*alpha*b2*q0 + a3**4*alpha*b2*q0 - 4*a3**2*a4*alpha*b2*q0 + - 2*a4**2*alpha*b2*q0 + 4*a3*a5*alpha*b2*q0 - 4*a6*alpha*b2*q0 - - 4*alpha*b3*q0 - a3**3*alpha*b3*q0 + 3*a3*a4*alpha*b3*q0 - - 3*a5*alpha*b3*q0 + 4*alpha*b4*q0 + a3**2*alpha*b4*q0 - - 2*a4*alpha*b4*q0 - 4*alpha*b5*q0 - a3*alpha*b5*q0 - - 6*alpha*b6*q0 - 2*alpha**2*b6*q0 + 10*ll1*q0 - 4*b2*q1 - - a3**3*b2*q1 + 3*a3*a4*b2*q1 - 3*a5*b2*q1 - 4*alpha*b2*q1 - - a3**3*alpha*b2*q1 + 3*a3*a4*alpha*b2*q1 - 3*a5*alpha*b2*q1 + - 4*b3*q1 + a3**2*b3*q1 - 2*a4*b3*q1 + 4*alpha*b3*q1 + - a3**2*alpha*b3*q1 - 2*a4*alpha*b3*q1 - 4*b4*q1 - a3*b4*q1 - - 4*alpha*b4*q1 - a3*alpha*b4*q1 - 7*b5*q1 - 9*alpha*b5*q1 - - 2*alpha**2*b5*q1 - 8*ll1*q1 + 8*b2*q2 + 2*a3**2*b2*q2 - - 4*a4*b2*q2 + 4*alpha*b2*q2 + a3**2*alpha*b2*q2 - - 2*a4*alpha*b2*q2 - 8*b3*q2 - 2*a3*b3*q2 - 4*alpha*b3*q2 - - a3*alpha*b3*q2 - 16*b4*q2 - 12*alpha*b4*q2 - 2*alpha**2*b4*q2 + - 6*ll1*q2 - 12*b2*q3 - 3*a3*b2*q3 - 4*alpha*b2*q3 - - a3*alpha*b2*q3 - 27*b3*q3 - 15*alpha*b3*q3 - 2*alpha**2*b3*q3 - - 4*ll1*q3 + 20*q0*xi + 4*b2*q0*xi - 4*a3*b2*q0*xi - - 4*a3**2*b2*q0*xi - 4*a3**3*b2*q0*xi + 9*a3**4*b2*q0*xi + - 8*a4*b2*q0*xi + 12*a3*a4*b2*q0*xi - 38*a3**2*a4*b2*q0*xi + - 20*a4**2*b2*q0*xi - 12*a5*b2*q0*xi + 42*a3*a5*b2*q0*xi - - 48*a6*b2*q0*xi - 12*b3*q0*xi + 4*a3*b3*q0*xi + 4*a3**2*b3*q0*xi - - 9*a3**3*b3*q0*xi - 8*a4*b3*q0*xi + 29*a3*a4*b3*q0*xi - - 33*a5*b3*q0*xi + 20*b4*q0*xi - 4*a3*b4*q0*xi + 9*a3**2*b4*q0*xi - -20*a4*b4*q0*xi - 28*b5*q0*xi - 9*a3*b5*q0*xi - 12*b6*q0*xi - - 16*q1*xi - 8*b2*q1*xi + 4*a3*b2*q1*xi + 4*a3**2*b2*q1*xi - - 8*a3**3*b2*q1*xi - 8*a4*b2*q1*xi + 26*a3*a4*b2*q1*xi - - 30*a5*b2*q1*xi + 16*b3*q1*xi - 4*a3*b3*q1*xi + 8*a3**2*b3*q1*xi - -18*a4*b3*q1*xi - 24*b4*q1*xi - 8*a3*b4*q1*xi - 10*b5*q1*xi + - 12*q2*xi + 12*b2*q2*xi - 4*a3*b2*q2*xi + 7*a3**2*b2*q2*xi - - 16*a4*b2*q2*xi - 20*b3*q2*xi - 7*a3*b3*q2*xi - 8*b4*q2*xi - - 8*q3*xi - 16*b2*q3*xi - 6*a3*b2*q3*xi - 6*b3*q3*xi)/ - (2.*(20*b2 + 9*alpha*b2 + alpha**2*b2 - ll1 - 2*xi + 2*b2*xi)) q5 = (-4*alpha*b2*q0 - a3**5*alpha*b2*q0 + 5*a3**3*a4*alpha*b2*q0 - - 5*a3*a4**2*alpha*b2*q0 - 5*a3**2*a5*alpha*b2*q0 + - 5*a4*a5*alpha*b2*q0 + 5*a3*a6*alpha*b2*q0 - 5*a7*alpha*b2*q0 + - 4*alpha*b3*q0 + a3**4*alpha*b3*q0 - 4*a3**2*a4*alpha*b3*q0 + - 2*a4**2*alpha*b3*q0 + 4*a3*a5*alpha*b3*q0 - 4*a6*alpha*b3*q0 - - 4*alpha*b4*q0 - a3**3*alpha*b4*q0 + 3*a3*a4*alpha*b4*q0 - - 3*a5*alpha*b4*q0 + 4*alpha*b5*q0 + a3**2*alpha*b5*q0 - - 2*a4*alpha*b5*q0 - 4*alpha*b6*q0 - a3*alpha*b6*q0 - - 7*alpha*b7*q0 - 2*alpha**2*b7*q0 - 12*ll1*q0 + 4*b2*q1 + - a3**4*b2*q1 - 4*a3**2*a4*b2*q1 + 2*a4**2*b2*q1 + 4*a3*a5*b2*q1 - - 4*a6*b2*q1 + 4*alpha*b2*q1 + a3**4*alpha*b2*q1 - - 4*a3**2*a4*alpha*b2*q1 + 2*a4**2*alpha*b2*q1 + - 4*a3*a5*alpha*b2*q1 - 4*a6*alpha*b2*q1 - 4*b3*q1 - a3**3*b3*q1 + - 3*a3*a4*b3*q1 - 3*a5*b3*q1 - 4*alpha*b3*q1 - a3**3*alpha*b3*q1 + - 3*a3*a4*alpha*b3*q1 - 3*a5*alpha*b3*q1 + 4*b4*q1 + a3**2*b4*q1 - - 2*a4*b4*q1 + 4*alpha*b4*q1 + a3**2*alpha*b4*q1 - - 2*a4*alpha*b4*q1 - 4*b5*q1 - a3*b5*q1 - 4*alpha*b5*q1 - - a3*alpha*b5*q1 - 8*b6*q1 - 10*alpha*b6*q1 - 2*alpha**2*b6*q1 + - 10*ll1*q1 - 8*b2*q2 - 2*a3**3*b2*q2 + 6*a3*a4*b2*q2 - 6*a5*b2*q2 - - 4*alpha*b2*q2 - a3**3*alpha*b2*q2 + 3*a3*a4*alpha*b2*q2 - - 3*a5*alpha*b2*q2 + 8*b3*q2 + 2*a3**2*b3*q2 - 4*a4*b3*q2 + - 4*alpha*b3*q2 + a3**2*alpha*b3*q2 - 2*a4*alpha*b3*q2 - 8*b4*q2 - - 2*a3*b4*q2 - 4*alpha*b4*q2 - a3*alpha*b4*q2 - 18*b5*q2 - - 13*alpha*b5*q2 - 2*alpha**2*b5*q2 - 8*ll1*q2 + 12*b2*q3 + - 3*a3**2*b2*q3 - 6*a4*b2*q3 + 4*alpha*b2*q3 + a3**2*alpha*b2*q3 - - 2*a4*alpha*b2*q3 - 12*b3*q3 - 3*a3*b3*q3 - 4*alpha*b3*q3 - - a3*alpha*b3*q3 - 30*b4*q3 - 16*alpha*b4*q3 - 2*alpha**2*b4*q3 + - 6*ll1*q3 - 16*b2*q4 - 4*a3*b2*q4 - 4*alpha*b2*q4 - - a3*alpha*b2*q4 - 44*b3*q4 - 19*alpha*b3*q4 - 2*alpha**2*b3*q4 - - 4*ll1*q4 - 24*q0*xi + 4*a3*b2*q0*xi + 4*a3**2*b2*q0*xi + - 4*a3**3*b2*q0*xi + 4*a3**4*b2*q0*xi - 10*a3**5*b2*q0*xi - - 8*a4*b2*q0*xi - 12*a3*a4*b2*q0*xi - 16*a3**2*a4*b2*q0*xi + - 52*a3**3*a4*b2*q0*xi + 8*a4**2*b2*q0*xi - 54*a3*a4**2*b2*q0*xi + - 12*a5*b2*q0*xi + 16*a3*a5*b2*q0*xi - 56*a3**2*a5*b2*q0*xi + - 58*a4*a5*b2*q0*xi - 16*a6*b2*q0*xi + 62*a3*a6*b2*q0*xi - - 70*a7*b2*q0*xi + 8*b3*q0*xi - 4*a3*b3*q0*xi - 4*a3**2*b3*q0*xi - - 4*a3**3*b3*q0*xi + 10*a3**4*b3*q0*xi + 8*a4*b3*q0*xi + - 12*a3*a4*b3*q0*xi - 42*a3**2*a4*b3*q0*xi + 22*a4**2*b3*q0*xi - - 12*a5*b3*q0*xi + 46*a3*a5*b3*q0*xi - 52*a6*b3*q0*xi - - 16*b4*q0*xi + 4*a3*b4*q0*xi + 4*a3**2*b4*q0*xi - - 10*a3**3*b4*q0*xi - 8*a4*b4*q0*xi + 32*a3*a4*b4*q0*xi - - 36*a5*b4*q0*xi + 24*b5*q0*xi - 4*a3*b5*q0*xi + 10*a3**2*b5*q0*xi - -22*a4*b5*q0*xi - 32*b6*q0*xi - 10*a3*b6*q0*xi - 14*b7*q0*xi + - 20*q1*xi + 4*b2*q1*xi - 4*a3*b2*q1*xi - 4*a3**2*b2*q1*xi - - 4*a3**3*b2*q1*xi + 9*a3**4*b2*q1*xi + 8*a4*b2*q1*xi + - 12*a3*a4*b2*q1*xi - 38*a3**2*a4*b2*q1*xi + 20*a4**2*b2*q1*xi - - 12*a5*b2*q1*xi + 42*a3*a5*b2*q1*xi - 48*a6*b2*q1*xi - - 12*b3*q1*xi + 4*a3*b3*q1*xi + 4*a3**2*b3*q1*xi - - 9*a3**3*b3*q1*xi - 8*a4*b3*q1*xi + 29*a3*a4*b3*q1*xi - - 33*a5*b3*q1*xi + 20*b4*q1*xi - 4*a3*b4*q1*xi + 9*a3**2*b4*q1*xi - -20*a4*b4*q1*xi - 28*b5*q1*xi - 9*a3*b5*q1*xi - 12*b6*q1*xi - - 16*q2*xi - 8*b2*q2*xi + 4*a3*b2*q2*xi + 4*a3**2*b2*q2*xi - - 8*a3**3*b2*q2*xi - 8*a4*b2*q2*xi + 26*a3*a4*b2*q2*xi - - 30*a5*b2*q2*xi + 16*b3*q2*xi - 4*a3*b3*q2*xi + 8*a3**2*b3*q2*xi - -18*a4*b3*q2*xi - 24*b4*q2*xi - 8*a3*b4*q2*xi - 10*b5*q2*xi + - 12*q3*xi + 12*b2*q3*xi - 4*a3*b2*q3*xi + 7*a3**2*b2*q3*xi - - 16*a4*b2*q3*xi - 20*b3*q3*xi - 7*a3*b3*q3*xi - 8*b4*q3*xi - - 8*q4*xi - 16*b2*q4*xi - 6*a3*b2*q4*xi - 6*b3*q4*xi)/ - (2.*(30*b2 + 11*alpha*b2 + alpha**2*b2 - ll1 - 2*xi + 2*b2*xi)) q = s**alpha*(q0+q1*s+q2*s**2+q3*s**3+q4*s**4+q5*s**5) qp = alpha*s**(alpha-1.d0)*(q0+q1*s+q2*s**2+q3*s**3+q4*s**4+q5*s**5) - + s**alpha*(q1+2.*q2*s+3.*q3*s**2+4.*q4*s**3+5.*q5*s**4) q2p = alpha*(alpha-1.d0)*s**(alpha-2.d0)*(q0+q1*s+q2*s**2+q3*s**3+q4*s**4+q5*s**5) - + 2.*alpha*s**(alpha-1.d0)*(q1+2.*q2*s+3.*q3*s**2+4.*q4*s**3+5.*q5*s**4) - + s**alpha*(2.*q2+6.*q3*s+12.*q4*s**2+20.*q5*s**3) Return End