perc <- c(5,10,25,50,75,90,95); tit <- c("g.l.",perc)
gl <- c(1,2,3,4,5,10,20,50,100)
write(tit,file="",sep="\t",ncolumns=8); for (g in gl) write( c(g,signif(qchisq(perc/100,df=g),3)),file="",sep="\t",ncolumns=8)