function diffeq_sin_osc

close all
tic
y0=[1.5 3 4.5 6];

options = odeset('RelTol',1e-4,'AbsTol',[1e-4],'stats','on','outputfcn','odeplot');

[t,y]=ode15s(@sin_osc_solve,[0 100],y0);

plot(t,y,'-o')

figure

for j=1:length(y0)
    plot(y(:,j),'k','linewidth',j),hold on
end
toc


%w+a12*sin(y(2)-y(1))+a13*sin(y(3)-y(1)) +a14*sin(y(4)-y(1));