% roty_.m  rotate about the Y axis
function [xd, yd, zd] = roty_(x, y, z, th)
cosf = cos(th*pi/180);
sinf = sin(th*pi/180);
yd = y;
xd =  cosf.*x + sinf.*z;
zd = -sinf.*x + cosf.*z;