--- a/src/g-builtins.cc
+++ b/src/g-builtins.cc
@@ -68,6 +68,7 @@
 #include "f-rand.h"
 #include "f-schur.h"
 #include "f-svd.h"
+#include "f-syl.h"
 
 #include "procstream.h"
 #include "error.h"
@@ -1573,6 +1574,23 @@
 }
 
 /*
+ * Sylvester equation solver.
+ */
+tree_constant *
+builtin_syl (tree_constant *args, int nargin, int nargout)
+{
+  tree_constant *retval = NULL_TREE_CONST;
+
+  if ((nargin == 4) && (nargout == 1))
+    DLD_BUILTIN (args, nargin, nargout, syl,
+		 retval = syl (args, nargin, nargout);)
+  else
+    usage ("X = syl (A,B,C)");
+
+  return retval;
+}
+
+/*
  * Schur Decomposition
  */
 tree_constant *