numeric integration of kdv partial differential