Change < to <= in getMu_vec. Add verification in laplacianSmoothing for new position.
Add calcRHS.f90, Euler only