Spatially discretise the K-S PDE in space using 'holistic discretisation' (that has very good performance for nonlinear systems) via my web service at http://www.maths.adelaide.edu.au/anthony.roberts/holistic1.php
Then integrate the spatial discrtisation in time using Matlab ode15s or Octave lsode routines.