function Ry = rotyd(theta)
% Creates rotation matrix about axis Y
% input in degrees
theta=theta*pi/180;
ct=cos(theta);
st=sin(theta);
Ry=[ct 0 st 0;0 1 0 0;-st 0 ct 0;0 0 0 1];