% PLOT PRIOR GAMMA DISTRIBUTION
function GammaPrior(alpha,beta)
% Paramerers
% alpha: shape parameter
% beta: scale parameter

% set defaults if parameters are not passed to the function
if ~exist('alpha','var')
   alpha = 5;
end
if ~exist('beta','var')
   beta = 6;
end

figure(1);clf;
theta = [ 0:0.01:100 ];
prior = gampdf(theta,alpha,beta);
plot(theta,prior,'LineWidth',3);
