div[x_, y_] = D[Cos[x + y^3], x] + D[Sin[y + x^3], y];
StreamDensityPlot[{{Cos[x + y^3], Sin[y + x^3]}, div[x, y]}, {x, -2, 2}, {y, 0, 2}, StreamPoints -> Fine, ColorFunction -> "MintColors", MaxRecursion -> 2, StreamStyle -> Black, ImageSize -> Large]