[watcher]

个人建议，直接离散化w，用for循环对每个w求解方程(solve())，最后应该得到与每个w对应的k序列，然后分列实部(real())虚部(imag())，分别绘图。
这样做虽然比较笨，但是应该可以解决问题，只是要注意，for循环变量不要用i、j，否则会与虚数(-1)^(1/2)冲突。