ceof1:={};

For[k=1,  k<10, k++, AppendTo[ceof1, Sin[N[Pi]/2*(-1+(2*I*Pi*k+N[Log[2]])/N[Log[3]])]*Gamma[(2*I*Pi*k+N[Log[2]])/N[Log[3]]-1]]];

Sin[N[Pi]/2*(-1+(2*I*Pi*100+N[Log[2]])/N[Log[3]])]*Gamma[(2*I*Pi*100+N[Log[2]])/N[Log[3]]-1]

d:=N[Log[2]]/N[Log[3]];

d

Clear[ceof1]

ceof1[[1]]

ceof1[1]

ceof2

For[k=1,  k<10, k++, AppendTo[ceof3, {3}]];

ceof3[[1]]

coef3[[1]]

For[i=1;t=x,i^2<10,i++,t=t^2+i;Print[t]]

For[i=0,i<4,i++,Print[i]]

ceof1

DiscretePlot[ceof1[k],{k,1,10}]

g[t_]:=Sin[N[Pi]/2*(-1+(2*I*Pi*t+N[Log[2]])/N[Log[3]])]*Gamma[(2*I*Pi*t+N[Log[2]])/N[Log[3]]-1];

g[100]

DiscretePlot[h[k],j[k],{k,1,1000}]

DiscretePlot[{h[k],j[k]},{k,1,1000}]

h[t_]:=Re[g[t]];

h[100]

j[t_]:=t^(d-3/2);

j[1]

coeftot[n_]:=g[n]epsi2[n,0][[2]];

coeftot[1]

g[2]

epsi2[2,0]

Abs[I]

Abs[1+I]

mod[n_]:=Abs[coeftot[n]]

mod[1]

DiscretePlot[{mod[10n], -0.017856881906666837*(10n)^{0.0010340834945649225}},{n,1,1000}]

cf:=Table[mod[10n],{n,1000}]

cf

FindFit[cf,d+b *x^(c),{d,b,c},x]

Plot[%,{x,0,10000}]

data:={}

For[h=1,h<1001, h++, AppendTo[data, {10h, mod[10h]}]];

data[[1]]

FindFit[data,x^(c),{c},x]

DiscretePlot[{mod[10n], (10n)^{-2.88991}},{n,1,1000}]

data2:={}

For[hh=1,hh<10001, hh++, AppendTo[data2,mod[hh]]]

data2[[1]]

FindFit[data2,x^(c),{b,c},x]

DiscretePlot[{mod[n],n^(-6)},{n,1000,5000}]