x <- seq(0,50,length=100)
fx <- df(x,df1=1, df2 = 9)
lb = 49.85
ub = 100

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)