function th = theta(t)

if(t<1)
    th = 0;
else
    th = 10*pi/180;
end