'''
Created on Nov 24, 2015
@author: ahanagrawal
'''
import numpy as np
from OperationsA import gaussElimin
if __name__ == '__main__':
a = np.array([[1,0,0],[1,1,0],[1,5/3,1]])
b = np.array([[1,2,4],[0,3,21],[0,0,0]])
print(np.linalg.det(a))
print(np.dot(a,b))
c = np.array([1,2,5])
x = np.linalg.solve(a, c)
print(x)
print(a,c)
print(np.dot(a,x))
x = gaussElimin(a, c)
print(x)
print(np.dot(a,x))
if __name__ == "__main__":
a = np.array([[1, 1.2, 1.4, 1.6, 1.8, 2.0]])
A = np.array([])
b = np.array([[0,1,0,1,0,1]])
copy = np.zeros((1,6))
print("Len is " + str(len(a)))
for i in range(getLen(a)):
for j in range(getLen(a)):
copy[0,j] = a[0,i]**(a.size-1-j)
if (A.size is 0):
A = np.copy(copy);
else:
A = np.append(A, copy, 0)
print("break_____" + "\n",A)
print("final")
print(A)
print(vandermode(np.array([1, 1.2, 1.4, 1.6, 1.8, 2.0])))
x = gaussElimin(A, b.T)
# print(A)
print(x)
# x = gaussian_elim.GEPP(A, b.T)
for i in range(6):
print(np.dot(A[i,:],x))
