tempP=[];
for k = 1:4785
P = [y(k);ydot(k);tetha(k);tethadot(k)];
tempP = [tempP P];
end
net= newff([-20 20;-20 20;-20 20;-20 20],[60 1],{'tansig','purelin'},'trainlm');
net.trainParam.epochs = 100;
net = train(net,tempP,u(1:4785)');