% generating of control actions function pol=mcontin1(mp,m,Ulower,Uupper) ul=Ulower; % lower bound of u uu=Uupper; % upper bound of u ur=(uu-ul)/m; cop=mp; rn=rand; pom=0; for f=1:m if (pomrn) uop=ul+(f-1)*ur; pol=uop; end pom=pom+cop(f); end