\[ \int \text {Si}(a+b x)^2 \, dx \]
Optimal antiderivative \[ \frac {2 \cos \left (b x +a \right ) \sinIntegral \left (b x +a \right )}{b}+\frac {\left (b x +a \right ) \sinIntegral \left (b x +a \right )^{2}}{b}-\frac {\sinIntegral \left (2 b x +2 a \right )}{b} \]
command
integrate(sin_integral(b*x+a)^2,x, algorithm="fricas")
Fricas 1.3.8 (sbcl 2.2.11.debian) via sagemath 9.6 output
\[ \frac {{\left (b x + a\right )} \operatorname {Si}\left (b x + a\right )^{2} + 2 \, \cos \left (b x + a\right ) \operatorname {Si}\left (b x + a\right ) - \operatorname {Si}\left (2 \, b x + 2 \, a\right )}{b} \]
Fricas 1.3.7 via sagemath 9.3 output
\[ {\rm integral}\left (\operatorname {Si}\left (b x + a\right )^{2}, x\right ) \]