1 files changed, 9 insertions, 0 deletions
diff --git a/R/SFORB.solution.R b/R/SFORB.solution.R
new file mode 100644
index 00000000..45a4533a
--- /dev/null
+++ b/R/SFORB.solution.R
@@ -0,0 +1,9 @@
+SFORB.solution = function(t, parent.0, k_12, k_21, k_1output) {
+  sqrt_exp = sqrt(1/4 * (k_12 + k_21 + k_1output)^2 + k_12 * k_21 - (k_12 + k_1output) * k_21)
+  b1 = 0.5 * (k_12 + k_21 + k_1output) + sqrt_exp
+  b2 = 0.5 * (k_12 + k_21 + k_1output) - sqrt_exp
+  
+  parent = parent.0 * 
+        (((k_12 + k_21 - b1)/(b2 - b1)) * exp(-b1 * t) +
+        ((k_12 + k_21 - b2)/(b1 - b2)) * exp(-b2 * t))
+}