棄却法のアルゴリズム作成

Question

### 前提・実現したいこと> 引用テキスト

上記の画像が課題として課されたのですが、どのように解けば良いのか調べても分かりませんでした。
考え方の指針、問題の解法などどなたか教えて頂けませんでしょうか！

q(x)の一様分布の乱数を用いて、uの一様乱数を作成する所の解き方が特に分かりません。

Comments

[Penpen7]

解き方は、画像の中に十分ありますがお読みになりましたか？
それでもわからなければ、
https://nineties.github.io/math-seminar/14.html#/29
を見てはいかがですか？

[tiitoi]

一様乱数を作成する関数は random モジュールで提供されているので、それを使えばいいのではないでしょうか
さすがに疑似乱数列を作るところから実装する課題ではないですよね

Answer 1

棄却法

from random import random
import matplotlib.pyplot as plt

def P(n):
    x = []
    i = 0
    while(i < n):
        xprime = random()
        u = random() * 1.6
        if(u <= xprime * (1 - xprime) * 6):
            x.append(xprime)
            i += 1
    return x

fig = plt.figure()
ax = fig.add_subplot(111)
ax.hist(P(100000), bins=25, ec='black')
fig.show()

Comments

[退会済みユーザー]

回答ありがとうございます。

考え違いをしていましたので、回答に対する質問を撤回します.
Zはp(x)の積分であり、計算結果が1/6
それを割っているので、*6としていたのですね。