x=linspace(0,2*pi,60);
y=sin(x);
z=cos(x);
t = sin(x) ./ (cos(x)+eps);