'Create the full matrix `self.fullmat` and
the corresponding band matrix `self.bandmat`.'

def create_bandmat(self):

    N = 10
    self.KL = 2
    self.KU = 2
    self.sym_mat = ((((diag((1.0 * ones(N))) + diag(((-1.0) * ones((N - 1))), (-1))) + diag(((-1.0) * ones((N - 1))), 1)) + diag(((-2.0) * ones((N - 2))), (-2))) + diag(((-2.0) * ones((N - 2))), 2))
    self.herm_mat = ((((diag(((-1.0) * ones(N))) + (1j * diag((1.0 * ones((N - 1))), (-1)))) - (1j * diag((1.0 * ones((N - 1))), 1))) + diag(((-2.0) * ones((N - 2))), (-2))) + diag(((-2.0) * ones((N - 2))), 2))
    self.real_mat = ((((diag((1.0 * ones(N))) + diag(((-1.0) * ones((N - 1))), (-1))) + diag(((-3.0) * ones((N - 1))), 1)) + diag((2.0 * ones((N - 2))), (-2))) + diag(((-2.0) * ones((N - 2))), 2))
    self.comp_mat = (((((1j * diag((1.0 * ones(N)))) + diag(((-1.0) * ones((N - 1))), (-1))) + (1j * diag(((-3.0) * ones((N - 1))), 1))) + diag((2.0 * ones((N - 2))), (-2))) + diag(((-2.0) * ones((N - 2))), 2))
    (ew, ev) = linalg.eig(self.sym_mat)
    ew = ew.real
    args = argsort(ew)
    self.w_sym_lin = ew[args]
    self.evec_sym_lin = ev[:, args]
    (ew, ev) = linalg.eig(self.herm_mat)
    ew = ew.real
    args = argsort(ew)
    self.w_herm_lin = ew[args]
    self.evec_herm_lin = ev[:, args]
    LDAB = (self.KU + 1)
    self.bandmat_sym = zeros((LDAB, N), dtype=float)
    self.bandmat_herm = zeros((LDAB, N), dtype=complex)
    for i in xrange(LDAB):
        self.bandmat_sym[((LDAB - i) - 1), i:N] = diag(self.sym_mat, i)
        self.bandmat_herm[((LDAB - i) - 1), i:N] = diag(self.herm_mat, i)
    LDAB = (((2 * self.KL) + self.KU) + 1)
    self.bandmat_real = zeros((LDAB, N), dtype=float)
    self.bandmat_real[(2 * self.KL), :] = diag(self.real_mat)
    for i in xrange(self.KL):
        self.bandmat_real[(((2 * self.KL) - 1) - i), (i + 1):N] = diag(self.real_mat, (i + 1))
        self.bandmat_real[(((2 * self.KL) + 1) + i), 0:((N - 1) - i)] = diag(self.real_mat, ((- i) - 1))
    self.bandmat_comp = zeros((LDAB, N), dtype=complex)
    self.bandmat_comp[(2 * self.KL), :] = diag(self.comp_mat)
    for i in xrange(self.KL):
        self.bandmat_comp[(((2 * self.KL) - 1) - i), (i + 1):N] = diag(self.comp_mat, (i + 1))
        self.bandmat_comp[(((2 * self.KL) + 1) + i), 0:((N - 1) - i)] = diag(self.comp_mat, ((- i) - 1))
    self.b = (1.0 * arange(N))
    self.bc = (self.b * (1 + 1j))
