# Flaeche vs. Saettigung
A <- seq(0.0,1,0.05) # Flaeche zwischen 0 und 1
d <- A*0.0           # Defizit initialisiert
im <- 1.0            # maximale Speicherung, 0 bis 1
b <- 0.5             # Formfaktor
io = 0.5             # aktuelle Feuchte
i <- im*(1-(1-A)^{1/b})
for (n in 1:(length(A))) {
    if (i[n]>io) {
        d[n] <- (i[n]-io)
    } else {
        d[n] <- 0.0
    }
}
plot(A,d,xlim=c(0:1),ylim=(0:1))
curve((1-(1-x)^(1/1.5))-io,0,1,
        add = TRUE, col = "blue")
curve((1-(1-x)^(1/1))-io,0,1,
        add = TRUE, col = "red")
curve((1-(1-x)^(1/b))-io,0,1,
        add = TRUE, col = "gray")
curve((1-(1-x)^(1/0.2))-io,0,1,
        add = TRUE, col = "green")