# Funktion zur Berechn. von BSB
# Variablen
ks = 0.07
kf = 0.30
BSB20 = 4
 
#Funktion definieren
BSB = function(k,t,BSB20){
BSB = BSB20*(1-10^(-k*t))
return(BSB)
}
 
# Funktion zur Temperaturkorrekur
tkorrektur = function(t=20,k20=0.07){
kc = k20*1.047^(t-20)
return(kc)
}
 
# Rechnen
t = seq(0,20,0.1)
Rs = BSB(ks,t,BSB20)
Rf = BSB(kf,t,BSB20)
 
# Temperaturkorrektur
kc = tkorrektur(t,kf)
plot(t,kc)