sim('opt021s');
plot(1-t,K1,'k', 1-t, K2, 'k--', 1-t, K3, 'k:')
legend('r = 0.01','r = 0.02', 'r = 0.10')
xlabel('t')
ylabel('P')
set(1, 'PaperUnits', 'centimeters')
set(1, 'PaperPosition', [3 10 11 8.25])
print -deps 'opt021.eps'
