function z=dtds(s,t)
%z=1/sqrt( 1^2*sin(t)^2+1^2*cos(t)^2);
z=1/sqrt( (1.0)^2*sin(t)^2+(6.0)^2*cos(t)^2);