Try x=4 tan theta. Then x = 4 sec^2theta d theta
So int x^3/sqrt{x^2+16} dx = int {64 tan^3 theta times 4 sec^2theta}/{4sectheta}d theta = 64 int tan^3theta sec theta d theta = 64 int tan^2theta tantheta sec theta d theta =64 int (sec^2theta -1) d(sec theta) = 64 ({sec^3theta}/3 -sec theta) = 64/3 sec theta (sec^2theta-3)=64/3 sqrt{1+x^2/16}(1+x^2/16-3) = 1/3 sqrt{x^2 +16}(x^2-32)