x = [1, 2, 3]
y = [1, 3, 7]
deg = 2
'''
x = [-2,0,1]
y = [-27,-1,0]
deg = 2
x = [0,1,2,3,4,5,6,7,8,9,10]
y = [0,1,2,3,4,5,6,7,8,9,10]
deg = 1
'''
print('Normal equn')
reg = Regression(x, y, deg)
s = reg.regNormalEq(x, y, deg)
yfit = reg.evalReg(s, x)
#print(yfit)
output.coeffPlotScreen(x, y, yfit, s)
print('Aug sys')
reg = Regression(x, y, deg)
s = reg.regAugSys(x, y, deg)
yfit = reg.evalReg(s, [2, 3])
print(yfit)
print('QR')
reg = Regression(x, y, deg)
s = reg.regOrthogonalTn(x, y)
yfit = reg.evalReg(s, [2, 3])
print(yfit)
output.coeffFile(s)
