How to use Gamfit function

Hi, I have written the code to get gamma function graph. Now i want to get the 95% confidence interval. So how to use Gamfit function to get it? Here's the code

X=1:50;

ALPHA=2.5;

BETA=0.6;

PDF=gampdf(X,ALPHA,BETA);

plot(X,PDF);