this is the code i have tried for the
Exploring the route to measure synchronization in non-linearly coupled Hamiltonian systems journal..but its not running..
any advice or suggestion is welcome.
%% %% plotting coupling measure schronization for N=2
clear all; close all; clc; %#ok<CLALL>
y=[0:0.01:40];
m=length(y);
%tspan=[0:0.1:100];
E=0.2;
p1=0.2;
p2=sqrt(2*E - (p1)^2)
for i=1:m
[t{i},x{i}]= ode45(@(t,x) exploring(t,x,y(i)),[0 100],[0 0 p1 p2]);
end
%EX=(1/2).*(y*X(:,2).^2 * X(:,1).^2)
%S= trapz(EX);
%%
figure(1)
plot(X(:,1),X(:,3),'b',X(:,2),X(:,4),'r ')%y=0 together
figure(2)
subplot(1,2,1)
plot(X(:,1),X(:,3),'b');%y=0 seperate plots
subplot(1,2,2)
plot(X(:,2),X(:,4),'r');
figure(3)%y=0.007 together
plot(Y(:,1),Y(:,3),'b',Y(:,2),Y(:,4),'r ')
figure(4)%y=0.007 seperate plots
subplot(1,2,1)
plot(Y(:,1),Y(:,3),'b');
subplot(1,2,2)
plot(Y(:,2),Y(:,4),'r');
figure(5)%y=0.15 together
plot(W(:,1),W(:,3),'b',W(:,2),W(:,4),'r ')
figure(6)%y=0.15 seperate plots
subplot(1,2,1)
plot(W(:,1),W(:,3),'b');
subplot(1,2,2)
plot(W(:,2),W(:,4),'r');
figure(7)%y=0.52 together
plot(Z(:,1),Z(:,3),'b',Z(:,2),Z(:,4),'r ')
figure(8)%y=0.52 seperate plots
subplot(1,2,1)
plot(Z(:,1),Z(:,3),'b');
subplot(1,2,2)
plot(Z(:,2),Z(:,4),'r');
