function u=nn_pid(e,f)

global um r ym s net en1 fn1

en=(e-um)/r;
fn=(f-ym)/s;
un=sim(net,[fn1;en;en1])';
en1=en;
fn1=fn;
u=ym+un*s;
