import numpy as np
import scipy.linalg as alg
import matplotlib.pyplot as pl

A = alg.pascal(4)
print(A)
print(np.linalg.matrix_power(A, -1))
print(alg.inv(A))
print(alg.lstsq(A, np.identity(4))[0])

x = np.linspace(-2*np.pi, 2*np.pi, 100)
y = np.cos(x)
pl.plot(x, y)

N = 6
M = np.zeros((N, N))
for i in range(1, N):
    for j in range (1, N):
        M[i, j] = np.round(np.sqrt(i*i + j*j), 2)

print(M)