*! 03-Aug-94.  STB-22: sg26.1
program define _jprfpgn
	local x `1'
	local lnx `2'
	local p `3'
	local np `4'
	local small 1e-6
/*
	Generate fractional powers.
	(Variables already exist in args 5, 6, ...)
*/
	quietly {
		local h0 1
		local plast 0
		local hlast "h0"
		local j 1
		while `j'<=`np' {
			local hj = `j'+4
			local pj = `p'[`j']
			if abs(`pj'-`plast')>`small' {
				if abs(`pj')<`small' {
					replace ``hj'' = `lnx'
				}
				else if abs(`pj'-1)<`small' {
					replace ``hj'' = `x'
				}
				else replace ``hj''=cond(`x'==0,0,`x'^`pj')
			}
			else {
				replace ``hj'' = `lnx'*``hlast''
			}
			local hlast `hj'
			local plast `pj'
			local j = `j'+1
		}
	}
end