Ok, simple one: keep B positive.
0001 function B=ConstrainControlParameterB(B,m,params); 0002 0003 %Ok, simple one: keep B positive. 0004 0005 pos=find(B<=0); 0006 B(pos)=10; 0007