clear

N=7;
options=odeset('Reltol',1e-6,'OutputFcn','odeplot','OutputSel',[1,N]);

t1=500;
dt=0.5;



for S=1:N,
Y(S)=-1.5*pi*sin(pi*S/(2*N));
Y(S+N)=0;
end


[t,Y]=ode45('pendulatest1',[0:dt:t1],Y,options,N);


for cc=1:N, FL(cc)=0; JJ(cc)=(cc);   %JJ(NN+2)=(NN+1)/sqrt(10);
   for H=1:(t1/dt), hh(H)=H*dt;  
       YY(H,cc)=Y(H,cc);
       end
end