a,bを自然数に絞った手抜き版。 > import matplotlib > import pylab > x = [ 0.01*i for i in range(100)] > def minigamma(a): b = a-1 if b > 1: return b * minigamma(b) else: return 1 > def dbeta(x, a, b): return minigamma(a+b)/(minigamma(a)*min…