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