|
MATHEMATICA
|
k[n_] = -(-1 + 2^(-n))^(-n)* (-2 + 2^(-n))^n *(-1 + 2^n) j[x_, n_] = (x^n - 2)^n/(k[n]*x^n*(x^n - 1)^(n - 1)) (* Farey like function*) f[x_] := 1/(j[x, 2]) /; 0 <= x <= 1/2 f[x_] := j[x, 2] /; 1/2 < x <= 2 ff[x_] = f[Mod[Abs[x], 2]] Plot[f[Mod[Abs[x], 2]], {x, 0, 2}] (*n=2 level*) b = 49*ReplacePart[Table[Coefficient[Series[ -3* (-2 + x^2)^2/(49* x^2 (-1 + x^2)), {x, 0, 30}], x^n], {n, -2, 30}], 3/49, 3] (* removing the zeros*) c = Flatten[Table[If[b[[n]] > 0, b[[n]], {}], {n, 1, Length[b]}]]
|