function outmat = QWP; outmat = [1 0;0 i];