diff --git a/BayesValidRox/BayesInference/BayesInference.py b/BayesValidRox/BayesInference/BayesInference.py
index 8c52958d320c8929efc684ddf26f0769597ddd49..2e7be04bd493a184b66181305132e652910e4007 100644
--- a/BayesValidRox/BayesInference/BayesInference.py
+++ b/BayesValidRox/BayesInference/BayesInference.py
@@ -711,24 +711,34 @@ class BayesInference:
         figPosterior = corner.corner(self.Posterior_df.to_numpy(), labels=parNames,
                                      quantiles=[0.15, 0.5, 0.85],show_titles=True,
                                      title_fmt=self.Corner_title_fmt,
-                                     title_kwargs={"fontsize": 12})
+                                     title_kwargs={"fontsize": 20})
         
-        # # Extract the axes
-        # axes = np.array(figPosterior.axes).reshape((NofPa, NofPa))
+        # Extract the axes
+        axes = np.array(figPosterior.axes).reshape((NofPa, NofPa))
         
         # # Loop over the diagonal
         # for i in range(NofPa):
         #     ax = axes[i, i]
         #     ax.axvline(MAP_theta[0,i], ls='--', color="r")
         
-        # # Loop over the histograms
+        for i in range(NofPa):
+            ax_x = axes[-1, i]
+            ax_y = axes[i, 0]
+            ax_x.set_xlabel(parNames[i],rotation=45)
+            if i != 0:
+                ax_y.set_ylabel(parNames[i],rotation=45)
+            
+        # Loop over the histograms
         # for yi in range(NofPa):
         #     for xi in range(yi):
         #         ax = axes[yi, xi]
+        #         ax.set_xlabel(parNames[xi],rotation=45)
+        #         ax.set_ylabel(parNames[yi],rotation=45)
         #         ax.axvline(MAP_theta[0,xi], ls='--', color="r")
         #         ax.axhline(MAP_theta[0,yi], ls='--', color="r")
         #         ax.plot(MAP_theta[0,xi], MAP_theta[0,yi], "sr")
-        
+
+        # plt.yticks(rotation=45, horizontalalignment='right')
         figPosterior.set_size_inches((24,16)) 
         
         if self.emulator: