
280 Green Chemical Engineering
if i == 4
n4 = n ;
E_thetahat4(1) = Ethetahat0 ;
end ;
if i == 5
n5 = n ;
E_thetahat5(1) = Ethetahat0 ;
end ;
for j = 1:500
thetahat = j*(thetahat_max/500) ;
theta_hat(j+1) = thetahat ;
Ethetahat = (n/factorial(n-1))*((n*thetahat)
∧
(n-1))*exp
(-n*thetahat) ;
if i == 1
E_thetahat1(j+1) = Ethetahat ;
end ;
if i == 2
E_thetahat2(j+1) = Ethetahat ;
end ;
if i == 3
E_thetahat3(j+1) = Ethetahat ;
end ;
if i == 4
E_thetahat4(j+1) = Ethetahat ;
end ;
if i == 5
E_thetahat5(j+1) = Ethetahat ;
end ;
end ;
end ;
plot(theta_hat,E_thetahat1,′g′,theta_hat,E_thetahat2,′r′,theta_hat,
E_thetahat3, ...