dQ/dx - dP/dy = 3y+1

Двойной интеграл от этой функции не равен 0, так как внутри этого треугольника эта функция только положительна.