--- a
+++ b/inst/example/Mixed-Poisson/MixedPoisson.ufl
@@ -0,0 +1,12 @@
+BDM = FiniteElement("BDM", triangle, 1)
+DG  = FiniteElement("DG", triangle, 0)
+W = BDM * DG
+
+(sigma, u) = TrialFunctions(W)
+(tau, v)   = TestFunctions(W)
+
+CG = FiniteElement("CG", triangle, 1)
+f = Coefficient(CG)
+
+a = (dot(sigma, tau) + div(tau)*u + div(sigma)*v)*dx
+L = - f*v*dx