function integrand=ps5p1_integrand(theta,z,xstar);

global GAMMA ALPHAX ALPHAZ ALPHA DELTA C;

integrand=exp(theta).*(1+(stdn_cdf(sqrt(ALPHAX)*(theta-xstar))).^GAMMA)...
   .*sqrt(ALPHA/(2*pi)).*exp(-0.5*ALPHA*(theta-DELTA*xstar-(1-DELTA)*z).^2);
          