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:![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/243139/image.png)
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