import numpy as np

# Evaluate the integral

result $=(4+(1/3)*$np$.$sqrt$(4)**3)-(1+(1/3)*$np$.$sqrt$(1)**3)$

# Print the result to $8$ digits past the decimal point

print$($f$\text{'}$Pi is $\{$result:$\mathrm{.}8$f$\}\text{'})$

# Calculate the difference from numpy.pi and print to $15$ digits past the decimal point

difference $=$ result $-$ np$.$pi

print$($f$\text{'}$Difference from numpy.pi is: $\{$difference:$\mathrm{.}15$f$\}\text{'})$