Как правильно выбрать функцию для фитдиста в Matlab

Question

Я хочу выбрать правильную функцию, соответствующую моим данным. Вот что я пробовал:

myTolerance=1e-12; % in erg units.
nbins=100;
for j=1:ntM/10:ntM
   data = Wkinet(abs(Wkinet(:,j))>myTolerance, j) * erg2eV;
   pd = fitdist(data,'Normal');
   x = linspace(min(data),max(data),1000);
   y = pdf(pd,x); 
   loglog(x,y);
   %%hold on;
   pause;
end

Вот что я получаю:

Я хочу, чтобы это выглядело так:

В начале моделирования распределение является максвелловским, а в конце - степенным законом, как вы можете видеть.