function fsch(x) ! for hot cathode
	common/cf/fmax
	fsch=fmax**(sqrt(x))
        return
        end
c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*