function dXdt = nni(t,x,lambda,mu,y0,y4) % nearest neighbour interaction (nni)
% y0 = linker äußerer Nachbar, konstant
% y4 = rechter äußerer Nachbar, konstant
 
dXdt = [exp(lambda*(X(2)-X(1)))+exp(mu*(X(1)-y0));exp(lambda*(X(3)-X(2)))+exp(mu*(X(2)-X(1))); exp(lambda*(y4-X(3)))+exp(X(3)-X(2)))];
 

% willkürliche Wahl: lambda=mu=1 und y0=0, y4=20 sowie y_1=5, y_2=10, y_3=15
% betrachte die ersten 100 Zeitschritte, d.h. <math>t\in [0,100]</math>
[t,X] = ode45(@(t,X) nni(t,X,1,1,0,20),[0 100], [5; 10; 15]);
 
% Plotte x_1=X(1), x(2)=X(2), x(3)=X(3) gegen t
 
plot(t,X(:,1),'-o',t,X(:,2),'-o',t,X(:,3),'-o')
 
xlabel('Time t');
ylabel('Solution x');
legend('x_1','x_2','x_3')