curve(df(x, 3, 15, ncp=0), col="blue", xlim=c(0,5), ylim=c(0,1), xlab="f", ylab="f(f)")
text(x=0.03,y=0.7, label="df=(3,15)", cex=0.7)
curve(df(x, 5, 15, ncp=0), add=T, col="red", xlim=c(0,5), ylim=c(0,1), xlab="f", ylab="f(f)")
text(x=0.77,y=0.72, label="df=(5,15)", cex=0.7)
curve(df(x, 10, 15, ncp=0), add=T, col="green", xlim=c(0,5), ylim=c(0,1), xlab="f", ylab="f(f)")
text(x=1.14,y=0.9, label="df=(10,15)", cex=0.7)
arrows(1.12,0.88,0.85,0.78, length=.05)
curve(df(x, 15, 15, ncp=0), add=T, col="black", xlim=c(0,5), ylim=c(0,1), xlab="f", ylab="f(f)")
text(x=1.14,y=0.8, label="df=(15,15)", cex=0.7)

df1_list<-c(3,5,10,15)
df2_list<-c(15,15,15,15)
par(mfrow=c(2,2))


for (i in 1:length(df1_list)) {
  curve(df(x, df1=df1_list[i], df2=df2_list[i], ncp=0), xlim=c(0,5), ylim=c(0,1), xlab="f", ylab="f(f)", main=paste("df1=", df1_list[i], "df2=", df2_list[i]))
}