clear
close all
attr=['r-';'g-'; 'b-'; 'c-'; 'm-'; 'k-'];

kxi=[0.05 0.1 1/sqrt(2) 1 1.2];
G=1;
omega0=1;
%Courbe de réponse
pul=[0:0.00001:1.8];
x=pul/omega0;

for ik=1 : length(kxi)
Ax(ik,:)=abs(G)./sqrt((1-x.*x).^2+4*kxi(ik)*kxi(ik).*x.*x);
end

for ik=1 : length(kxi)
figure
hold on
plot(pul, Ax(ik,:))
grid on
xlabel('Pulsation')
ylabel('H_F(w)')
legend(['kxi= ' num2str(kxi(ik))])
end

disp('Calcul de la fonction de transfert pour une pulsation donnée');
disp('Indiquer la valeur du facteur d''amortissement puis celle de la pulsation');
disp('Entrer 0 pour terminer les calculs');

encore=1;
while encore ~=0
    val_kxi=input('Facteur d''amortissement : ');
    if val_kxi~=0
        omega = input('Valeur de la pulsation : ') ;
        if omega~=0
            x=omega/omega0;
            out=abs(G)./sqrt((1-x.*x).^2+4*val_kxi*val_kxi.*x.*x);
           sprintf('Amplitude  %0.5g',out) 
        else
            encore=0;
        end    
    else
       encore =0; 
    end    
end