(PARI) x='x+O('x^50); v=Vec(serlaplace(atanh(sin(x)*sin(x)))); concat([0], vector(#v\2, n, v[2*n-1])) \\ G. C. Greubel, Oct 25 2018
(Magma) m:=30; R<x>:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( Argtanh(Sin(x)^2) )); [0] cat [Factorial(2*n+2)*b[2*n+1]: n in [0..Floor((m-4)/2)]]; // G. C. Greubel, Oct 25 2018