#!/usr/bin/python
# -*- coding: utf-8 -*-


import numpy as np
from scipy import sparse

# Número de puntos de colocación
N = 100
dl = du = np.ones(N - 1)
d0 = -2 * np.ones(N - 1)
d = np.vstack((dl, d0, du))
A = sparse.spdiags(d, (-1, 0, 1), N - 1, N - 1)
print(A.todense())

"""
Metodos numericos para resolver sistemas:
   https://docs.scipy.org/doc/scipy/reference/sparse.html
   https://docs.scipy.org/doc/scipy/reference/generated/scipy.sparse.csr_matrix.html
   http://www.scipy-lectures.org/advanced/scipy_sparse/storage_schemes.html
   http://pysparse.sourceforge.net/spmatrix.html
   https://docs.scipy.org/doc/scipy/reference/tutorial/linalg.html
   https://docs.scipy.org/doc/scipy-0.13.0/reference/linalg.html
   https://docs.scipy.org/doc/numpy-1.13.0/reference/routines.linalg.html


"""