title: I can't plot the kernel density estimation of an Inverse gamma(0.001, 0.001) in matlab
body: <p>I am trying to plot the ksdensity of an Inverse <code>gamma(0.001, 0.001)</code> but the plot has only one point. The commands I used are</p>
  
  <pre><code> alpha1 = 0.001;
      beta1 = 0.001;
      n = 1e+5;
      r=1./gamrnd(alpha1,1/beta1,n,1);  
     [f,xi] = ksdensity(r);
      plot(xi,f,'--m');
  </code></pre>
  
  <p>The first term of <code>f</code> is a real number and all the others are <code>NaN</code>. The first term of <code>xi</code> is a real number and all the others are <code>Inf</code>.</p>
  
  <p>Could you please me help me with this.</p>
  
  <p>Thank you very much.</p>