coefficients for case 12 was corrected

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];