library(MASS);
n <- 2 # Anzahl Bieter
vbar <- 10 # Höchstmögliche Wertschätzung der Bieter
## X~U[0,vbar]
F <- function(r) punif(r, min = 0, max = vbar, lower.tail = TRUE, log.p = FALSE)
f <- function(r) dunif(r, min = 0, max = vbar, log = FALSE)
G <- function(r) F(r)^(n-1)
g <- function(r) (n-1)*F(r)^(n-2)*f(r)
payment <- function(r) r*(1-F(r))*G(r) + area(function(y) y*(1-F(y))*g(y),r,vbar)
revenue <- function(r) n*payment(r)+F(r)^n*0 # Wertschätzung des Verkäufers: 0
revenue2 <- function(r) n*payment(r)+F(r)^n*2 # Wertschätzung des Verkäufers: 2
curve(revenue,xlab="r",ylab="Erlös",lwd=3,col="blue",xlim=c(0,11),ylim=c(0,5))
curve(revenue2,lwd=3,col="red",add=TRUE)