x <- seq(0,8,length=100)
fx <- df(x,df1=1, df2 = 6)
lb = 2.26
ub = 10
#changing the UB to larger numbers doesnt do much after 10.

plot(x, fx, type="l", xlab="", ylab="")

i <- x >= lb & x <= ub
lines(x, fx)
polygon(c(lb,x[i],ub), c(0,fx[i],0), col="red")

area <- pf(ub,df1=4, df2=2) - pf(lb,df1=4, df2 = 2)
result <- paste("P(X >",lb,") =", signif(area, digits=3))

mtext(result,3)