diff options
Diffstat (limited to 'R/HS.solution.R')
| -rw-r--r-- | R/HS.solution.R | 6 |
1 files changed, 6 insertions, 0 deletions
diff --git a/R/HS.solution.R b/R/HS.solution.R new file mode 100644 index 00000000..4651a6a8 --- /dev/null +++ b/R/HS.solution.R @@ -0,0 +1,6 @@ +HS.solution <- function(t, parent.0, k1, k2, tb)
+{
+ parent = ifelse(t <= tb,
+ parent.0 * exp(-k1 * t),
+ parent.0 * exp(-k1 * tb) * exp(-k2 * (t - tb)))
+}
|