varified marching_global.f90 with Matlab (remember the periodic boundary along x)