--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/vandermonde.m	Tue Feb 20 15:00:30 2018 +0100
@@ -0,0 +1,10 @@
+% Create vandermonde matrix for points x and polynomials of order p
+% x and p are vectors
+% v is a length(x) by length(p) matrix
+function V = vandermonde(x, p)
+    V = sym(zeros(length(x), length(p))); % Is there a way to make this work for both double and sym
+
+    for i = 1:length(p)
+        V(:, i) = mononomial(x,p(i));
+    end
+end