program se
       implicit none
       real pi, n0, k0, il, z0, w0
       real c, fr, sd, freq, w, aseh, ste, asee, ff, fm
       integer i, ii
       complex j, sth
!
       pi = 4.*atan(1.)
       j = (0., 1.)
       c = 3.e+08
! not used:       e0 = (1.e-09)/36./pi
       n0 = 120.*pi
       fr = 1.e+06
       print *, 'Inseration Loss in dB>0 ?'
       read (*, *) il

       print *, 'Separation distance (m) between Magnetic dipoles?'
       read (*, *) z0

       print *, 'Separation distance (m) between Electric dipoles'
       read (*, *) w0
       sd = 2.*(10**(il/20.)-1.)/n0
       ff = 20.*alog10(1.+5*n0*sd)
       do i = 1, 3
          do ii = 1, 10
             if (ii==1 .and. i>1) cycle
             freq = ii*fr
             w = 2.*pi*freq
             k0 = w/c
             sth = j*k0*z0*n0*sd/6.
             aseh = 20.*alog10(cabs(1.+sth))
             ste = n0*sd/k0/w0
             asee = 20*alog10(ste)
             fm = freq/(1.0e+06)
             print '(1x, 4f9.2)', fm, aseh, asee, ff
          end do
          fr = 10.*fr
       end do
    end program se