Do you mean something like this
%plot n vs x n vs y
forn=1:10
[t,r]=ode45(@fn,[0 14],[0.01+n 0.01 ]);
xlast(n) = r(end,1);
ylast(n) = r(end,2);
end
plot(1:n,xlast,'o',1:n,ylast,'s'),grid
xlabel('n'),ylabel('x,y')
legend('last x','last y')
functiondr=fn(~,r)
x=r(1); y=r(2);
dr=zeros(2,1);
dr(1)=x+y;
dr(2)=x-y;
end
