x<-(0:2500)/1000 clip<-function(x,y) ifelse(x>2,0,ifelse(x<0,0,y)) plot(x, clip(x,(3/8)*(2-x)^2), type="l", xlab="x", ylab="dens(X=x)", xlim=c(0,2)) points(x, clip(x,(2-x)/2), type="l") points(x, clip(x,0.5), pch=".")