I want to produce samples from a double Schechter luminosity function.

Ideally, I am looking for a python implementation of this problem but am open to other tools if necessary.

Do you know of any such package that will facilitate this?

  • $\begingroup$ Do you mean you want to draw a random luminosity from a PDF given by a double Schechter? The rather old package CosmoloPy can return the cumulative distribution of a regular Schechter and other LFs, so perhaps you can modify and put in the extra term? $\endgroup$
    – pela
    Sep 23, 2020 at 0:14
  • 1
    $\begingroup$ That is what I meant, yes. Thank you for the link - looks promising! $\endgroup$
    – jm22b
    Sep 23, 2020 at 17:18


You must log in to answer this question.

Browse other questions tagged .