diff --git a/wiki/pages/model.html b/wiki/pages/model.html
index cd9cb92dee34605bab893576ab8129a47faf0875..b03d8457e20da757b78a66aa57de4e51f0f21642 100644
--- a/wiki/pages/model.html
+++ b/wiki/pages/model.html
@@ -304,7 +304,7 @@
           muskone in the intestine can be obtained in combination with experimental determination. Because
           there is no animal experimental support, the data are manually drafted, and the calculation method
           is more meaningful than the calculation results.</p>
-        <button onclick="toggleCodeSnippet()">Expand the code</button>
+        <button id="Button1" onclick="toggleCodeSnippet()">Expand the code</button>
 
         <div id="codeSnippet" class="code-snippet">
           ```metlab
@@ -318,8 +318,85 @@
           k_dist = 0.001;
           k_excrete = 0.05;
           k_move = 0.02;
+
+          % Define the time range
+          tspan = [0 300]; % From 0 to 5 minutes
+          initial_conditions = [0 0 0 0 0]; % The initial condition is 0
+
+          % solve ODE
+          [t, y] = ode45(@(t,y) odefun(t, y, Q_inhale, k_exhale, k_perm, k_adh, k_diffMC, k_diffLC, k_dist,
+          k_excrete, k_move), tspan, initial_conditions);
+
+          % calculate V_inhale
+          V_inhale = Q_inhale / 5 * (heaviside(t) - heaviside(t-5));
+
+          figure('Position', [100, 100, 1200, 1000]);
+
+          % V_inhale(t)
+          subplot(3,2,1)
+          plot(t, V_inhale)
+          title('V_{inhale}(t)')
+          xlabel('Time (s)')
+          ylabel('V_{inhale}')
+
+          % Q_A(t)
+          subplot(3,2,2)
+          plot(t, y(:,1))
+          title('Q_A(t)')
+          xlabel('Time (s)')
+          ylabel('Q_A')
+
+          % Q_L(t)
+          subplot(3,2,3)
+          plot(t, y(:,2))
+          title('Q_L(t)')
+          xlabel('Time (s)')
+          ylabel('Q_L')
+
+          % Q_M(t)
+          subplot(3,2,4)
+          plot(t, y(:,3))
+          title('Q_M(t)')
+          xlabel('Time (s)')
+          ylabel('Q_M')
+
+          % Q_C(t)
+          subplot(3,2,5)
+          plot(t, y(:,4))
+          title('Q_C(t)')
+          xlabel('Time (s)')
+          ylabel('Q_C')
+
+          % Q_I(t)
+          subplot(3,2,6)
+          plot(t, y(:,5))
+          title('Q_I(t)')
+          xlabel('Time (s)')
+          ylabel('Q_I')
+
+          sgtitle('Simulation Results')
+
+          % ODE
+          function dydt = odefun(t, y, Q_inhale, k_exhale, k_perm, k_adh, k_diffMC, k_diffLC, k_dist,
+          k_excrete, k_move)
+          V_inhale = Q_inhale / 5 * (heaviside(t) - heaviside(t-5));
+          dydt = zeros(5,1);
+          dydt(1) = V_inhale - (k_exhale + k_perm) * y(1); % dQ_A/dt
           ```
         </div>
+        <script>
+          function toggleCodeSnippet() {
+                        var codeSnippet = document.getElementById("codeSnippet");
+                        var button = document.getElementById("Button1"); // 注意变量名通常使用小写开头
+                        if (codeSnippet.style.display === "none") {
+                            codeSnippet.style.display = "block";
+                            button.textContent = "Collapse the code"; // 使用之前选中的按钮元素
+                        } else {
+                            codeSnippet.style.display = "none";
+                            button.textContent = "Expand the code"; // 使用之前选中的按钮元素
+                        }
+                    }
+        </script>
     </div>
   </div>