/* LIKELIHOOD FUNCTION FOR EXPONENTIAL DISTRIBUTION */

theta = seqa(0.05,0.05,150);
L = zeros(150,1);
loops = 1;

do while loops <= 150;
 	beta = (1/theta[loops,1]);
	L[loops,1] = (beta^5)*exp((-21.0)*beta);
  	loops = loops + 1;
endo;

library pgraph;
graphset;
title("Likelihood Function for Exponentially Distributed Data (n=5)"); 
xlabel("Beta");
xy(theta,L);