Numericalintegrate(x (( 2*Si(x) - 1 )^2 + ( 2*Ci(x) - 1 )^2 )^2 , x = 0 .. infinity )