load nelm
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,8);
%ui(2:Q,1:4)=un(1:(Q-1),:);
%ui(3:Q,5:8)=un(1:(Q-2),:);
ui=un;

%ys=sim(net,ui')';
ys=seq2con(sim(net,con2seq(ui')));
ys=ys{1}';
plot(t,yn,'.',t,ys);
