(PARI) x='x+O('x^50); v=Vec(serlaplace(asinh(sin(x)^2))); 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!( Argsinh(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