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 ,
Amplification matrix and its eigenvalues C=zeros ,
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 ,
Amplification matrix and its eigenvalues C=zeros ,
Technopôle de Nancy-Brabois -Campus scientifique 615, rue du Jardin Botanique -BP 101 -54602 Villers-lès ,
BP 105 -78153 Le Chesnay Cedex (France) ?tt??????????r????r ,