diff --git a/matlab/BeijingPose2Graph.m b/matlab/BeijingPose2Graph.m
new file mode 100644
index 000000000..3f093494e
--- /dev/null
+++ b/matlab/BeijingPose2Graph.m
@@ -0,0 +1,60 @@
+load beijing.mat;
+load beijing_angles.mat;
+load beijing_graph.mat;
+%load beijing_config.mat;
+
+cov = [ 0.25, 0, 0; 0, 0.25, 0; 0, 0, 0.01];
+
+
+factors = Pose2Graph;
+factors2 = Pose2Graph;
+ord = Ordering();
+
+for i=1:50%length(ways)
+    if mod(i,50) == 0
+        fprintf(1, 'processing way %d\n', i);
+    end
+    for j=1:length(ways{i})-1
+        id1 = ways{i}(j);
+        id2 = ways{i}(j+1);
+        key1 = sprintf('x%d', id1);
+        key2 = sprintf('x%d', id2);
+
+        delta_x = points(id1,:) - points(id2,:);
+        delta_angle = angles(id1) - angles(id2);
+        measured = Pose2(delta_x(1), delta_x(2), delta_angle);
+
+        if pred(id1) == id2 || pred(id2) == id1 %% in the spanning tree
+            factor=Pose2Factor(key1,key2,measured, cov);
+            factors.push_back(factor);
+            ord.push_back(key1);
+            ord.push_back(key2);
+        else %% not in the spanning tree
+            factors2.push_back(Pose2Factor(key1,key2,measured, cov));
+        end
+    end
+end
+
+ord.unique();
+config=Pose2Config();
+n=size(points,1);
+for j=1:n
+    pose=Pose2(points(j,1),points(j,2),angles(j));
+     key = sprintf('x%d', j);
+    config.insert(key,pose);
+       if mod(j,50) == 0
+           key
+    end
+end
+save('beijing_config.mat','config');
+
+ord2=factors.getOrdering_();
+fprintf(1,'linearize...');
+LFG=factors.linearize_(config);
+fprintf(1,'linearize done!');
+ijs=LFG.sparse(ord2);
+fprintf(1,'sparse done');
+A=sparse(ijs(1,:),ijs(2,:),ijs(3,:)); 
+spy(A);
+%save('beijing_factors.mat', 'factors');
+
diff --git a/matlab/beijing.mat b/matlab/beijing.mat
new file mode 100644
index 000000000..6c2ae4b76
Binary files /dev/null and b/matlab/beijing.mat differ
diff --git a/matlab/beijing_angles.mat b/matlab/beijing_angles.mat
new file mode 100644
index 000000000..8128a433d
Binary files /dev/null and b/matlab/beijing_angles.mat differ
diff --git a/matlab/beijing_config.mat b/matlab/beijing_config.mat
new file mode 100644
index 000000000..4fb99d905
Binary files /dev/null and b/matlab/beijing_config.mat differ
diff --git a/matlab/beijing_graph.mat b/matlab/beijing_graph.mat
new file mode 100644
index 000000000..226272a0a
Binary files /dev/null and b/matlab/beijing_graph.mat differ
diff --git a/matlab/testPose2Factor.m b/matlab/testPose2Factor.m
index b549cd369..b0c45b24b 100644
--- a/matlab/testPose2Factor.m
+++ b/matlab/testPose2Factor.m
@@ -1,3 +1,4 @@
+% test linearize Pose2Factor
 cov = [ 0.25, 0, 0; 0, 0.25, 0; 0, 0, 0.01];
 key1='x1';
 key2='x2';
@@ -6,8 +7,20 @@ measured.print('Pose');
 factor=Pose2Factor(key1,key2,measured, cov);
 factor.print('Factor');
 
+
+% test linearize Pose2Graph 
+
 p1=Pose2(1.1,2,pi/2); % robot at (1.1,2) looking towards y (ground truth is at 1,2, see testPose2)
 p2=Pose2(-1,4.1,pi);  % robot at (-1,4) looking at negative (ground truth is at 4.1,2)
 config= Pose2Config() ;
+config.insert('x1',p1);
+config.insert('x2',p2);
 
-%fg = Pose2Graph;
\ No newline at end of file
+lf = factor.linearize(config);
+lf.print('lf ');
+
+factors = Pose2Graph;
+factors.push_back(factor);
+
+lfg=factors.linearize_(config);
+lfg.print('lfg');
\ No newline at end of file