clear;
set (gca, 'fontsize', 24)

function xdot = f(x,t) 
  xdot(1) = 2*x(1) - x(2) + t;
  xdot(2) = x(1) - x(2)^2 - 1; 
end

x0 = [1, 1]';
t = linspace(0, 1, 51);
x = lsode('f', x0, t);

plot(t, x(:,1), 'b', 'Linewidth', 8, t, x(:,2), 'r', 'Linewidth', 4);
box off;
grid on;
xlabel 't';
ylabel 'x';
legend ('x1', 'x2', 'Location', 'NorthWest');

print -depsc2 -landscape lsode2.eps
print -djpg lsode2.jpg