summaryrefslogtreecommitdiff
path: root/regress.m
blob: 89a4e0b793cf4f345b83c11356d9ed665c3363ba (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
function polynom = regress(x, y, degree)
	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