src/quadpkd/dqc25f.f

104      2  D1MACH,D2,ESTC,ESTS,F,FVAL,HLGTH,OFLOW,OMEGA,PARINT,PAR2,PAR22,
187       PAR2 = PARINT*PARINT
188       PAR22 = PAR2+0.2D+01
195       V(2) = (0.8D+01*COSPAR+(PAR2+PAR2-0.8D+01)*SINPAR/PARINT)/PAR2
196       V(3) = (0.32D+02*(PAR2-0.12D+02)*COSPAR+(0.2D+01*
197      1  ((PAR2-0.80D+02)*PAR2+0.192D+03)*SINPAR)/PARINT)/(PAR2*PAR2)
212         D2(K) = (AN-0.1D+01)*(AN-0.2D+01)*PAR2
213         D1(K+1) = (AN+0.3D+01)*(AN+0.4D+01)*PAR2
220       V(4) = V(4)-0.56D+02*PAR2*V(3)
222       ASAP = (((((0.210D+03*PAR2-0.1D+01)*COSPAR-(0.105D+03*PAR2
223      1  -0.63D+02)*ASS)/AN2-(0.1D+01-0.15D+02*PAR2)*COSPAR
225       V(NOEQU+3) = V(NOEQU+3)-0.2D+01*ASAP*PAR2*(AN-0.1D+01)*
244      1  +AS-PAR2*(AN+0.1D+01)*(AN+0.2D+01)*V(I-2))/
245      2  (PAR2*(AN-0.1D+01)*(AN-0.2D+01))
254       V(1) = 0.2D+01*(SINPAR-PARINT*COSPAR)/PAR2
255       V(2) = (0.18D+02-0.48D+02/PAR2)*SINPAR/PAR2
256      1  +(-0.2D+01+0.48D+02/PAR2)*COSPAR/PARINT
269         D2(K) = (AN-0.1D+01)*(AN-0.2D+01)*PAR2
270         D1(K+1) = (AN+0.3D+01)*(AN+0.4D+01)*PAR2
277       V(3) = V(3)-0.42D+02*PAR2*V(2)
279       ASAP = (((((0.105D+03*PAR2-0.63D+02)*ASS+(0.210D+03*PAR2
280      1  -0.1D+01)*SINPAR)/AN2+(0.15D+02*PAR2-0.1D+01)*SINPAR-
282       V(NOEQU+2) = V(NOEQU+2)-0.2D+01*ASAP*PAR2*(AN-0.1D+01)
300      1  +AC-PAR2*(AN+0.1D+01)*(AN+0.2D+01)*V(I-2))
301      2  /(PAR2*(AN-0.1D+01)*(AN-0.2D+01))