t=0:0.1:4;
x0=1;
r=1;
gamma = sqrt(1 + 1 / r);
x= x0 * exp( - gamma*t);
u = - x0/(r*(1 + gamma )) * exp( - gamma* t);
plot(t,u,'k-',t,x,'k--')
legend('u^*','x^*')
xlabel('t')
ylabel('u^*, x^*')
set(1, 'PaperUnits', 'centimeters')
set(1, 'PaperPosition', [3 10 11 8.25])
print -deps 'opt012.eps'
