function outmat = HWP; outmat = [1 0;0 -1];