function tr=tps_reponse(valeur, temps, val_final, taux)
%entree : valeur : tableau des amplitudes
%         temps : tableau des pas de temps
%         val_final : valeur finale
%         taux pour le cacul du temps de reponse (0.05)

seuilh= val_final*(1+taux); 
seuilb= val_final*(1-taux); 

pos=find((valeur < seuilh)&(valeur > seuilb));
ref=[length(temps)-length(pos)+1:1:length(temps)]';
comp=ref-pos;
pos_comp=find(comp >0);
tr = temps(pos(length(pos_comp)+1));
