Calculate integral by CAS
Using default functions for sanity check
Created with the Wolfram Language