function rk2_test2

% Initial conditions
x0 = [0.003,0.000];  % State is expected to be a row vector.
t0 = 0;

% Set up for the integration with coarse steps
tfinal = 100;
Nsteps = 10000;
mu = [1e-1, 1e-6, 1e-8];
plot_colors = ['m','b','g']'; % Note that this is a column vector.

% Integrate for each viscosity
figure(1)
hold off
for n = 1:numel(mu)
   [x,t] = ode_fixed(@rk2, @(x,t) goofy(x,t,mu(n)), x0, t0, tfinal, Nsteps);
   plot(t, x(:,1), plot_colors(n))
   hold on
end
axis normal
legend('mu = 10^{-1}', 'mu = 10^{-6}', 'mu = 10^{-8}')
legend('location','southwest')
title('Nonlinear Damped Oscillator','fontsize', 20)
xlabel('t', 'fontsize', 20)
ylabel('x', 'fontsize', 20)
hold off