--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/laplace.m	Tue Jul 02 22:52:19 2013 +0200
@@ -0,0 +1,28 @@
+#Example of how we expect the Octave interface to look like
+#This file is our test for the TDD
+
+# Create mesh 
+init_dolfin ();
+
+x = y = linspace (0, 1, 32) ;
+msho = msh2m_structured_mesh (x, y, 1, [ 1 1 2 2 ])
+mshd = dolfin_mesh (msho);
+
+#Compile the problem defined in the .ufl file
+ffc ("Laplace.ufl");
+
+#Import the problem inside Octave
+V  = FunctionSpace (mshd, "Laplace");
+sd = SubDomain (mshd, @(x, y) abs (x - 1.0) < eps, "Laplace");
+bc = DirichletBC (V, @(x, y) 0.0, sd, "Laplace");
+
+#Get the matrices
+A  = BilinearForm (V, V, bc, "Laplace");
+L  = LinearForm (V, bc, "Laplace");
+
+#Solve it
+u = A\L;
+
+#Plot
+uf = Function (V, u);
+dolfin_plot (uf, "Solution")