getf('simulineaire.sci');
function [A,B]=circuitRLC(R,L,C)
  N=4
  A=zeros(N,N)
  B=zeros(N,1)
  A=ajouter_resistance(A,1,2,R)
  A=ajouter_inductance(A,2,3,L)
  A=ajouter_capacite(A,3,4,C)
  [A,B]=ajouter_masse(A,B,4)
  [A,B]=definir_entree(A,B,1)
endfunction
            
[A,B]=circuitRLC(10,1e-3,1e-6);
[w,db,phi]=reponse(A,B,3,2,6,200);
            
plot1=scf();
plot2d(w,db,style=5,logflag='ln');
xtitle("Circuit RLC","log f","GdB");
            
plot2=scf();
plot2d(w,phi,style=5,logflag='ln');
xtitle("Circuit RLC","log f","phi");
            
