Tính tích phân \(I = \int_0^2 {{x^2}\sqrt {{x^3} + 1} }dx \)

A \(I =  - {{16} \over 9}\)

B \(I =  {{52} \over 9}\)

C \(I =  {{16} \over 9}\)

D \(I =  - {{52} \over 9}\)