from sympy import *
from spb import *
x, y, a, b = symbols("x, y, a, b")
v = [-sin(a * y), cos(b * x)]
plot_vector(
    v, (x, -3, 3), (y, -3, 3),
    params={a: (1, -2, 2), b: (1, -2, 2)},
    quiver_kw=dict(color="black", scale=30, headwidth=5),
    contour_kw={"cmap": "Blues_r", "levels": 15},
    grid=False, xlabel="x", ylabel="y", use_latex=False)