diff options
Diffstat (limited to 'src/ksq_sergio.x')
-rw-r--r-- | src/ksq_sergio.x | 49 |
1 files changed, 49 insertions, 0 deletions
diff --git a/src/ksq_sergio.x b/src/ksq_sergio.x new file mode 100644 index 0000000..8a6180b --- /dev/null +++ b/src/ksq_sergio.x @@ -0,0 +1,49 @@ + o1 = a**2 + o2 = -2.00000000000000d0*mass*rBL + o3 = rBL**2 + o4 = o1 + o2 + o3 + o5 = sqrt(o4) + o6 = cos(qgrd(j)) + o7 = o6**2 + o8 = o1*o7 + o9 = o3 + o8 + o10 = o9**2 + o11 = 1/o10 + o12 = 1.0!psi2dv(eta,q) + o13 = o12**2 + o14 = 1/o13 + o15 = sin(qgrd(j)) + o16 = o15**2 + o17 = -o1 + o18 = mass**2 + o19 = o17 + o18 + o20 = sqrt(o19) + o21 = 1/o3 + o22 = 2.50000000000000d-1*o19*o21 + o23 = 1.00000000000000d0 + o22 + o24 = o1 + o3 + o25 = 2.00000000000000d0*o24*o3 + o26 = o17 + o3 + o27 = o26*o9 + o28 = o25 + o27 + o29 = o23**2 + o30 = o10**2 + o31 = 1/o30 + o32 = o28**2 + o33 = o13**2 + o34 = o33**2 + o35 = o33*o34 + o36 = 1/o35 + o37 = o1**2 + o38 = o16**2 +c exc33(i,j) = 0 +c exc32(i,j) = -2.00000000000000d0*a*mass*o1*o11*o14*o15*o16*o5*o6 +c & *rBL +c exc31(i,j) = 5.0000000000000d-1*a*mass*o11*o14*o16*o20*o23*o28*e +c & xp(eta) +c exc22(i,j) = 0 +c exc21(i,j) = 0 +c exc11(i,j) = 0 + ksq(i,j) = 8.0000000000000d0*o1*o18*o3*o31*o36*o37*o38*o4*o7 + 5 + & .0000000000000d-1*o1*o16*o18*o19*o29*o31*o32*o36*exp(2.000000000 + & 00000d0*etagrd(i)) |