function dy=trajectory2_eq0(t,y,wn,te); 

%---2006.3.31----
%%%  nonlinear System

%---------------definition---------------------------------------
%     q'' +  2 te w q' + w^2 q=0
%                  
%     x1:=q   x2:=q'
%     x'1=x2
%     x'2=-2 te w x2 - w^2 x1                  
%-----------------------------------------------------------------
%----------------definition---   
%  y(1,1)   y(2,1)        
%           dy(1,1)        
%-----------------------------------------------------------------

global ni
dy=zeros(2,1);

dy(1,1)=y(2,1);
dy(2,1)= -2*te*wn*y(2,1) - wn^2 * y(1,1) ;







