diff --git a/lib/srfi/144.sld b/lib/srfi/144.sld
index 2323045e..f2384be7 100644
--- a/lib/srfi/144.sld
+++ b/lib/srfi/144.sld
@@ -20,13 +20,14 @@
                   (balanced-quotient flbalanced-quotient)
                   (balanced-remainder flbalanced-remainder)))
   (export
-   fl-e fl-e-2 fl-log2-e fl-log10-e fl-log-2 fl-log-10 fl-pi
-   fl-1/pi fl-2/pi fl-pi/2 fl-pi/4 fl-sqrt-pi fl-2/sqrt-pi
-   fl-sqrt-2 fl-sqrt-3 fl-sqrt-5 fl-sqrt-10 fl-1/sqrt-2
-   fl-cbrt-2 fl-cbrt-3 fl-4thrt-2 fl-phi fl-log-phi fl-1/log-phi
-   fl-euler fl-e-euler fl-sin-1 fl-cos-1 fl-greatest fl-least
-   fl-epsilon fl-integer-exponent-zero fl-integer-exponent-nan 
-   fl-fast-+*
+   fl-e fl-1/e fl-e-2 fl-e-pi/4 fl-log2-e fl-log10-e fl-log-2
+   fl-1/log-2 fl-log-3 fl-log-pi fl-log-10 fl-1/log-10 fl-pi
+   fl-1/pi fl-2pi fl-2/pi fl-pi/2 fl-pi/4 fl-pi-squared fl-degree
+   fl-sqrt-pi fl-2/sqrt-pi fl-sqrt-2 fl-sqrt-3 fl-sqrt-5 fl-sqrt-10
+   fl-1/sqrt-2 fl-cbrt-2 fl-cbrt-3 fl-4thrt-2 fl-phi fl-log-phi
+   fl-1/log-phi fl-euler fl-e-euler fl-sin-1 fl-cos-1 fl-gamma-1/2
+   fl-gamma-1/3 fl-gamma-2/3 fl-greatest fl-least fl-epsilon
+   fl-integer-exponent-zero fl-integer-exponent-nan fl-fast-+*
 
    flonum flonum? fl=? fl<? fl>? fl<=? fl>=? flodd? fleven?
    flunordered? flinteger? flzero? flpositive? flnegative?
diff --git a/lib/srfi/144/math.stub b/lib/srfi/144/math.stub
index d0abbcf0..a7f16703 100644
--- a/lib/srfi/144/math.stub
+++ b/lib/srfi/144/math.stub
@@ -3,18 +3,27 @@
 
 (define-c-const double
   (fl-e "M_E")
+  (fl-1/e "1.0/M_E")
   (fl-e-2 "logb(2.0)")
+  (fl-e-pi/4 "exp(M_PI/4.0)")
   (fl-log2-e "M_LOG2E")
   (fl-log10-e "M_LOG10E")
   (fl-log-2 "M_LN2")
+  (fl-1/log-2 "1.0/M_LN2")
+  (fl-log-3 "log(3.0)")
+  (fl-log-pi "log(M_PI)")
   (fl-log-10 "M_LN10")
+  (fl-1/log-10 "1.0/M_LN10")
   (fl-pi "M_PI")
   (fl-1/pi "M_1_PI")
-  (fl-2/pi "M_2_PI")
+  (fl-2pi "2.0*M_PI")
   (fl-pi/2 "M_PI/2")
   (fl-pi/4 "M_PI/4")
-  (fl-sqrt-pi "sqrt(M_PI)")
+  (fl-pi-squared "M_PI*M_PI")
+  (fl-degree "M_PI/180")
+  (fl-2/pi "M_2_PI")
   (fl-2/sqrt-pi "M_2_SQRTPI")
+  (fl-sqrt-pi "sqrt(M_PI)")
   (fl-sqrt-2 "M_SQRT2")
   (fl-sqrt-3 "sqrt(3.0)")
   (fl-sqrt-5 "sqrt(5.0)")
@@ -30,6 +39,9 @@
   (fl-e-euler "exp(0.57721566490153286060651209008240243)")
   (fl-sin-1 "sin(1.0)")
   (fl-cos-1 "cos(1.0)")
+  (fl-gamma-1/2 "tgamma(1.0/2.0)")
+  (fl-gamma-1/3 "tgamma(1.0/3.0)")
+  (fl-gamma-2/3 "tgamma(2.0/3.0)")
   (fl-greatest "DBL_MAX")
   (fl-least "-DBL_MAX")
   (fl-integer-exponent-zero "FP_ILOGB0")