add yule simon rng

[lingium]

Description

This issue is a feature proposal.

As a continuation of #3219, we will add rng for yule simon distribution.

Expected Output

ints yule_simon_rng(reals alpha)

Checklist

This will require a yule_simon_rng with analytical derivatives.

• rng
• test of rng

Current Version:

v5.1.0