--- a/example/Evolution.m	Wed Aug 07 18:25:27 2013 +0200
+++ b/example/Evolution.m	Wed Aug 07 18:26:26 2013 +0200
@@ -21,7 +21,10 @@
 
 u0 = Expression ('u0', @(x,y) 10*exp(-((x - 0.5)^2 + (y - 0.5)^2) / 0.02));
 
-A = Evolution_BilinearForm (V, bc, k);
+a = Evolution_BilinearForm (V, k);
+b = Evolution_LinearForm (V, k, u0);
+
+A = assemble (a, bc);
 
 # solve the problem for each time step
 # We need to update only the lhs
@@ -30,7 +33,7 @@
 
   # we can pass u0 to the lhs indifferently as a fem_coeff or
   # as a fem_func
-  b = Evolution_LinearForm (V, k, u0, bc);
+  b = assemble (L, bc);
 
   u = A \ b;