CtCn+1Cnt\frac{\partial \Psi_C}{\partial t} \approx \frac{\Psi_C^{n+1} - \Psi_C^n}{\Delta t} 2CCi+12Ci+Ci1(x)2\nabla^2 \Psi_C \approx \frac{\Psi_C^{i+1} - 2\Psi_C^i + \Psi_C^{i-1}}{(\Delta x)^2}
dengan t\Delta t sebagai langkah waktu dan x\Delta x sebagai langkah spasial.
2. Implementasi Algoritma Split-Step Fourier (SSF)
SSF adalah metode numerik yang cocok untuk sistem gelombang karena dapat menangani operator diferensial dengan transformasi Fourier.
Inisialisasi kondisi awal C(x,0)\Psi_C(x,0).
Evolusi dalam domain Fourier menggunakan eksponensiasi operator diferensial:
 C(k,t+t)=eiH^t/C(k,t)\Psi_C (k, t + \Delta t) = e^{-i \hat{H} \Delta t / \hbar} \Psi_C (k,t)Transformasi balik ke domain spasial menggunakan inversi Fourier.
3. Validasi Model dengan Eksperimen
Untuk membuktikan validitas model ini, kita dapat membandingkan hasil simulasi numerik dengan data dari: