# Flaeche vs. Saettigung
A <- seq(0.0,1,0.05) # Flaeche hier als Anteil zwischen 0 und 1
d <- A*0.0           # Defizit initialisiert

im <- 1.0            # maximale Speicherung, hier als Bereich 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)