load narx
u=data(:,3:end);
y=data(:,1:2);
un=(u-meshgrid(um,1:size(u,1)))./meshgrid(r,1:size(u,1));
yn=(y-meshgrid(ym,1:size(y,1)))./meshgrid(s,1:size(y,1));

Q=size(u,1);
ui=zeros(Q,10);
ui(2:Q,1:2)=yn(1:(Q-1),:);
ui(2:Q,3:6)=un(1:(Q-1),:);
ui(3:Q,7:10)=un(1:(Q-2),:);

ys=sim(net,ui')';
plot(t,yn,'.',t,ys);