returns the rotation matrix around the y-axis with angle beta. This uses the same real matrix formual as in Elk - see Elk’s manual for ylmroty description