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 '20 at 0:14
  • 1
    $\begingroup$ That is what I meant, yes. Thank you for the link - looks promising! $\endgroup$ – jm22b Sep 23 '20 at 17:18

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.