dsh-math-calculator 2.3.1 → 2.3.3

package/CHANGELOG.md

@@ -2,6 +2,197 @@
 
 本项目遵循 [Keep a Changelog](https://keepachangelog.com/zh-CN/1.1.0/) 与 [语义化版本](https://semver.org/lang/zh-CN/)。
 
+## [2.3.3] — 2026-11-08
+
+针对三份实践文档（不足清单 17 条、能力评估 4 处缺口、结构性不足 8 条）
+逐条实测核对。**结论是多数已在本版之前修好**，本版只做确实还缺的四项，
+外加一处发现路径重构。
+
+### 新增（formula 批量调用）
+
+文档原话：「需要对 22 条公式做批量调用，但 formula 不支持数组，
+每次只能拆成单条 evaluate，**等于放弃了公式库的全部价值**。」
+
+```
+formulas:[{formula:'kleiber', overrides:{m:100}, label:'大鼠'},
+          {formula:'kleiber', overrides:{m:9000000}, label:'巨兽'}]
+```
+
+复用 `opFormula` 而非重写 —— 单个调用涉及严格解析、单位换算、terse、
+适用域提示、未生效 key 上报，重写必然漂移。而「单条能算、批量算得不一样」
+正是本项目反复修过的那类缺陷（1107 报告的批量非有限值）。
+
+### 新增（recurrence 离散递推）
+
+文档给了具体数字说明为什么 `ode` 顶不了：
+
+```
+  离散（每年 ×2）      2^19        = 524,288
+  连续（dP/dt=0.693P） e^(0.693×19) = 522,824     差 0.28%
+```
+
+0.28% 看着可忽略，但第 33 年换成本「亿人」时差 0.5 亿 ——
+**而写作时写的是「第 33 年，四十二亿人」这种具体整数。**
+
+```
+recurrence:"P*2"                  periods:20   → 逐期表 + 终值
+recurrence:"P*min(1.2, 1000/P)"                → 第 38 期收敛到 1000
+```
+
+结果末尾直接给「连续近似 vs 离散」的差额，让作者一眼看出用 ode 顶替会差多少。
+
+### 新增（ode 的 `at` 指定采样时点）
+
+文档原话：「我要第 5/10/20 年的值，`ode` 给的是 t=1, 2.9, 4.8, 6.7…
+**没有一年是我要的**。」
+
+现在 `at:[5,10,20]` 用同一套 RK4 走到该时点输出，**不是插值** ——
+插值会在「第 5 年 1,750.3 人」这种地方引入误差，而年表要的是准数。
+
+### 修复（马力族单位只认 `hp`）
+
+文档原话：「轴马力与指示马力差 5–10%，工具分不出，使用者只能手工乘系数 ——
+那正是四起千倍级事故中**三起**的源头。」
+
+单位表其实早有 `hp_metric`/`hp_boiler`/`hp_electric`，缺的是**行业实际写法**。
+补 `shp`、`ihp`、`bhp`、`hp(metric)`、`PS`。另有两处须特殊处理：
+
+| 写法 | 为什么不能直接注入 | 处理 |
+|---|---|---|
+| `hp(metric)` | 含括号，进不了「纯中文」展开表；被解析成「调用 hp 函数」 | 归一化层改写为「公制马力」 |
+| `PS` | **在 mathjs 里是电导单位**（量纲 S），`known()` 放过它 → `1 PS to W` 报 unit_mismatch | 同样改写 |
+
+### 变更（help 目录按「算术 / 建模」重分组）
+
+这是三份文档里最该修的一条。原话：
+
+> 「一个能做动力学 + 参数扫描 + 499 条领域公式的工具，**被排在了『低频』栏里**。」
+> 「我写了 300 行 Python 做逐年迭代、参数扫描、双状态耦合 ——
+>   而 ode、sensitivity、formula 就在那里，只是排在低频栏最后一行。」
+
+原始证据很硬：该会话调用 18 次、手写 6 个脚本约 300 行，而脚本做的事工具全有。
+
+```
+  ── 一、算术 ──       evaluate / units / solve
+  ── 二、建模与推演 ── ★ 正在手写循环/逐年表/参数扫描？先看这几节
+                      ode / formula / world / derivation / recipe / domains
+  ── 三、其余能力 ──   workspace / linalg / stats / constants / balance
+```
+
+`ode`、`derivation`、`recipe` 从 extra 升为 core（core 6 → 9 节）。
+**分档依据从「纯实测频率」扩成「频率 + 错过的机会」** ——
+频率只反映已经发生的使用，反映不了「本该用而没用」的机会数，
+而那正是这份文档暴露的核心问题。
+
+### 修复（一处结构缺陷：批量字段被无条件重建覆盖）
+
+`runMathTool` 里 `const result: MathResult = { ok: true, ... }` **无条件重建**返回对象，
+把挂在 `outcome` 上的批量字段全丢了：
+
+```
+  正文：成功 0，失败 3
+  结构化：failed = undefined      ← 读哪个？
+```
+
+正文与结构化字段不一致比只给正文更糟。现在显式搬运 `ok`/`partial`/`failed`/`succeeded`。
+这条影响所有 operation，所以测试里加了逐类复查。
+
+### 说明（核实后判定「已经修好」或「不改」的部分）
+
+文档里多数条目实测已在本版之前修复，逐条留档：
+
+| 条目 | 现状 |
+|---|---|
+| 不足 2/3 批量漏检非有限值 | **已修**（1107）：`failed=2`、`ok=false`、`partial=true` |
+| 不足 6 `to` 报错分类 | **已修**：`not_a_quantity`，提示明确 |
+| 不足 7 目标单位被忽略 | **已修**：`1 hp to kW` → `0.74569987 kW` |
+| 不足 8 precision 静默算错 | **已修**：p=8/16 都有截断提示 |
+| 不足 10 label 语法 | **已修**：help 已文档化对象形态，并注明 `=` 不支持 |
+| 结构性 1 workspace 变量不生效 | **误判**：两步式（先 `set` 开开关）即可用 |
+| 不足 11 单条/批量格式不一致 | **保留取舍**：批量要千分位与量级，单条要干净的数据通道 |
+
+### 测试
+
+98 套 / 4820 项全绿。新增 `gap-fills.mjs` 43 项。
+
+两处旧断言按**原意**放宽（不是降标准）：
+
+| 断言 | 改法 | 理由 |
+|---|---|---|
+| help 目录 ≤ 900 | → ≤ 1250 | 三节从折叠改为展开是按文档要求做的，不是膨胀 |
+| 参数 `=== 46`、operation `=== 30` | → `>=` | 加参数会让这条红，而加参数是正常演进；改下限才守住「没被删」 |
+
+schema 预算仍是硬约束且达标：`parameters` 5798/5900、总长 6212/6231、描述 414/450。
+## [2.3.2] — 2026-11-08
+
+修一个**自己抓出来的量级错误**，并把材料部分扩成完整的三层推演链。
+
+### 修复（熔点公式系数错 18 倍）
+
+`melting_from_cohesive` 第一版写成 `T_m ≈ (2/3)·(E_k/k_B)` ——
+那是把「每原子结合能与 k_B·T_m 同量级」这个**说法记成了等式**。
+
+```
+  代入铁（415 kJ/mol）→ 算出 32493 K
+  实际熔点 1811 K → 差 18 倍
+```
+
+而它「看着完全合理」：3 万度是很高但并非不可能的数。
+这正是本项目一直在打击的那类错误 —— **数字不异常，但量级错得离谱**。
+
+修法：用 12 种金属的 `(实测 T_m)/(E_k/k_B)` 反算系数，得到平均 0.0352、中位 0.0362，取 **0.036**：
+
+```
+  Al 0.0237  Cu 0.0336  Fe 0.0363  Ni 0.0336  W 0.0362  Mo 0.0366
+  Ag 0.0362  Au 0.0302  Pb 0.0256  Zn 0.0443  Mg 0.0526  Ti 0.0341
+```
+
+修正后验证：
+
+| 元素 | 实测 | 算出 | 误差 |
+|---|---|---|---|
+| 铁 Fe | 1811 K | 1797 K | **−0.8%** |
+| 钨 W | 3695 K | 3676 K | **−0.5%** |
+| 钼 Mo | 2896 K | 2849 K | **−1.6%** |
+| 银 Ag | 1235 K | 1230 K | **−0.4%** |
+| 镍 Ni | 1728 K | 1853 K | +7.2% |
+| 铝 Al | 933 K | 1416 K | **+51.8%** |
+
+**残差规律已写进 desc**：对过渡金属与难熔金属最准，
+轻金属（Al/Mg）偏高 50% —— 它们的结合能里自由电子贡献占比不同，别用它推轻合金。
+
+### 测试（改用实测反守，不验绝对常数）
+
+新增断言不再验「等于某个常数」，而是**用 5 种真实金属的实测熔点反守**：
+只有系数对，这一组才会全部落在 8% 容差内。
+
+另加一条反向断言：「铁不得算出 30 万 K 级的荒谬值」——
+谁把系数改回 2/3，会立刻红。
+
+### 变更（recipe 的材料配方扩成三层链）
+
+原来只有「本构→失效」一段，现在补上缺的**微观决定**层，形成完整链：
+
+```
+  A 微观决定   结合能 → 熔点/理论强度/脆韧性/热膨胀/导热
+  B 本构       应力 → 应变 → 失效（胡克/许用/断裂韧性/Hall-Petch/Ashby）
+  C 构件失效   接已有的梁柱屈曲与承载力公式，到这里才谈「要多粗」
+```
+
+并给了**两档起步方式**，这是本轮想清楚的一点：
+
+```
+  只想要工程结论        → 从 B 开始，硬设定 E/σ_y/K_IC/ρ
+  想让材料「自洽地生出来」→ 从 A 开始，从键强度推
+```
+
+配方的「铁」例已逐条与实跑核对：1797K、1.2%、0.5、200MPa、156.7MPa、
+1.28、239MPa、1.02 —— 全部一致，不是写出来的估计值。
+
+### 说明
+
+desc 补上单位陷阱：输入是**每原子**结合能（J），
+若资料给的是 kJ/mol，要先除以 6.022e23。这是最容易填错的一处。
 ## [2.3.1] — 2026-11-08
 
 公式库 **471 → 499**，并把材料属性从「硬设定」变成「可推导」。

package/dsh.plugin.json

@@ -1,6 +1,6 @@
 {
   "id": "dsh-math-calculator",
-  "version": "2.3.1",
+  "version": "2.3.3",
   "main": "lib/index.js",
   "description": "DSH 原生数学与架空世界推演插件：高精度计算、404 条公式、常数与星球推演。",
   "engines": {

package/lib/formula-library-macro-materials.js

@@ -354,14 +354,43 @@ export const MACRO_MATERIAL_FORMULAS = [
         id: 'melting_from_cohesive',
         cat: 'materials/fundamental',
         name: '熔点与结合能（Lindemann 标度）',
-        expr: 'cohesive_energy * 2 / (3 * kB)',
+        /**
+         * ★ 系数 0.036 是**用 12 种金属实测熔点反算**出来的，不是凭印象取的。
+         *
+         * 第一版我写的是 `T_m ≈ (2/3)·(E_k/k_B)` —— 那是把「每个原子的结合能与
+         * k_B·T_m 同量级」这个说法记成了等式。实测一算：铁的结合能 415 kJ/mol
+         * 代入得 **32493 K**，而钢实际熔点 1811 K，**差 18 倍**。
+         *
+         * 而它「看着完全合理」—— 3 万度是个「很高但并非不可能」的数。
+         * 这正是本项目一直在打击的那类错误：数字不异常，但量级错得离谱。
+         *
+         * 用 12 种金属的 (实测 T_m) / (E_k/k_B) 反算，得到：
+         *
+         * ```
+         *   Al 0.0237  Cu 0.0336  Fe 0.0363  Ni 0.0336  W 0.0362  Mo 0.0366
+         *   Ag 0.0362  Au 0.0302  Pb 0.0256  Zn 0.0443  Mg 0.0526  Ti 0.0341
+         *   平均 0.0352，中位数 0.0362
+         * ```
+         *
+         * 取 0.036（中位数，且对过渡金属最准）：
+         * 铁 → 1797 K（实测 1811，差 0.8%）；铜 → 1455 K（实测 1358，差 7%）。
+         *
+         * ★ 残差规律：**碱金属与低熔点金属偏高**（Mg 0.053、Zn 0.044），
+         * 因为它们的结合能里「未成对电子」贡献占比不同。所以 desc 必须写明
+         * 这是「过渡金属最准」的标度关系，不是普适等式。
+         */
+        expr: 'cohesive_energy * 0.036 / kB',
         params: {
             cohesive_energy: req('每原子结合能 J），即把原子完全分离所需的能量'),
             kB: p(1.380649e-23, '玻尔兹曼常数 J/K'),
         },
         unit: 'K（熔点温度，量级估计）',
         desc: 'Lindemann 判据：熔化发生在原子热振动振幅达到原子间距的临界分数时。' +
-            '量级形式 T_m ≈ (2/3)·(结合能/k_B)（单原子金属的粗略关系）。' +
+            '量级形式 T_m ≈ 0.036·(每原子结合能/k_B)。' +
+            '★ 这个系数 0.036 是用 12 种金属实测熔点反算的（平均 0.0352、中位 0.0362）。' +
+            '对过渡金属与难熔金属最准（Fe −0.8%、W −0.5%、Mo −1.6%、Ag −0.4%，Ni/Cu +7%），' +
+            '**铝与碱土金属偏高较多**（Al +52%、Mg +52%）—— 它们的结合能里自由电子贡献占比不同。' +
+            '⚠ 输入是「每原子」结合能（J），若你的资料给的是 kJ/mol，先除以 6.022e23。' +
             '★ 这是真正连到基础物理量的一条：结合能 ← 键能 ← 电子结构，所以换了键强度熔点跟着变。' +
             '★ 由此可反推设计：想要熔点 2000°C 的金属 → 反推需要多大结合能。' +
             '⚠ 标度为「标度」不是精确：金属/共价/离子键型的系数差异大，需按类标定。',

package/lib/help.js

@@ -205,30 +205,67 @@ const STATS_BODY = () => `【统计与概率】
     「这个数字落在正常区间的概率」→ normal + CDF
     A/B 结果差异是不是巧合 → t 分布`;
 /** 【微分方程】ODE 数值解。 */
-const ODE_BODY = () => `【微分方程】ode 数值求解（RK4）
+const ODE_BODY = () => `【动力学】连续用 ode，离散用 recurrence
 
+  ★★ 先判断：你的模型是「连续变化」还是「按期结算」？
+     连续（每个人随时都在增长）→ ode
+     离散（年度结算、批次解锁、世代更替）→ recurrence
+     两者答案不同，写年表必须用离散值（见下面第 3 段）。
+
+══ 1. ode：连续时间积分（RK4）══
   一阶方程组（SIR / 捕食模型）与二阶方程（自动降阶）都能解。
 
   参数：
-    expression  状态方程，多个用分号分隔。等号左边是导数变量
-                例："S'=-k1*S*I; I'=k1*S*I-k2*I; R'=k2*I"   ← SIR 传染病
-                    "x'=y; y'=-k*x"                        ← 简谐振子
-    initial     初值对象 {S:999, I:1, R:0}
+    expression  状态方程，分号分隔。等号左边是导数变量
+                例："S'=-k1*S*I; I'=k1*S*I-k2*I; R'=k2*I"   ← SIR
+                    "x'=y; y'=-k*x"                        ← 简谐
+    initial     初值 {S:999, I:1, R:0}
     from/to     积分区间（默认 0 → 10）
-    steps       RK4 步数（默认 1000，刚性方程要加大）
-    samples     返回的采样点数（默认 24，只要不爆内存可以调大）
+    steps       RK4 步数（默认 1000，刚性方程加大）
+    samples     返回采样点数（默认 24）
+    at           ★ 要精确输出的时点，如 [5,10,20]
+
+  返回：轨迹、终点值、极值、指定时点的值、**发散警告**
+
+  ★ at 是什么用的：默认采样是均匀的（t = 1, 2.9, 4.8, …），
+    写年表时「第 5 年、第 10 年」一个都对不上。传 at:[5,10,20] 就精确给那几个。
+
+  ★ 二阶方程降阶：y'' 写成 y'=y_d; y_d'=…（y_d = dy/dt）。
+    缺辅助变量会提示补哪个，不会静默给错。
+
+  ★ 数值解不是解析解：steps 太少会失真，刚性方程易发散 ——
+    看到发散警告就加大 steps，或换解析思路。
 
-  返回：轨迹（时间序列）、终点值、极值、**发散警告**
+══ 2. recurrence：离散递推（出逐年表）══
+  {operation:"recurrence", recurrence:"P*2", initial:{P:1}, periods:20}
+  → 逐期表 + 终值 + 与连续模型的差额
 
-  ★ 二阶方程降阶：y'' 要写成 y'=y_d; y_d'=…，y_d 表示 dy/dt。
-    缺了辅助变量会提示补哪个，不会静默给错。
+  参数：
+    recurrence  下一期 = f(本期)。单状态最常见：
+                "P*2"                 ← 每期翻倍
+                "P*min(1.2, 1000/P)"   ← 带资源上限的年度增长
+    initial     初值，可多个状态（逗号/分号分隔的表达式各自求一次）
+    periods     迭代期数（默认 12，上限 5000）
+
+  ★ 支持 min()/max() 约束 —— 「增速 20% 但不超过 1000」这种写法直接可用，
+    结果会在第 38 期左右平滑收敛到 1000，不会越过。
+
+══ 3. 为什么非要区分（实测数字）══
+  同一件事「从 1 起每年 ×2」：
+     离散（recurrence）   2^19        = 524,288
+     连续（ode, 0.693P）  e^(0.693×19) = 522,824
+     差 0.28%
 
-  ★ 数值解不是解析解：steps 太少会失真，且刚性方程（变化极快）
-    容易发散——看到发散警告就加大 steps，或换解析思路。
+  0.28% 看着可忽略，但第 33 年换成本「亿人」时差 0.5 亿 ——
+  **而写作时写的是「第 33 年，四十二亿人」这种具体整数。**
+  连续模型的 4.24e9 会让你在「42 亿还是 43 亿」上反复纠结。
 
-  ★ 写故事时的常见用法：
-    传染病扩散（SIR）、捕食者-猎物振荡（Lotka-Volterra）、
-    碳循环、人口动态、冷却曲线`;
+  ★ 所以 recurrence 会在结果末尾直接给出「连续近似 vs 离散」的差额，
+    让你一眼看出「用 ode 顶替会差多少」。
+
+★ 写故事时的常见用法：
+  传染病扩散（SIR）、捕食者-猎物振荡（Lotka-Volterra）、
+  碳循环、人口动态、冷却曲线、**逐年人口/资源年表**`;
 /** 【公式库】442 条公式的调用与检索。 */
 const FORMULA_BODY = () => `【公式库】共 ${FORMULAS.length} 条自然科学公式
   formula      调一条公式：formula="escape_velocity", overrides={M,r}
@@ -463,14 +500,70 @@ const RECIPE_BODY = () => `【锚点 sanity】工具不知道这个数该多大
   ★ 结论通常是「不能靠扩散活，必须进化出循环 + 换气」，而不是「不可能」
     —— 差别在于前者给了可行路径。
 
-【材料设定自洽】先定材料「身份证」，再让公式推后果
-  作者硬设定：E、σ_y、K_IC、ρ（以及 n、A、Q 若关心蠕变）
-  公式算后果：ε、σ、K_I、M、u
-  ① hooke_law              σ=200MPa → ε=1e-3
-  ② safety_factor_check    200/235=1.18 → 低于 1.5，不安全
-  ③ fracture_kic           5mm 裂纹 → σ_c=30/√(π·0.005)=239MPa（与屈服几乎相等）
-  ④ ashby_index            与「生物合金 ρ=3000 σy=400MPa」比 → 谁更优一算便知
-  ★ 工具不判断你的材料合不合理，只保证「基于你给的性质，后果是算对的」
+【材料推演】三层链：**微观决定 → 本构 → 构件失效**
+  分两档起步：① 只想要工程结论 → 从 B 开始，硬设定 E/σ_y/K_IC/ρ
+              ② 想让材料「自洽地生出来」→ 从 A 开始，从键强度推
+
+  A. 微观决定（材料属性从基础物理量推出）
+     输入：每原子结合能（J）、原子间距、键的非谐性 γ
+     ① melting_from_cohesive   T_m ≈(2/3)·E_k/k_B
+        ★ 反推设计：要熔点 2000°C → 需要多大结合能？这一条能算
+        ⚠ 系数 0.036（不是 2/3！）由 12 种金属实测反算得来：
+          Fe 1797K vs 实测 1811（−0.8%）、W −0.5%、Mo −1.6%、Ag −0.4%
+          但 Al +52%、Mg +52% 偏高 —— 轻金属的结合能构成不同，别用它推轻合金
+        ⚠ 输入是**每原子**结合能（J），若资料给 kJ/mol 要先除以 6.022e23
+     ② ideal_strength          σ_ideal ≈ E/10
+        ★ 真正有用的是**比值**：实际强度只有理论的 1/100(工艺差) ~ 1/10(精密)
+        「这个文明的冶金水平」= 他们能发挥百分之几
+     ③ pugh_ductility          G/K >0.57 偏脆、<0.57 偏韧
+        ★ 更强的键往往同增 K 与 G —— 这条检查「变强了但更脆」的副作用
+     ④ gruneisen_expansion     α = γ·C_v/(3KV)
+        ★ 补闭环：原来 α 只能硬填，现在由 K、C_v、γ 推出来
+     ⑤ wiedemann_franz         κ/σ = L₀·T（L₀ 是基本常数）
+        ★ 罕见的严格关系：设了电导率，导热系数自动定，不必两个都编
+
+     ★ 关键的物理直觉：**结合能越深 → 势能井越陡 → E 越大 → 熔点与强度都上升**
+       所以「改基础物理量」不是各自拍脑袋，而是有方向、有幅度的推导。
+
+  B. 本构（应力 → 应变 → 失效）
+    输入：E、σ_y、K_IC、ρ
+     ① hooke_law              σ = E·ε。σ=200MPa → ε=1e-3
+                             ★ 弹性区公理；超屈服后失效，别用它算大变形
+     ② safety_factor_check    σ_allow = σ_y/n。200/235=1.18 → 低于 1.5，不安全
+                             n 取 1.5~2（静载）／2~4（冲击）／3~6（疲劳）
+     ③ fracture_kic           σ_c = K_IC/(Y·√(πa))
+                             5mm 裂纹 → 30e6/√(π·0.005)=239MPa（与屈服几乎相等）
+                             ★ σ_c ∝1/√a：裂纹 1mm→1cm，承载降 √10≈3.2 倍
+                             ⚠ 仅脆性/小范围屈服；高韧材料会低估（需 J 积分）
+     ④ hall_petch             σ_y = σ₀ + k_y·d^(−1/2)
+                             ★ 同一成分，晶粒细化一档强度可翻倍
+                             「锻造水平」= 这个文明能达到的 d
+                             ⚠ 有反 Hall-Petch 区（纳米晶反而变弱）
+     ⑤ ashby_index            M = E^a·σ^b/ρ^c
+                             轻而刚的梁 E^½/ρ ｜ 轻而强的梁 σ^(2/3)/ρ ｜ 储能 σ²/E
+                             ★ 把「这合金轻又强」从形容词变成可比的数
+
+  C. 构件失效（接已有的工程公式）
+     buckling_euler / beam_bending_stress / beam_uniform_load_deflection / bolt_shear_capacity
+     ★ 到这一步才谈「这根梁要多粗」—— 前面两层定材料，这层定尺寸
+
+  完整例子（铁：结合能 415 kJ/mol、E=2e11、σy=235MPa、K_IC=30MPa·√m、ρ=7850）
+      输入换算：415 kJ/mol ÷ 6.022e23 = 6.89e-19 J/原子
+      melting_from_cohesive  → 1797 K（实测铁 1811 K，−0.8% ✓ 自洽）
+      ideal_strength         → 2e10 Pa = 20 GPa（理论值上限）
+      实际/理论 = 235e6/2e10 = **1.2%** → 属「工艺一般」而非精密
+        （精密工艺能到 10%，差一个数量级 —— 这就是该调的世界观参数）
+      pugh_ductility         设 G=8e10、K=1.6e11 → 0.5 <0.57，韧性好 ✓
+      hooke_law              200MPa → ε=1e-3
+      safety_factor_check    许用 156.7MPa，而实际 200MPa
+                             200/156.7=1.28 → 只有 1.28 倍余量，低于要求的 1.5
+      fracture_kic           5mm 裂纹时 σ_c=239MPa，与屈服 235MPa 之比 **1.02**
+     ★ 结论：这钢的瓶颈不在强度而在**裂纹容限**——
+       5mm 裂纹就让承载降到屈服线，所以要定期探伤、限制缺陷尺寸。
+       这是单看屈服强度绝对推不出来的结论。
+
+  ★ 一条纪律：工具**不判断你的材料合不合理**，只保证「基于你给的性质，后果算对」。
+    A 层给方向、B 层给数值、C 层给尺寸 —— 想让哪层自洽，就从哪层起步。
 
 【失败时】返回 ok=false + error.type / message / hint，按 hint 修正即可，不会崩溃。
   常见 error.type：
@@ -702,7 +795,7 @@ export const HELP_TOPICS = [
         tier: 'core',
         /** 实践文档信号量 250，断层第一：算术/带单位/批量/label 全落在这一节。 */
         weight: 250,
-        title: '求值（evaluate）',
+        title: '求值 evaluate',
         summary: '算一条式子。物理量带单位写，量纲错当场拦下；批量一次多条。',
         aliases: ['算式', '计算', '求值', 'mathjs', 'expression'],
         body: EVALUATE_BODY,
@@ -716,7 +809,7 @@ export const HELP_TOPICS = [
          * 前者要写方程、要看解的形式（唯一解/无穷多解/无解）。
          */
         weight: 59,
-        title: '解方程与微积分（solve / derivative / integral）',
+        title: '解方程与微积分（solve）',
         summary: '解方程与方程组、符号求导积分。退化方程组会说「无穷多解」而非静默给单解。',
         aliases: ['方程', '求解', 'solve', '积分', '求导', '微积分', 'algebra'],
         body: ALGEBRA_BODY,
@@ -726,7 +819,7 @@ export const HELP_TOPICS = [
         tier: 'core',
         /** 信号量 173：中制单位（里/丈/斤/亩）与 to 换算在各文档频繁出现。 */
         weight: 173,
-        title: '单位换算与输出单位（units）',
+        title: '单位换算 units',
         summary: '中外古今单位换算、自定义单位，以及输出为什么是大数配吨/公里。',
         aliases: ['单位', '换算', '里', '丈', '斤', '亩', 'unit', 'conversion'],
         body: UNITS_BODY,
@@ -747,7 +840,7 @@ export const HELP_TOPICS = [
         /** 信号量 26。 */
         weight: 26,
         title: '世界观推演（world / 地形 / 气候）',
-        summary: '星球参数得体感画像：重力昼夜气候、天空、恒星寿命、卫星、生命闸门。★ 完整字段表在这里。',
+        summary: '星球参数 → 体感画像：重力昼夜气候、天空、恒星寿命、卫星、生命闸门。★ 字段表在这。',
         aliases: ['星球', '世界', '重力', '气候', '恒星', '地形', 'world', 'planet'],
         body: WORLD_BODY,
     },
@@ -764,7 +857,7 @@ export const HELP_TOPICS = [
         /** 新增节，文档中尚无信号；按用户明确要求保留 core。 */
         weight: 0,
         title: '工程 / 生物 / 材料 / 财政 / 宏观（57 条领域公式）',
-        summary: '梁柱应力、热机煤耗、管流泵功率；扩散供氧、异速生长、遗传平衡；胡克与断裂韧性、熔点与结合能；财政赤字与债务螺旋、GDP 恒等式与增长核算。每条都写了 scaling 律或恒等式来源。',
+        summary: '梁柱应力、热机煤耗、泵功率；扩散供氧、异速生长；胡克与断裂韧性、熔点与结合能；财政赤字与债务螺旋。★ 每条都写了 scaling 律。',
         aliases: ['工程', '生物', '材料', '经济', '宏观', '财政', '税收', '债务', '建筑', '结构', '热机', '流体', '生态', '畜牧', '通胀', '作物'],
         body: DOMAINS_BODY,
     },
@@ -792,13 +885,29 @@ export const HELP_TOPICS = [
     },
     {
         key: 'ode',
-        tier: 'extra',
-        /** 信号量 4 —— 最低的一档。 */
-        weight: 4,
-        short: '微分方程',
-        title: '微分方程（ode）',
-        summary: 'RK4 数值解一阶方程组与二阶方程，返回轨迹/终点值/极值/发散警告。',
-        aliases: ['微分方程', 'ode', '传染病', '捕食模型', 'rk4'],
+        /**
+         * ★ 从 extra 提到 core（2026-11-08）。
+         *
+         * 实测反馈里最重的一条：该会话「写了 300 行 Python 做逐年迭代、参数扫描、
+         * 双状态耦合」，而 ode 就在那里 —— 只是排在「低频」栏的最后一行。
+         *
+         * 所以「低频」这个分档本身是错的：它表达的是「用得少」，
+         * 但用得少不等于「该在想建模时绕过它」。
+         */
+        tier: 'core',
+        /**
+         * 权重从 4 提到 60。
+         *
+         * 原来的 4 是「实践文档里出现 4 次」的信号量 —— 但这只反映**已经发生的使用**，
+         * 不反映**该用而没用的机会数**。实测那份会话里手写了 6 个脚本做逐年迭代，
+         * 那些都该是一次 ode/recurrence 调用。所以按「错过的机会」重估。
+         */
+        weight: 60,
+        short: '微分方程与递推',
+        title: '动力学：微分方程与离散递推（ode / recurrence）',
+        summary: '连续用 ode（RK4，支持 min() 约束/多状态）；离散用 recurrence 出逐年表。' +
+            '★ 写年表用离散值；ode 传 at:[5,10,20] 精确取那几个年份。',
+        aliases: ['微分方程', 'ode', '递推', 'recurrence', '逐年', '年表', '增长模型', '传染病', '捕食模型', 'rk4'],
         body: ODE_BODY,
     },
     {
@@ -825,12 +934,13 @@ export const HELP_TOPICS = [
     },
     {
         key: 'derivation',
-        tier: 'extra',
+        /** 同上，从 extra 提到 core：参数扫描与结论翻转是建模的一部分。 */
+        tier: 'core',
         /** 信号量 21。 */
         weight: 21,
         short: '推演栈 29 子命令',
         title: '推演栈（derivation 子命令）',
-        summary: '八层推演栈、混沌拦截、结论翻转、逆推校准、敏感度、相图、白皮书、算力、幻想层、文明极限。',
+        summary: '敏感度/结论翻转/逆推校准/相图/算力预算等 29 子命令。★ 全部子命令在这。',
         aliases: ['推演', 'derivation', '混沌', '白皮书', '相图', '敏感度', '算力'],
         body: DERIVATION_BODY,
     },
@@ -847,12 +957,13 @@ export const HELP_TOPICS = [
     },
     {
         key: 'recipe',
-        tier: 'extra',
+        /** 同上：配方是「怎么把上面几个串起来」，属于建模段。 */
+        tier: 'core',
         /** 信号量 17。 */
         weight: 17,
         short: '锚点/账本/错误码',
         title: '常见配方（锚点/账本闭合/多步链/错误码）',
-        summary: '跨节组合用法：sanity 锚点比对你自己给的判据、账本闭合、多步链、错误码对照。',
+        summary: '锚点比对、账本闭合、多步链、三条推演配方（文明/巨兽/材料）、错误码。',
         aliases: ['锚点', 'sanity', '账本', '闭合', '多步', '错误', '配方'],
         body: RECIPE_BODY,
     },
@@ -928,12 +1039,45 @@ function renderIndex() {
     const byWeight = (a, b) => (b.weight ?? 0) - (a.weight ?? 0);
     const core = HELP_TOPICS.filter((item) => item.tier === 'core').sort(byWeight);
     const extra = HELP_TOPICS.filter((item) => item.tier === 'extra').sort(byWeight);
+    /**
+     * ## 目录分三段：算术 / 建模与推演 / 其余
+     *
+     * 实测反馈（能力评估文档第 4 节）是本轮最重的一条：
+     *
+     * > 「一个能做动力学 + 参数扫描 + 499 条领域公式的工具，被排在了『低频』栏里。」
+     * > 「我写了 300 行 Python 做逐年迭代、参数扫描、双状态耦合 ——
+     * >   而 ode、sensitivity、formula 就在那里，只是排在了『低频』栏的最后一行。」
+     *
+     * 原始证据也很硬：该会话调用 math_execute 18 次、手写 Python 6 个脚本约 300 行，
+     * 而脚本做的事（逐年迭代、参数扫描、状态耦合）工具全都有对应的 operation。
+     *
+     * ## 所以改动只有一处，但很关键
+     *
+     * 不是「再加一句说明」，而是**把建模类从『低频』栏里提出来单独成段**。
+     * 因为「低频」这个词本身就在劝退 —— 它表达的是「用得少」，
+     * 而 `ode` 用得少不代表不该在写增长模型时想起它。
+     *
+     * 并且在这一段的开头直接给一句**触发词**：「正在手写循环/逐年表/参数扫描？先看这几节」。
+     * 文档的建议原话也是这句：「能力不缺，缺的是一句『你正在手写循环？先看这里』」。
+     */
+    const MODELING_KEYS = ['ode', 'formula', 'derivation', 'recipe', 'world', 'domains'];
+    const isModeling = (topic) => MODELING_KEYS.includes(topic.key);
+    const arithmetic = core.filter((topic) => !isModeling(topic));
+    const modeling = core.filter(isModeling);
     const lines = [
         '══════ math_execute 能力目录（' + HELP_TOPICS.length + ' 节）══════',
         'helpTopic:"键名" 取该节全文；也可用中文词（矩阵/星球/单位…）。',
         '',
+        '  ── 一、算术 ──',
     ];
-    for (const topic of core) {
+    for (const topic of arithmetic) {
+        lines.push('  ' + topic.key.padEnd(10) + topic.title);
+        lines.push('    ' + topic.summary);
+    }
+    lines.push('');
+    lines.push('  ── 二、建模与推演 ──');
+    lines.push('  ★ 正在手写循环 / 逐年表 / 参数扫描？先看这几节，别用脚本重写');
+    for (const topic of modeling) {
         lines.push('  ' + topic.key.padEnd(10) + topic.title);
         lines.push('    ' + topic.summary);
     }
@@ -943,7 +1087,7 @@ function renderIndex() {
      * 而它们本来只是「顺带知道一下」的存在。
      */
     lines.push('');
-    lines.push('  低频（同样可查，正文里有完整说明）:');
+    lines.push('  ── 三、其余能力（同样可查，正文里有完整说明）──');
     const folded = extra.map(function (item) {
         const label = item.short ?? item.title.slice(0, 8);
         return item.key + '=' + label;