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