I'm trying to figure out how to get the entropy of a probability distribution. It's the standard differential entropy and the formula is: , where is the probability denstiy function.
I have referred to integral function q = integral(fun,xmin,xmax) I'm confused how to define this fun? For instance, I'm dealing with Beta distribution and I have used the following code for its pdf. However, this y is just the evaluated values of x, not exactly a fun. I suppose I can't use this y directly. So, is there a way that I can create or transfer y to a fun that I can apply with integral? Please help me with it. Thanks in advance!
pd = makedist('beta','a',50,'b',40)
x = linspace(0,1,10000);
y = pdf(pd,x);