fig=figure;
N=1000;
initPos=[];
for i=1:N/2
    initPos(i)=sin(i*4*pi/N);
end
for i=1:N/2
    initPos(i+N/2)=0;
end
initPos=initPos';
initVel=zeros(N,1);
IC=[initPos;initVel];
[t x]=RungeKuttaMethod(@string2,0,.01,2000,IC);
i=(1:N)';
[r c]=size(x);
for j=1:c
    plot(i,x(1:N,j))
    axis([0,N,-2,2]);
    drawnow;
end
    
