diff --git a/samples/matlab_experimental/flame.m b/samples/matlab_experimental/flame.m
index e80c50caf..a7f5b2363 100644
--- a/samples/matlab_experimental/flame.m
+++ b/samples/matlab_experimental/flame.m
@@ -48,7 +48,8 @@ function f = flame(gas, left, flow, right)
     f.setProfile(2, {'u', 'V'}, [0.0            1.0
                                  mdot0/rho0     -mdot1/rho0
                                  0.0            0.0]);
-    f.setProfile(2, 'T', [0.0 z1 1.0; t0 2000.0 t1]);
+    f.setProfile(2, 'T', [0.0, 1.0
+                          t0, t1]);
 
     for n = 1:gas.nSpecies
       nm = gas.speciesName(n);
@@ -63,11 +64,11 @@ function f = flame(gas, left, flow, right)
       else
         y1 = yeq(n);
       end
-      f.setProfile(2, nm, [0, z1, 1
-                           left.massFraction(n), yint, y1]);
+      f.setProfile(2, nm, [0, 1.0
+                           left.massFraction(n), y1]);
     end
 
     % set minimal grid refinement criteria
     f.setRefineCriteria(2, 10.0, 0.8, 0.8);
-
 end
+
diff --git a/samples/matlab_experimental/periodic_cstr.m b/samples/matlab_experimental/periodic_cstr.m
index 1d7e90ccb..363eb2552 100644
--- a/samples/matlab_experimental/periodic_cstr.m
+++ b/samples/matlab_experimental/periodic_cstr.m
@@ -23,7 +23,7 @@ function periodic_cstr
     help periodic_cstr
 
     % create the gas mixture
-    gas = Solution('h2o2.yaml');
+    gas = Solution('h2o2.yaml','ohmech');
 
     % pressure = 60 Torr, T = 770 K
     p = 60.0 * 133.3;
diff --git a/samples/matlab_experimental/surfreactor.m b/samples/matlab_experimental/surfreactor.m
index a94675016..699982e41 100644
--- a/samples/matlab_experimental/surfreactor.m
+++ b/samples/matlab_experimental/surfreactor.m
@@ -73,7 +73,7 @@ for n = 1:nSteps
   tim(n) = t;
   temp(n) = r.T;
   pres(n) = r.P - p0;
-  cov(n,:) = surf.coverages';
+  cov(n,:) = surf.X';
   x(n,:) = gas.moleFraction(names);
 end
 disp(['CPU time = ' num2str(cputime - t0)]);