From d70c2904dd2e2d9407bc4b1393a6e9fdf2619cdf Mon Sep 17 00:00:00 2001
From: Neda Darbeheshti <neda.darbeheshti@aei.mpg.de>
Date: Wed, 18 Sep 2019 16:17:27 +0200
Subject: [PATCH] coefficients for case 12 was corrected
---
function_gfr/odeint_abm.m | 5 ++---
1 file changed, 2 insertions(+), 3 deletions(-)
diff --git a/function_gfr/odeint_abm.m b/function_gfr/odeint_abm.m
index d763496..4eee3d8 100644
--- a/function_gfr/odeint_abm.m
+++ b/function_gfr/odeint_abm.m
@@ -71,9 +71,8 @@ function [ y, dydt ] = odeint_abm(deriv, y0, dydt, t0, h, abm_order, deriv_para
beta_j = [0.280189596443937 -3.08887141086794 15.4861788582752 -46.6170361852654 93.6472204560486 -131.891420554754 133.019135965308 -96.2690500741543 49.2504906142276 -17.2688256657180 4.45198840045628];
beta_j_star = [-0.00678584998463470 0.0757510538586927 -0.385752772015793 1.18465362954946 -2.44382699765512 3.57154240820908 -3.80648324765512 3.01920720097803 -1.85839786130151 1.36990283957298 0.280189596443937];
case 12
- beta_j = [-18.1734761126059 60.7399929634890 -137.124664395693 220.357899950647 -258.602100049353 223.526764175736 -141.522864179368 64.3350953159655 -20.2857466060656 4.72625394048788];
- beta_j_star = [0.00592405641233766 -0.0719504705203490 0.401574156537264 -1.36322208005151 3.13959224562089 -5.18074106015512 6.30845647070907 -5.76142186372655 3.99667650901375 -2.18422096398008 1.43506746010869 0.274265540031599];
-
+ beta_j = [-0.274265540031599 3.29711053679153 -18.1734761126059 60.7399929634891 -137.124664395693 220.357899950647 -258.602100049353 223.526764175736 -141.522864179368 64.3350953159656 -20.2857466060656 4.72625394048788];
+ beta_j_star = [0.00592405641233766 -0.0719504705203490 0.401574156537264 -1.36322208005151 3.13959224562089 -5.18074106015512 6.30845647070908 -5.76142186372655 3.99667650901375 -2.18422096398008 1.43506746010869 0.274265540031599];
otherwise % use 8th order
beta_j = [-0.304224537037037 2.44516369047619 -8.61212797619048 17.3796544312169 -22.0277529761905 18.0545386904762 -9.52520667989418 3.58995535714286];
beta_j_star = [0.0113673941798942 -0.0938409391534392 0.343080357142857 -0.732035383597884 1.01796461640212 -1.00691964285714 1.15615906084656 0.304224537037037];
--
GitLab