diff --git a/interfaces/matlab/toolbox/@Kinetics/fwdRateConstants.m b/interfaces/matlab/toolbox/@Kinetics/fwdRateConstants.m
new file mode 100644
index 000000000..5749913ab
--- /dev/null
+++ b/interfaces/matlab/toolbox/@Kinetics/fwdRateConstants.m
@@ -0,0 +1,9 @@
+function kf = fwdRateConstants(a)
+%FWDRATECONSTANTS Forward reaction rate constants for all the reactions.
+%
+%    kf = fwdRateConstants(a)
+%
+%        Returns a column vector of the forward rate constants of
+%        all of the reactions.
+%
+kf = kinetics_get(a.id,15,0);
diff --git a/interfaces/matlab/toolbox/@Kinetics/revRateConstants.m b/interfaces/matlab/toolbox/@Kinetics/revRateConstants.m
new file mode 100644
index 000000000..2a066b01e
--- /dev/null
+++ b/interfaces/matlab/toolbox/@Kinetics/revRateConstants.m
@@ -0,0 +1,9 @@
+function kr = revRateConstants(a)
+%REVRATECONSTANTS Reverse reaction rate constants for all the reactions.
+%
+%    kr = revRateConstants(a)
+%
+%        Returns a column vector of the reverse rate constants of
+%        all of the reactions.
+%
+kr = kinetics_get(a.id,16,0);
diff --git a/src/matlab/kineticsmethods.cpp b/src/matlab/kineticsmethods.cpp
index 43ac2e28b..c79bc2165 100644
--- a/src/matlab/kineticsmethods.cpp
+++ b/src/matlab/kineticsmethods.cpp
@@ -93,6 +93,12 @@ void kineticsmethods(int nlhs, mxArray* plhs[],
             case 14:
                 ok = kin_getEquilibriumConstants(kin,nr,h);
                 break;
+            case 15:
+                ok = kin_getFwdRateConstants(kin,nr,h);
+                break;
+            case 16:
+                ok = kin_getRevRateConstants(kin,1,nr,h);
+                break;
             default:
                 ;
             }