K. Np and K. Np, size of an element M=21; % number of time points T=10; % time of integration % Initialization dt=T/(M-1); % time step size t=0:dt:T; % time discretization X=0:h:L; % spatial interval partitioning into elements x=zeros(K,Np); % grid matrix; first index stands for an element number, second for a point inside it, electric field vector E_=zeros(M,K*Np); % vector of coefficients of electric field expansion B=zeros(M,K*Np); % magnetic field vector B_=zeros(M,K*Np); % vector of coefficients of magnetic field expansion w=zeros(M,2*K*Np); % combined electric and magnetic fields vector w_=zeros(M,2*K*Np); % combined vector of coefficients S_=zeros(Np,Np); % stiffness matrix % Discretization matrices for electric and magnetic fields A_E=zeros

*. Np, Amplification matrix and its eigenvalues C=zeros

K. Np and K. Np, size of an element M=21; % number of time points T=0.56; % 10*0.9/0.05=10/18~=0.56; time of integration % Initialization dt=T/(M-1); % time step size t=0:dt:T; % time discretization X=0:h:L; % spatial interval partitioning into elements x=zeros(K,Np); % grid matrix; first index stands for an element number, second for a point inside it E=zeros(M,K*Np); % electric field vector E_=zeros(M,K*Np); % vector of coefficients of electric field expansion B=zeros(M,K*Np); % magnetic field vector B_=zeros(M, Np); % vector of coefficients of magnetic field expansion w=zeros(M,2*K*Np); % combined electric and magnetic fields vector w_=zeros(M,2*K*Np); % combined vector of coefficients S_=zeros(Np,Np); % stiffness matrix % Discretization matrices for electric and magnetic fields A_E=zeros

*. Np, Amplification matrix and its eigenvalues C=zeros

I. Unité-de-recherche and . Lorraine, Technopôle de Nancy-Brabois -Campus scientifique 615, rue du Jardin Botanique -BP 101 -54602 Villers-lès

I. De-voluceau-rocquencourt, BP 105 -78153 Le Chesnay Cedex (France) ?tt??????????r????r