Calculating own distribution of data set

I have data collected in some 150x1 vector. My goal is to calculate distribution of this values according to this formula, where α, β and B0 are some input parameters selected by user:

I wrote function to calculate distribution for whole data set, but the result is vector of 150x1 size filled with complex numbers, which is totally incoherent with my expected distribution. May I ask for help with correct calculation of this distribution?

function [y] = distribution(B,B0,alfa,beta)
GM = gamma(alfa);
for i = 1:length(B(:,1))
    y(i,1) = (((beta^alfa)/Y)*((B(i,1)-B0).^(-alfa-1)))*(exp(-(beta/(B(i,1)-B0))));
end

[darova]

Maybe is float number and is negative. What do you think?

Answer 1

[Sri Harsha Kondapalli]

Use histogram for obtaining a distribution function of the data set and then use curve fitting toolbox to estimate the parameters in your distribution.

https://www.mathworks.com/help/matlab/ref/hist.html

https://www.mathworks.com/help/curvefit/curve-fitting.html