Zunächst wird das package XAJ installiert.

```{r xajinstall}
install.packages("remotes")
remotes::install_github("Sibada/XAJ")
```

```{r xajinstall}
library(XAJ)
attach(SB)
area = 12
param.names = c('KC', 'IM', 'WUM', 'WLM', 'WDM', 'C', 'B', 'SM', 'EX', 'KI', 'KG', 'CI', 'CG', 'N', 'NK')

params=data.frame(lower = c(0.20, 0.00, 5.0, 10.0, 10.0, 0.05, 0.1, 10.0, 0.50, 0.01, 0.01, 0.50, 0.95, 0.1, 1.0), upper = c(1.50, 0.05, 20., 90.0, 60.0, 0.20, 0.6, 60.0, 2.00, 0.70, 0.70, 0.90, 0.998,5.0, 6.0), row.names = param.names)
params = (params$lower + params$upper)/2
Q <- XAJ(PREC, EVAP, params, area, dt = 24, full.UH = FALSE)
plot(Q$Q)
```