Updated likelihood: now includes IMF term

- loop to calculate prob for star i now calls salpeter (or chabrier)
- include optional plots to diagnose variation of rho term, exp term
  and rho*exp term as a function of mass.
- at the end of the loop for star i, np.sum is used to add all the
  components k=0,...,K.