%Programme pour une seule valeur de k (5 lignes qui suivent, sans "%") :
%global k
%k=0.9;
%[t,Y] = ode45('equa_diffexo2',[0 100],[1 0]);
%plot(t,Y(:,1))
%title(['Vibrations quand  k=',num2str(k)])

%Programme amélioré :
figure(1), clf
global k %la variable k prendra une valeur externe dans equa_diffexo2
parametres=[0,0.3,0.9,1.05,2]; %ensemble des valeurs décrites par k

hold on %les courbes vont être superposées
couleur=hsv(length(parametres)); %matrices 1x3 des couleurs des courbes
for i=1:length(parametres) %début de boucle
k=parametres(i);
[t,Y] = ode45('equa_diffexo2',[0 100],[1 0]);
plot(t,Y(:,1),'Color',couleur(i,:)) %la solution y est la colonne 1 de Y
text(1,5*i-2,['k=',num2str(parametres(i))],'Color',couleur(i,:))
end
hold off %les dessins suivants effaceront la fenêtre

title('Vibrations à l''approche de la résonnance (k=1)','FontSize',14)