--- a/test/laplace.m	Thu Jul 18 14:31:40 2013 +0200
+++ b/test/laplace.m	Thu Jul 18 14:46:49 2013 +0200
@@ -16,12 +16,16 @@
 #Set BC equal to zero on sides with label 2 and 4
 bc = fem_bc (V, @(x, y) 0.0, [2,4]);
 
-##Set the rhs
-#g = fem_rhs (V, @(x, y));
+##Set the coefficient for the rhs
+fv = @(x,y) 10*exp(-((x - 0.5).^2 + (y - 0.5).^2) / 0.02);
+f = fem_coeff ('f', fv);
+
+gv = @(x,y) sin (5.0 * x);
+g = fem_coeff ('g', gv);
 
 # #Get the matrices
-# A  = BilinearForm (V, V, bc, "Laplace");
-# L  = LinearForm (V, bc, "Laplace");
+# A  = BilinearForm (V, V, bc);
+# L  = LinearForm (V, bc, f, g);
 
 # #Solve it
 # u = A\L;