```
# NOT RUN {
library(LaplacesDemon)
x <- dinvchisq(1,1,1)
x <- rinvchisq(10,1)
#Plot Probability Functions
x <- seq(from=0.1, to=5, by=0.01)
plot(x, dinvchisq(x,0.5,1), ylim=c(0,1), type="l", main="Probability Function",
ylab="density", col="red")
lines(x, dinvchisq(x,1,1), type="l", col="green")
lines(x, dinvchisq(x,5,1), type="l", col="blue")
legend(3, 0.9, expression(paste(nu==0.5, ", ", lambda==1),
paste(nu==1, ", ", lambda==1), paste(nu==5, ", ", lambda==1)),
lty=c(1,1,1), col=c("red","green","blue"))
# }
```

Run the code above in your browser using DataLab