ttt=[];
order=3;

for theta=0:0.1:3*pi
tt=[0 1];
for order_n=2:order
    tt(order_n+1)=(2*order-1)/(order-1)*cos(theta)*tt(order_n)-order/(order-1)*tt(order_n-1);
end
t=order*cos(theta)*tt(order+1)-(order+1)*tt(order);
tt=tt(end);

ttt=[ttt t];
end