In linear quantum optics we consider phase shifters, beam splitters, displacement operations, squeezing operations etc. The evolution can be described by unitary operators using Bose creation and annihilation operators. This evolution can be reduced to matrix multiplication using unitary matrices. We derive these evolutions for the different unitary operators. Finally a computer algebra implementation is provided.