# -*- coding: utf-8 -*-
"""
Created on Sun Nov 06 08:49:00 2016

@author: wderousse
"""

import sys
import xlrd
import numpy as np
import matplotlib.pyplot as plt
from scipy.interpolate import interp1d

np.set_printoptions(threshold=np.nan)
print sys.version
print __name__

filename='linear_4.xlsx'
print
print 'Opening',filename
ss=xlrd.open_workbook(filename,on_demand=True)
sheet_index=0
print ' Worksheet',ss.sheet_names()[sheet_index]
sh=ss.sheet_by_index(sheet_index)

print ' Reading data values'
xa=[]
for row_index in range(0,sh.nrows):
    row_single=[]
    for col_index in range(0,sh.ncols):
        row_single.append(sh.cell_value(row_index,col_index))
    xa.append(row_single)
xa=np.matrix(xa)
print ' Design x array size:',len(xa)
print
print 'Adjacency Matrix:'
print xa

k=2
nt=xa**2
print
print 'For ',k,' points:'
print nt