已知kdV方程：du/dt  + u(du/dx) +d^3u/dx^3 = 0    即：     (ut +uux+uxxx = 0)

用拟谱法对方程进行编程：即把函数u（x,t）定义区间标准化为[0,2pi]，在[0,2pi]区间取N个点，取空间步长为2pi/N在这些点上把函数u（x,t）变换到离散的傅里叶空间，然后再通过傅里叶逆变换把函数变回真实的物理空间。