c=======================================================================

      function plkavg(wvnmlo,wvnmhi,t) 3,1
      wl=20000./(wvnmlo+wvnmhi)
      plkavg=planck(wl,t)
      return
      end