今天是贝贝报告给你的第 242 天
2019-10-28 星期一

今天重构的代码跑出了结果，我又把结果画出了图，看起来一片混乱。

有这个进展，我高兴死了，理由是三点。

第一，我得出了结果。过去几天，我每次运行求解命令，有时马上，有时等半小时，得到的都是空集。我知道我的几十个方程中肯定有错误，但是我不确定它们在哪里。重构代码时，我重新思考了自己设定的各个条件之间的关系，排除了实际上是重复的约束，增加必要的条目，现在它一下子能给出我所有 36 个多项式系数的数值了，我能不高兴嘛。

第二，我能存储结果。重构代码后，我可以逐个运行各个函数，这样就知道在哪一步需要花费比较长的时间。除了改善算法，我还可以把不需要变化的中间结果存储到文件里。而从文件中读回来，是眨眼之间的事，节省了再次运算的二三十分钟。这样我可以为后续函数作更多次的试运行，我能不高兴嘛。

第三，我会理解结果。运算出来的结果只是数字，我需要把它们代回方程，再把表达式画出图像，才能理解结果是否靠谱。这次我没有转用数值的方式绘图，而是直接用符号运算原生的绘图工具，这样可以避免把复杂表达式函数化的时间。能画图，快画图，我就可以调整方程，向靠谱结果逼近了，我能不高兴嘛。

总之，现在的结果虽然错乱，过程却是可喜，我当然高兴啦。

图 1：明显有错啊，头尾不归零

图 2：有变化了吧，头尾归零了

图 3：调整了连续性条件，更不靠谱了

图 4：这是我要靠的那个谱 ​