Index: /issm/trunk-jpl/examples/Inversion/CheatSheet.m
===================================================================
--- /issm/trunk-jpl/examples/Inversion/CheatSheet.m	(revision 20711)
+++ /issm/trunk-jpl/examples/Inversion/CheatSheet.m	(revision 20712)
@@ -1,5 +1,5 @@
 step=1;
 if step==1
-	%Generate observation
+	%Generate observations
 	md = model;
 	md = triangle(md,'DomainOutline.exp',100000);
@@ -23,9 +23,8 @@
 	md.friction.coefficient(find(md.mesh.x<600000 & md.mesh.x>400000))=10;
 
-	md = solve(md,StressbalanceSolutionEnum);
+	md = solve(md,StressbalanceSolutionEnum());
 
-	plotmodel(md,'data',md.friction.coefficient,'figure',1);
-	plotmodel(md,'data',md.results.StressbalanceSolution.Vel,'figure',2);
-
+	plotmodel(md,'axis#all','tight','data',md.friction.coefficient,'caxis',[0 100],'title','"True" \alpha',...
+		'data',md.results.StressbalanceSolution.Vel,'title','"observed velocities"')
 	save model1 md
 end
@@ -42,8 +41,8 @@
 	md.inversion.vel_obs = md.results.StressbalanceSolution.Vel;
 
-	md = solve(md,StressbalanceSolutionEnum);
+	md = solve(md,StressbalanceSolutionEnum());
 
-	plotmodel(md,'data',md.friction.coefficient,'figure',1);
-	plotmodel(md,'data',md.results.StressbalanceSolution.Vel,'figure',2);
+	plotmodel(md,'axis#all','tight','data',md.friction.coefficient,'caxis',[0 100],'title','\alpha first guess',...
+		'data',md.results.StressbalanceSolution.Vel,'title','modeled velocities')
 
 	save model2 md
@@ -62,6 +61,6 @@
 
 	%Hands on: STEP 2
-	md.inversion.cost_functions = 101;
-	md.inversion.cost_functions_coefficients = ones(md.mesh.numberofvertices,1);
+	md.inversion.cost_functions = [101 103];
+	md.inversion.cost_functions_coefficients = ones(md.mesh.numberofvertices,2);
 
 	%Hands on: STEP 3
@@ -71,7 +70,7 @@
 	%Go solve!
 	md.verbose=verbose(0);
-	md=solve(md,StressbalanceSolutionEnum);
+	md=solve(md,StressbalanceSolutionEnum());
 
-	plotmodel(md,'data',md.results.StressbalanceSolution.FrictionCoefficient,'figure',1,'caxis',[10 50]);
-	plotmodel(md,'data',md.results.StressbalanceSolution.Vel,'figure',2);  
+	plotmodel(md,'axis#all','tight','data',md.results.StressbalanceSolution.FrictionCoefficient,'caxis',[0 100],'title','inferred \alpha',...
+		'data',md.results.StressbalanceSolution.Vel,'title','modeled velocities')
 end