The gamma density function.
f(x)=x.^(a-1)exp(-x)./Gamma(a)1_{x>=0}
x=0:0.1:10; y=dgamma(x,3); plot(x,y)