From 3ba4a2ea40c8ee2c3019778658c662ed30829c9f Mon Sep 17 00:00:00 2001 From: Benjamin Franzke Date: Mon, 14 May 2012 08:44:08 +0200 Subject: regress: Add an ugly workaround hack for old octave versions --- regress.m | 13 +++++++++++-- 1 file changed, 11 insertions(+), 2 deletions(-) (limited to 'regress.m') diff --git a/regress.m b/regress.m index 4d13c62..89a4e0b 100644 --- a/regress.m +++ b/regress.m @@ -1,5 +1,14 @@ function polynom = regress(x, y, degree) - pot = @(x) x.^[degree:-1:0]; - M = reshape(cell2mat(arrayfun(pot, x(:)', 'UniformOutput', false)), degree+1, [])'; + assert(degree == 1 || degree == 2); + % doesnt work with old octave versions currently + %pot = @(x) x.^[degree:-1:0]; + %M = reshape(cell2mat(arrayfun(pot, x(:)', 'UniformOutput', false)), degree+1, [])'; + M = []; + if degree == 1 + M = [ x.^1 x.^0 ]; + end + if degree == 2 + M = [ x.^2 x.^1 x.^0 ]; + end polynom = M \ y(:); endfunction -- cgit