% coefficient for FDTD formula
function [ee,eh,hh,he]=Coefficient(eps,u,segema,dt)
temp=segema*dt./(2*eps);
ee=(1.-temp)./(1.+temp);
eh=(dt./eps)./(1.+temp);
temp2=segema*dt./(2*u);
hh=(1.-temp2)./(1.+temp2);
he=dt./u./(1.+temp2);