% 
% Initial condition variation example

clear all;
close all;

tspan = [0 4];  % time interval we would like to solve for

figure(1)
hold on
for i=1:80
    % Set our initial condition
    y0 = i/10-4;
    
    % Solve our ode using 4th Order Runge-Kutta
    [t,y] = ode45('ic_variation_rhs',tspan,y0);
    plot(t,y)
end

title('Solutions y(t)')
xlabel('t'); ylabel('y(t)')
hold off