em_edipole_con_e_impulse_cylinder.m

function s = em_edipole_con_e_impulse_cylinder(r, z, t, con)
% author: Ainray
% date: 20210826
% email: wwzhang0421@163.com
% introduction: h filed of electrical dipole in cylinder coordinates in
%               conducting media
%   funciton name:  em, electromagnetic
%                   edipole, electric dipole
%                   con, conductive medium
%                   h, magnetic field
%                   impulse, impulse response
%                   cylinder, cylinder coordinates
    if(t(1) <eps)
        t(1)=eps;
    end
    mu0 = 4*pi*1e-7;
    eta = sqrt(mu0*con./t/4);
    rr2 = r.*r + z.*z;

    eta2rr2 = eta.^2*rr2;
    eta2r2 = eta.^2*r*r;
    eta2rz = eta.^2*r*z;
    eta3 = eta.^3;
    
    coef = eta3*pi^(-1.5)/con./t;
    coefe = coef.*exp(-eta2rr2);
    
    er = coefe.* eta2rz;
    ez = coefe.*(1-eta2r2);
    s(:,1) = er;
    s(:,2) = ez;
    s(1,:) = 0;