	function factor(DIS)
        factor = 0.0000001179*DIS**2+(-0.00003905)*DIS+0.984609
      return
      end
