I would let:
#x = 4sintheta#
#dx = 4costhetad theta#
#sqrt(16 - x^2) = 4costheta#
#x^3 = 64sin^3theta#
#=>1024int sin^3thetacos^2thetad theta#
#=1024int sin^2thetasinthetacos^2thetad theta#
#=1024int (1-cos^2theta)cos^2thetasinthetad theta#
At this point you can let:
#u = costheta#
#du = -sinthetad theta#
#=> 1024int (cos^2theta - 1)cos^2theta(-sintheta)d theta#
#= 1024int (u^2 - 1)u^2du#
#= 1024int u^4 - u^2du#
#= 1024(u^5/5 - u^3/3)#
#= 1024(cos^5theta/5 - cos^3theta/3)#
Since #costheta = sqrt(16-x^2)/4#:
#= 1024((sqrt(16-x^2)/4)^5/5 - (sqrt(16-x^2)/4)^3/3)#
#= 1024(((16-x^2)^(5/2)/1024)/5 - ((16-x^2)^(3/2)/64)/3)#
#= 1024((16-x^2)^(5/2)/(5*1024) - (16-x^2)^(3/2)/(3*64))#
#= cancel(1024)((16-x^2)^(5/2)/(5*cancel(1024)) - (16-x^2)^(3/2)/(3/16*cancel(1024)))#
#= color(green)(1/5(16-x^2)^("5/2") - 16/3(16-x^2)^("3/2") + C)#
Normally here would be okay, but we can go a bit further.
#= (16-x^2)^(3/2)(1/5(16-x^2) - 16/3)#
#= (16-x^2)^(3/2)(3/15(16-x^2) - 80/15)#
#= 1/15(16-x^2)^(3/2)(3(16-x^2) - 80)#
#= 1/15(16-x^2)^(3/2)(48-3x^2 - 80)#
#= 1/15(16-x^2)^(3/2)(-32-3x^2)#
#= color(blue)(-1/15(16-x^2)^("3/2")(3x^2 + 32) + C)#