--- a/test/laplace.m	Thu Jul 11 22:44:38 2013 +0200
+++ b/test/laplace.m	Thu Jul 11 22:46:12 2013 +0200
@@ -3,18 +3,18 @@
 
 # Create mesh 
 fem_init_env ();
+pkg load msh;
 
-x = y = linspace (0, 1, 32) ;
-msho = msh2m_structured_mesh (x, y, 1, [ 1 1 2 2 ]);
+msho = msh2m_structured_mesh (0:0.05:1, 0:0.05:1, 1, 1:4);
 mshd = fem_init_mesh (msho);
 
 #Compile the problem defined in the .ufl file
-fem_ffc ("Laplace.ufl");
+#fem_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");
+V  = fem_fs (mshd);
+#Set BC equal to zero on sides with label 2 and 4
+#bc = fem_bc (V, @(x, y) 0.0, [2,4]);
 
 # #Get the matrices
 # A  = BilinearForm (V, V, bc, "Laplace");