Get["simulineaire.m"]; 
circuitRLC[R_,L_,C_]:=Module[{A,B,n},
    n=4;
    A=Table[0,{n},{n}];
    B=Table[0,{n}];
    A=ajouterResistance[A,1,2,R];
    A=ajouterInductance[A,2,3,L];
    A=ajouterCapacite[A,3,4,C];
    {A,B}=ajouterMasse[A,B,4];
    {A,B}=definirEntree[A,B,1];
    Return[{A,B}];
]
{A,B}=circuitRLC[10,10^-3,10^-6];
                
h=transfert[A,B,3]
bodeGain[h,2,6,-80,20]
bodePhase[h,2,6]
reponseEchelon[h,2*10^-3,0,2]

{A,B}=circuitRLC[R,L,C];
h=transfert[A,B,3]/.s->I*omega 
                
