diff --git a/matlab/main_car.m b/matlab/main_car.m
index ac57c2b2d11b78d1b8c619df37e37b1dfe1df0cb..f6bf76f1e3eeb7edadb645b1b6d0bff1b605c83d 100644
--- a/matlab/main_car.m
+++ b/matlab/main_car.m
@@ -18,4 +18,4 @@ out_data=GPSaidedINS_car(in_data,settings);
 
 %% Plot the data 
 disp('Plot data')
-plot_data
\ No newline at end of file
+plot_single_trace
diff --git a/matlab/main_cycling.m b/matlab/main_cycling.m
index 1d21fcedbcec17fd9940bdd2559076602139d375..a0ba614425d7d87c2717e2aa0e3cc699b8318888 100644
--- a/matlab/main_cycling.m
+++ b/matlab/main_cycling.m
@@ -18,4 +18,4 @@ out_data=GPSaidedINS_cycling(in_data,settings);
 
 %% Plot the data 
 disp('Plot data')
-plot_data
+plot_single_trace
diff --git a/matlab/plot_data.m b/matlab/plot_data.m
deleted file mode 100644
index 03f5d00c90e7b3878291b9dcfa109f5ca442fa0b..0000000000000000000000000000000000000000
--- a/matlab/plot_data.m
+++ /dev/null
@@ -1,32 +0,0 @@
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-% Script for data plot:
-% plots the trajecotry and the stimation accuracy with the GPS turn on
-% signal
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
-% display the overall error (root mean square)
-out_data.error    
-
-% Display trajectory estimated with sensor fusion against trajectory for the
-% GPS alone always on
-h=zeros(1,2);
-figure(1)
-hold on
-plot(in_data.GNSS.pos_ned(2,:),in_data.GNSS.pos_ned(1,:),'b-');
-h(1)=plot(in_data.GNSS.pos_ned(2,:),in_data.GNSS.pos_ned(1,:),'b.');
-h(2)=plot(out_data.x_h(2,:),out_data.x_h(1,:),'-');
-title('Trajectory')
-ylabel('North [m]')
-xlabel('East [m]')
-legend(h,'GNSS position estimate','GNSS aided INS trajectory')
-axis equal
-grid on
-
-%plot control signal P and turn on signal of the GPS
-h=zeros(1,2);
-figure(2)
-hold on
-h(1)=plot(out_data.P_store);
-h(2)=plot(out_data.turn);
-legend(h,'sum of the trace of P','antenna turned ON')
-grid on
\ No newline at end of file
diff --git a/matlab/plot_tradeoff.m b/matlab/plot_tradeoff.m
index 95645e0bafe817b3a70b801db9fdbc65b4675982..b2560c33ab1b168be25598660ddfa1c8bf42d307 100644
--- a/matlab/plot_tradeoff.m
+++ b/matlab/plot_tradeoff.m
@@ -8,7 +8,7 @@ xlabel('energy')
 ylabel('error')
 title('cycling energy accuracy tracking trade-off')
 grid
-
+%iterate and plot the results for the different P thresholds
 for j=1:length(P_treshold)
     scatter(energy(:,j), error(:,j),20,color(j,:),'o','filled')
 end
\ No newline at end of file
diff --git a/matlab/run_many_cycling.m b/matlab/run_many_cycling.m
index d29e06041623c6f49e9c9a929681704819f8c843..d98bbb1854c19c329a9bf5a5ea0fe91163abda60 100644
--- a/matlab/run_many_cycling.m
+++ b/matlab/run_many_cycling.m
@@ -25,7 +25,7 @@ for j=1:length(P_treshold)
     disp(['Starting simulations with threshold for P equal to:', num2str(P_treshold(j))]) %display reached treshold
     settings.P_treshold=P_treshold(j);                   %set P reference
     for k=1:n_run                                        %iterate over the number of simulations per threshold
-        disp(['Simulation number:',num2str(k)])                                       %display simulation number
+        disp(['Simulation number:',num2str(k)])          %display simulation number
         out_data=GPSaidedINS_cycling(in_data,settings);  % Run the GNSS-aided INS
         energy(k,j)=out_data.energy;                     % store energy
         error(k,j)=out_data.error;                       % store positioning error