[t y] = ode23('f501', [0 60], 100);
plot(t, y, '+', 'linewidth', 3);
xlabel('time'); ylabel('y value');
hold on;
plot(t, 90*exp(-0.1*t)+10, 'o', 'linewidth', 3); % Exact solution.
grid
hold off;
