Calculate the double integral \(\iint_R \frac{xy^2}{x^2 + 1} dA\), where \(R = \{(x, y) | 0 \le x \le 1, -3 \le y \le 3\}\).