diff --git a/identification/index.org b/identification/index.org
index 620f1a7..545ebf3 100644
--- a/identification/index.org
+++ b/identification/index.org
@@ -385,6 +385,29 @@ initializeAxisc();
 #+end_src
 
 ** Bis
+#+begin_src matlab :results none
+  %% Initialize Ground
+  initializeGround();
+
+  %% Initialize Granite
+  initializeGranite(struct('rigid', false));
+
+  %% Initialize Translation stage
+  initializeTy(struct('rigid', false));
+
+  %% Initialize Tilt Stage
+  initializeRy(struct('rigid', false));
+
+  %% Initialize Spindle
+  initializeRz(struct('rigid', false));
+
+  %% Initialize Hexapod Symétrie
+  initializeMicroHexapod(struct('rigid', false));
+
+  %% Initialize Center of gravity compensation
+  initializeAxisc();
+#+end_src
+
 For Granite Fx to Tx/Ty/Tz/Rx/Ry/Rz
 #+begin_src matlab
   dirs = {'x', 'y', 'z', 'rx', 'ry', 'rz'};
diff --git a/simscape/index.org b/simscape/index.org
index e8b8651..66b112e 100644
--- a/simscape/index.org
+++ b/simscape/index.org
@@ -1069,7 +1069,7 @@ end
 This Matlab function is accessible [[file:../src/initializeGranite.m][here]].
 
 #+begin_src matlab
-  function [granite] = initializeGranite()
+  function [granite] = initializeGranite(opts_param)
       %% Default values for opts
       opts = struct('rigid', false);
 
@@ -1399,7 +1399,11 @@ This Matlab function is accessible [[file:../src/initializeMicroHexapod.m][here]
       Leg = struct();
 
       Leg.stroke     = 10e-3; % Maximum Stroke of each leg [m]
-      Leg.k.ax       = 5e7;   % Stiffness of each leg [N/m]
+      if opts.rigid
+        Leg.k.ax     = 1e12; % Stiffness of each leg [N/m]
+      else
+        Leg.k.ax     = 2e7; % Stiffness of each leg [N/m]
+      end
       Leg.ksi.ax     = 0.1;     % Modal damping ksi = 1/2*c/sqrt(km) []
       Leg.rad.bottom = 25;    % Radius of the cylinder of the bottom part [mm]
       Leg.rad.top    = 17;    % Radius of the cylinder of the top part [mm]
@@ -1585,8 +1589,8 @@ This Matlab function is accessible [[file:../src/initializeAxisc.m][here]].
       axisc.m      = 40; % TODO [kg]
 
       %% Axis Compensator - Dynamical Properties
-      axisc.k.ax   = 1; % TODO [N*m/deg)]
-      axisc.c.ax   = (1/1)*sqrt(axisc.k.ax/axisc.m);
+      % axisc.k.ax   = 1; % TODO [N*m/deg)]
+      % axisc.c.ax   = (1/1)*sqrt(axisc.k.ax/axisc.m);
 
       %% Save
       save('./mat/stages.mat', 'axisc', '-append');
diff --git a/simscape/sim_micro_station_modal_analysis_com.slx b/simscape/sim_micro_station_modal_analysis_com.slx
index 3a0bf32..5a0166b 100644
Binary files a/simscape/sim_micro_station_modal_analysis_com.slx and b/simscape/sim_micro_station_modal_analysis_com.slx differ