text: Title: Get costs at any given iteration
  Body: I want to visualise the convergence rate of k-modes (and k-prototypes, eventually) by plotting the cost function at each iteration. I've had a look through `kmodes.py` and I can't see a way of extracting that information directly. Could you suggest a way of going about this? 
  
  Thanks,
  
  Henry