Find the partial derivative of f(x, y) = x³ + 2x³y³ – 3y² + x + y, at (2,1)