from sympy import *
from sympy.physics.quantum.state import Ket, Bra
from sympy.physics.quantum import TensorProduct
import numpy as np
init_printing(use_unicode=True)

plus=Matrix([1/sqrt(2), 1/sqrt(2)])
zero=Matrix([1, 0])
minus=Matrix([1/sqrt(2), -1/sqrt(2)])

plus, zero, minus

TensorProduct(plus, zero, minus)