getting probability distribution function of output

hi,

how could i get probability disribution of qult instead of frequency on a Y-axis??for example by hist(qult)command u have that but i would like to have PDF of qult.is there anyway to get that probability distribution of qult/

my code:

B=1000;L=2000;Df=200;Sc=1.1;Sq=1.1;Sgama=.8;dgama=1;dc=1.08;
nsamples=10000;
   for i=1:nsamples
       CS=C(1,i);
       WS=G(1,i); 
       F=fi(1,i);
       Nq=tan((pi/4)+(pi*F/360))*tan((pi/4)+(pi*F/360))*2.718^(pi*tan(F*pi/180));
       Nc=(Nq-1)*cot(F*pi/180);
       Ngama=2*(Nq+1)*tan(F*pi/180);
       dq=1+(0.4)*(tan(F*pi/180))*(1-sin(F*pi/180))^2;
       qult(i)=(CS*Nc*Sc*dc)+(384*Nq*Sq*dq)+(980*Ngama*Sgama*dgama);
   end