fix: plot_function补上跟calculator一致的安全沙箱

之前plot_function直接把expressions丢给sp.sympify，没有走_check_expr_safe
白名单/黑名单校验，也没有超时保护——跟calculator用的是同一个sympify攻击面，
却是唯一没设防的一条。现在每条表达式先过_check_expr_safe，通过后在原有的
_CALC_POOL子进程池里求值（15秒超时），跟calculator复用同一套沙箱机制。
被拒绝/超时/解析失败的表达式会在图上标注对应文字，不影响其他正常表达式
继续画出来。已用真实的正常/恶意/混合三种输入验证过，68个既有测试全过。

Author: heliotrope

tools.py

@@ -805,6 +805,21 @@ def _run_step_decomposer(problem_type: str, problem: str) -> str:
 # 4. 可视化工具实现
 # ─────────────────────────────────────────
 
+def _plot_eval_impl(expr_str: str, xmin: float, xmax: float) -> tuple:
+    """在子进程中解析并求值单条曲线——跟 calculator 走同一个安全沙箱/进程池/超时，
+    避免 plot_function 成为绕开 _check_expr_safe 的另一个攻击面。"""
+    import numpy as np
+    import sympy as sp
+    x_sym = sp.Symbol("x")
+    sym_expr = sp.sympify(expr_str.replace("^", "**"))
+    fn = sp.lambdify(x_sym, sym_expr, "numpy")
+    x_arr = np.linspace(xmin, xmax, 2000)
+    y_arr = np.array(fn(x_arr), dtype=complex)
+    y_arr = np.real(y_arr)
+    y_arr[np.abs(y_arr) > 1e5] = np.nan
+    return x_arr.tolist(), y_arr.tolist(), sp.latex(sym_expr)
+
+
 def _run_plot_function(expressions, xmin=-10, xmax=10, ymin=None, ymax=None,
                        title="", labels=None) -> str:
     try:
@@ -839,22 +854,27 @@ def _run_plot_function(expressions, xmin=-10, xmax=10, ymin=None, ymax=None,
         # 多条曲线用实线/虚线/点划线区分，颜色只用深色系
         line_styles = ["-", "--", "-.", ":"]
         colors = ["#1a1a1a", "#2255aa", "#cc2200", "#228833", "#aa44bb"]
-        x_sym = sp.Symbol("x")
-        x_arr = np.linspace(float(xmin), float(xmax), 2000)
+        _xmin, _xmax = float(xmin), float(xmax)
 
         for i, expr_str in enumerate(expressions):
+            err = _check_expr_safe(expr_str)
+            if err:
+                ax.text(0.05, 0.95 - i * 0.07, f"表达式被拒绝: {expr_str}（{err}）",
+                        transform=ax.transAxes, fontsize=9, color="red")
+                continue
             try:
-                sym_expr = sp.sympify(expr_str.replace("^", "**"))
-                fn = sp.lambdify(x_sym, sym_expr, "numpy")
-                y_arr = np.array(fn(x_arr), dtype=complex)
-                y_arr = np.real(y_arr)
-                y_arr[np.abs(y_arr) > 1e5] = np.nan
+                fut = _CALC_POOL.submit(_plot_eval_impl, expr_str, _xmin, _xmax)
+                x_list, y_list, latex_expr = fut.result(timeout=_CALC_TIMEOUT)
+                x_arr, y_arr = np.array(x_list), np.array(y_list)
                 lbl = (labels[i] if labels and i < len(labels)
-                       else f"$y = {sp.latex(sym_expr)}$")
+                       else f"$y = {latex_expr}$")
                 ax.plot(x_arr, y_arr,
                         color=colors[i % len(colors)],
                         linestyle=line_styles[i % len(line_styles)],
                         lw=1.8, label=lbl)
+            except _FutTimeout:
+                ax.text(0.05, 0.95 - i * 0.07, f"计算超时: {expr_str}",
+                        transform=ax.transAxes, fontsize=9, color="red")
             except Exception as e:
                 ax.text(0.05, 0.95 - i * 0.07, f"解析失败: {expr_str}",
                         transform=ax.transAxes, fontsize=9, color="red")