Moved linspace function into module

include/linalg.h

@@ -0,0 +1,10 @@
+#include <gsl/gsl_matrix.h>
+
+#ifndef LINALG_H
+#define LINALG_H
+
+int
+linspace
+( double start, double end, int length, gsl_matrix *out );
+
+#endif

src/linalg.c

@@ -0,0 +1,12 @@
+#include "linalg.h"
+
+int linspace ( double start, double end, int length, gsl_matrix *out ) {
+  if ( out->size2 != length )
+    return 1;
+
+  double by = end / length;
+  for ( int i = 0 ; i < length ; i++ )
+    (void)gsl_matrix_set(out, 0, i, i*by);
+
+  return 0;
+}