global ALPHAX ALPHAZ ALPHA BETA DELTA ALPHATILDEZ B C X1STAR SELECT OPTIONS...
   THETA2LOWER THETA2UPPER THETA1STAR Z1STAR ZUNDERBAR ZUPPERBAR;

OPTIONS=optimset('Display','off');
ALPHAX=1;
ALPHAZ=3;
ALPHATILDEZ=ALPHAZ*ALPHAX;
ALPHA=ALPHAX*(1+ALPHAZ);
BETA=0.8;
DELTA=1/(1+ALPHAZ);
B=1;
C=0.5;
X1STAR=0;

if ALPHAZ*sqrt(ALPHAX)<=sqrt(2*pi)
   ZUNDERBAR=0.5;
   ZUPPERBAR=0.5;
else
   ZUNDERBAR=ps5p2_zunderbar;
   ZUPPERBAR=ps5p2_zupperbar;
end

SELECT=0;
x1star=fzero('ps5p2_stage1_integral',[0 0.5]);
X1STAR=x1star;
THETA1STAR=ps5p2_Theta(X1STAR);

figure(1);
clf;
hold on;
fplot('ps5p2_stage2_x2',[ZUNDERBAR-0.01 ZUPPERBAR+0.01]);
SELECT=1;
fplot('ps5p2_stage2_x2',[ZUNDERBAR-0.01 ZUPPERBAR+0.01],':');
SELECT=0;
xlabel('z');
ylabel('x_2*(z)');
YLIM=get(gca,'Ylim');
plot([ZUNDERBAR ZUNDERBAR],YLIM,'--');
plot([ZUPPERBAR ZUPPERBAR],YLIM,'--');
XLIM=get(gca,'Xlim');
plot(XLIM,[X1STAR X1STAR],'--');

filename=['c:\teaching\14.462 Spring 2004\problem sets\ps5_fig3'];
orient portrait;
print('-depsc',filename);
delete_student([filename '.eps']);

figure(2);
clf;
hold on;
fplot('ps5p2_stage2_xp',[ZUNDERBAR-0.01 ZUPPERBAR+0.01]);
SELECT=1;
fplot('ps5p2_stage2_xp',[ZUNDERBAR-0.01 ZUPPERBAR+0.01],':');
SELECT=0;
xlabel('z');
ylabel('x_p*(z)');
YLIM=get(gca,'Ylim');
plot([ZUNDERBAR ZUNDERBAR],YLIM,'--');
plot([ZUPPERBAR ZUPPERBAR],YLIM,'--');
XLIM=get(gca,'Xlim');
plot(XLIM,[X1STAR X1STAR],'--');

filename=['c:\teaching\14.462 Spring 2004\problem sets\ps5_fig4'];
orient portrait;
print('-depsc',filename);
delete_student([filename '.eps']);



SELECT=1;
x1star=fzero('ps5p2_stage1_integral',[0 0.5]);
X1STAR=x1star;
THETA1STAR=ps5p2_Theta(X1STAR);


figure(3);
clf;
hold on;
fplot('ps5p2_stage2_x2',[ZUNDERBAR-0.01 ZUPPERBAR+0.01]);
SELECT=0;
fplot('ps5p2_stage2_x2',[ZUNDERBAR-0.01 ZUPPERBAR+0.01],':');
SELECT=1;
xlabel('z');
ylabel('x_2*(z)');
YLIM=get(gca,'Ylim');
plot([ZUNDERBAR ZUNDERBAR],YLIM,'--');
plot([ZUPPERBAR ZUPPERBAR],YLIM,'--');
XLIM=get(gca,'Xlim');
plot(XLIM,[X1STAR X1STAR],'--');

filename=['c:\teaching\14.462 Spring 2004\problem sets\ps5_fig5'];
orient portrait;
print('-depsc',filename);
delete_student([filename '.eps']);


figure(4);
clf;
hold on;
fplot('ps5p2_stage2_xp',[ZUNDERBAR-0.01 ZUPPERBAR+0.01]);
SELECT=0;
fplot('ps5p2_stage2_xp',[ZUNDERBAR-0.01 ZUPPERBAR+0.01],':');
SELECT=1;
xlabel('z');
ylabel('x_p*(z)');
YLIM=get(gca,'Ylim');
plot([ZUNDERBAR ZUNDERBAR],YLIM,'--');
plot([ZUPPERBAR ZUPPERBAR],YLIM,'--');
XLIM=get(gca,'Xlim');
plot(XLIM,[X1STAR X1STAR],'--');

filename=['c:\teaching\14.462 Spring 2004\problem sets\ps5_fig6'];
orient portrait;
print('-depsc',filename);
delete_student([filename '.eps']);





