marginalekansY <- function(y){
  return((1/6) * (2+y) )
}

xwaarden <- seq(0,2,0.1)
ywaarden <- sapply(xwaarden, marginalekansY)

plot(xwaarden,ywaarden)



functieI <- function(n){

  xwaarden <- runif(n,min=0,max=2)
  ywaarden <- sapply(xwaarden,marginalekansY)
  return(ywaarden)

}