def METHOD_NAME(mesh_quad, family, degree):
    x, y = SpatialCoordinate(mesh_quad)
    if family in ['CG', 'DG']:
        V = VectorFunctionSpace(mesh_quad, family, degree)
        f = Function(V).interpolate(as_vector([0.2 + y, 0.8*x + 0.2*y]))
    else:
        V = FunctionSpace(mesh_quad, family, degree)
        f = Function(V).project(as_vector([0.2 + y, 0.8*x + 0.2*y]))
    assert np.allclose([0.6, 0.56], f([0.6, 0.4]))
    assert np.allclose([1.1, 0.18], f([0.0, 0.9]))