constrain M45 by keeping X25 to sane values
better outputs for M45
add plot of cluster centroids