1 files changed, 4 insertions, 3 deletions
diff --git a/regress.m b/regress.m
index ad06796..4d13c62 100644
--- a/regress.m
+++ b/regress.m
@@ -1,4 +1,5 @@
-function mn = regress(x, y)
-  mn = [ x ones(length(x), 1) ] \ y;
+function polynom = regress(x, y, degree)
+	pot = @(x) x.^[degree:-1:0];
+	M = reshape(cell2mat(arrayfun(pot, x(:)', 'UniformOutput', false)), degree+1, [])';
+	polynom = M \ y(:);
 endfunction
-