t = linspace(0,10,1000);
c = -5:5;
m = length(c);

figure(1)
for i=1:m,
    plot(t,c(i)-exp(-t));
    hold on
end
title('Parallel solutions: yprime = exp(-t)')
xlabel('t')
ylabel('y')
figure(2)
t = linspace(0,2,1000);
for i=1:m,
    plot(t,c(i)*exp(-t));
    hold on
end
title('Converging solutions: yprime = -y')
xlabel('t')
ylabel('y')

figure(3)
for i=1:m,
    plot(t,c(i)*exp(t))
    hold on
end
title('Diverging solutions: yprime = y')
xlabel('t')
ylabel('y')