Ok, simple one: keep drag positive.
0001 function drag=ConstrainControlParameterdrag(drag,m,params); 0002 %Ok, simple one: keep drag positive. 0003 0004 pos=find(drag<=0); 0005 drag(pos)=5;