vector  = VectorElement("Lagrange", triangle, 1)
scalar  = FiniteElement("Discontinuous Lagrange", triangle, 0)
element = vector * scalar

(v, q) = TestFunctions(element)
(u, p) = TrialFunctions(element)

h = Coefficient(scalar)
delta = 0.2*h*h

a = (inner(grad(v), grad(u)) - div(v)*p + q*div(u) + delta*inner(grad(q), grad(p)))*dx
