Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- Clear[t, m, s, M, r, b, M, v, d, x, k, q, z, nn];
- (*b=N[Sum[Sum[(BernoulliB[2*r]/((2*r)!))*(-d k^(1-2 r) E^(I x) \
- Abs[StirlingS1[2*r-1,m]] Gamma[1+m,s Log[d*k]] \
- Log[d*k]^(-1-m)),{m,1,2*r-1}],{r,1,q-1}],30];*)
- b = 4;
- M = 7;
- v = 2;
- (*d=11;*)
- k = 8;
- q = 4;
- c = 1;
- nn = 10;
- (*r=19;
- m=5;*)
- h[t_] = Sum[(I d^(1/2 - I t) E^(I t x) n^(-(1/2) - I t))/(
- 2 (-x + Log[d] + Log[n])) - (
- I d^(1/2 + I t) E^(I t x) n^(-(1/2) + I t))/(
- 2 (x + Log[d] + Log[n])), {n, 1, k}] + -(1/2) I E^(x/2)
- ExpIntegralEi[1/2 I (I + 2 t) (x - Log[d k])] +
- 1/2 I E^(-x/2) ExpIntegralEi[1/2 (1 + 2 I t) (x + Log[d k])] - (
- I d^(1/2 - I t) E^(I t x) k^(-(1/2) - I t))/(
- 4 (-x + Log[d] + Log[k])) + (
- I d^(1/2 + I t) E^(I t N[x]) k^(-(1/2) + I t))/(
- 4 (x + Log[d] + Log[k]));
- t1 = 1/200;
- t2 = 100;
- g1 = ListLinePlot[
- Table[-Im[
- N[Sum[Sum[-I*(h[t2] - h[t1])*MoebiusMu[d], {d, Divisors[z]}]/
- z^c, {z, 1, nn}] -
- I*(N[Sum[
- Sum[(Sum[
- Sum[-(d/v) k^(1 - 2 r)
- E^(-(x/
- v) (1 -
- 2 (1/2 + I*t2))) ((BernoulliB[2*r]/((2*r)!))*
- Abs[StirlingS1[2*r - 1, m]]) Gamma[
- 1 + m, (1/2 + I*t2) Log[d*k]] Log[d*k]^(-1 -
- m) + (d*x)/v k^(
- 1 - 2 r) ((BernoulliB[2*r]/((2*r)!))*
- Abs[StirlingS1[2*r - 1, m]]) (1/
- x E^(-(x/v) + x (1/2 + I*t2))
- Gamma[1 + m, (1/2 + I*t2) Log[d*k]] +
- 1/x (Sum[
- Sum[(-1)^(m - kk)*
- m!/kk!/(N[x] - Log[d*k])^(m + 1)*
- E^(((1/2 + I*t2) - 1/v)*x)*(1/2 + I*t2)^(kk +
- p)*(d*k)^(-(1/2 + I*t2))*x^kk*Log[d*k]^p/p!, {kk,
- 0, m - p}], {p, 0, m}]) Log[d*k]^(1 + m)) Log[
- d*k]^(-1 - m), {m, 1, 2*r - 1}], {r, 1,
- q - 1}] -
- Sum[Sum[-(d/v) k^(1 - 2 r)
- E^(-(x/
- v) (1 -
- 2 (1/2 + I*t1))) ((BernoulliB[2*r]/((2*r)!))*
- Abs[StirlingS1[2*r - 1, m]]) Gamma[
- 1 + m, (1/2 + I*t1) Log[d*k]] Log[d*k]^(-1 -
- m) + (d*x)/v k^(
- 1 - 2 r) ((BernoulliB[2*r]/((2*r)!))*
- Abs[StirlingS1[2*r - 1, m]]) (1/
- x E^(-(x/v) + x (1/2 + I*t1))
- Gamma[1 + m, (1/2 + I*t1) Log[d*k]] +
- 1/x (Sum[
- Sum[(-1)^(m - kk)*m!/kk!/(x - Log[d*k])^(m + 1)*
- E^(((1/2 + I*t1) - 1/v)*x)*(1/2 + I*t1)^(kk +
- p)*(d*k)^(-(1/2 + I*t1))*x^kk*Log[d*k]^p/p!, {kk,
- 0, m - p}], {p, 0, m}]) Log[d*k]^(1 + m)) Log[
- d*k]^(-1 - m), {m, 1, 2*r - 1}], {r, 1, q - 1}])*
- MoebiusMu[d], {d, Divisors[z]}]/z^c, {z, 1, nn}],
- b]) + -I*(N[
- Sum[Sum[(Sum[
- Sum[-(-(d/v) k^(1 - 2 r)
- E^(-(x/
- v) (1 -
- 2 (1/2 - I*t2))) ((BernoulliB[2*r]/((2*r)!))*
- Abs[StirlingS1[2*r - 1, m]]) Gamma[
- 1 + m, (1/2 - I*t2) Log[d*k]] Log[d*k]^(-1 -
- m) + (d*x)/v k^(
- 1 - 2 r) ((BernoulliB[2*r]/((2*r)!))*
- Abs[StirlingS1[2*r - 1, m]]) (1/
- x E^(-(x/v) + x (1/2 - I*t2))
- Gamma[1 + m, (1/2 - I*t2) Log[d*k]] +
- 1/x (Sum[
- Sum[-(+1)^(m - kk)*
- m!/kk!/(N[x] + Log[d*k])^(m + 1)*
- E^((-(1/2 - I*t2) + 1/v)*x)*(1/2 - I*t2)^(kk +
- p)*(d*k)^(-(1/2 - I*t2))*x^kk*Log[d*k]^p/p!, {kk,
- 0, m - p}], {p, 0, m}]) Log[d*k]^(1 + m)) Log[
- d*k]^(-1 - m)), {m, 1, 2*r - 1}], {r, 1,
- q - 1}] -
- Sum[Sum[-(-(d/v) k^(1 - 2 r)
- E^(-(x/
- v) (1 -
- 2 (1/2 - I*t1))) ((BernoulliB[2*r]/((2*r)!))*
- Abs[StirlingS1[2*r - 1, m]]) Gamma[
- 1 + m, (1/2 - I*t1) Log[d*k]] Log[d*k]^(-1 -
- m) + (d*x)/v k^(
- 1 - 2 r) ((BernoulliB[2*r]/((2*r)!))*
- Abs[StirlingS1[2*r - 1, m]]) (1/
- x E^(-(x/v) + x (1/2 - I*t1))
- Gamma[1 + m, (1/2 - I*t1) Log[d*k]] +
- 1/x (Sum[
- Sum[-(+1)^(m - kk)*m!/kk!/(x + Log[d*k])^(m + 1)*
- E^((-(1/2 - I*t1) + 1/v)*x)*(1/2 - I*t1)^(kk +
- p)*(d*k)^(-(1/2 - I*t1))*x^kk*Log[d*k]^p/p!, {kk,
- 0, m - p}], {p, 0, m}]) Log[d*k]^(1 + m)) Log[
- d*k]^(-1 - m)), {m, 1, 2*r - 1}], {r, 1, q - 1}])*
- MoebiusMu[d], {d, Divisors[z]}]/z^c, {z, 1, nn}], b]),
- b]], {x, 1/200, 2, 1/200}], DataRange -> {1/200, 2},
- ImageSize -> Large, PlotStyle -> {Thickness[0.004]}];
- g2 = Table[Graphics[Line[{{Log[n], -45}, {Log[n], 0}}]], {n, 1, 7}];
- g3 = Table[Graphics[Text[Log[n], {Log[n], -50}]], {n, 1, 7}];
- Show[g1, g2, g3]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement