看起来你可能想从这个讨论中窃取一些关于在 Julia 中实现类似东西的想法。因为有很多不同的细节，所以很难清楚地总结，但本质上它是一种算法，它做一个简单的网格，然后自适应地在线性近似较差的区域中递归地添加点。这就像有限元方法中的课程网格自适应性，您需要更多的网格点并查看您的插值在网格点处偏离了多远，如果它很多，您接受新的网格点并递归地重新应用粗化操作. 然后是细节，但它都是开源的，有相关的讨论，所以你可以查看所有血淋淋的细节。从 1D 到 2D 可能涉及新代码，但原理相同。

答案可能是这个https://scicomp.stackexchange.com/a/13048/25606

听起来您想要 SymPy ** 的 ** 多精度绘图函数，它能够在给定范围内绘制任意黑盒函数。