mata:
function hf${i}(tnodes,expxb,tdexb,hvars,lt,time0)
{
	hf = ${matahazard${i}}
	return(hf)
}
end