mimo_dat
u=data(:,3:end);
y=data(:,1:2);
um=mean(u);
u_max=max(u);
u_min=min(u);
r=max(u_max-um,um-u_min);
un=(u-meshgrid(um,1:size(u,1)))./meshgrid(r,1:size(u,1));

ym=mean(y);
y_max=max(y);
y_min=min(y);
s=max(y_max-ym,ym-y_min);
yn=(y-meshgrid(ym,1:size(y,1)))./meshgrid(s,1:size(y,1));

ui=[yn(2:end-1,:) yn(3:end,:) un(2:end-1,:) un(3:end,:)];
yo=yn(1:end-2,:);
% net=newff(ones(12,1)*[-1 1],[5 2],{'tansig','tansig','purelin'},'trainlm');
% net=train(net,ui',yo');

ys=sim(net,ui')';
plot(t(1:end-2),[yo,ys]);
