▲ 另一组用于模拟计算的参数, L = 1，m1 = m2 = 1，k1 = k2 =1，G = 0.01
当引力较小时，近似于简谐运动。
▲ G = 0.01，k1 = 1，k2 = 5，m1 = m2 = 1，L = 1
当两个人魅力值相当时，劲度系数小的被吸引的距离会更远，但劲度系数大的运动周期会更短。
▲ G = 0.1，k1 = k2 = 1，m1 = 1, m2 = 0.5，L = 1
当引力较大时，可以看出此时已经不是简谐运动。
关于这个问题更多的细节和延伸，我们将其 remain as an open issue，欢迎大家积极讨论！
1.1.5
爱情摩擦力
两个物体接触面具有相对滑动趋势或已有相对滑动时，接触面上会产生阻碍相对运动趋势或阻碍相对滑动的力，这就是摩擦力。前者称为静摩擦力，后者称为滑动摩擦力。
固态物理间的静摩擦力大小是可变的，取决于受到的其他力的大小。但静摩擦力有上限。实验表明，最大静摩擦力与两物体接触面间的法向弹力（例如正压力与支持力）N 成正比，即
其中 μs 是静摩擦因数。
固态物体间的滑动摩擦力是一个定值，大小 f 与 N 成正比，有
