         double precision function   fe02rf
      *           (time ,uhigh ,ulow ,r0,delta ,r9,vrf,plz)
c time : aktuelle Simulationszeit
c uhigh : Ansprechspannung der Funkenstrecke
c ulow : Mindestspannung der Funkenstrecke (Kathodenfall)
c r0     : Anfangswiderstand beim Z"unden
c delta : Abklingkonstante
c r9     : Widerstand nicht gez"undet
C vrf    : Spannung , die an der Funkenstrecke anliegt
c plz : Zeitpunkt , zu dem Funkenstrecke gez"undet hat, sonst 0
         implicit real*8 (a-h,o-z),integer*4 (i-n)
         if (plz.gt.0.d0) then ! Funkenstrecke bereits gez"undet
            if (abs (vrf).lt.ulow)  then
               fe02rf = r9       ! Haltespannung unterschritten
               plz = 0.d0        ! Funkenstrecke gel"oscht
             else
               fe02rf = r0* exp ((plz-time)/delta)
            end if
         else                    ! noch nicht gez"undet
            if (abs (vrf).gt.uhigh)  then
               fe02rf = r0       ! Funkenstrecke spricht an
               plz = time        ! gez"undet
             else
               fe02rf = r9
            end if
         end if
         return
         end
