function y=ali1(Aeq) pp=Aeq; clc v=zeros(10,9); for i=1:10 if pp(i,1)==1 v(i,1)=[1]; v(i,2)=[1]; v(i,3)=[1]; elseif pp(i,1)==2 v(i,1)=[-1]; v(i,2)=[1]; v(i,3)=[1]; elseif pp(i,1)==3 v(i,1)=[1]; v(i,2)=[-1]; v(i,3)=[1]; elseif pp(i,1)==4 v(i,1)=[1]; v(i,2)=[1]; v(i,3)=[-1]; elseif pp(i,1)==5 v(i,1)=[-1]; v(i,2)=[-1]; v(i,3)=[1]; elseif pp(i,1)==6 v(i,1)=[-1]; v(i,2)=[1]; v(i,3)=[-1]; elseif pp(i,1)==7 v(i,1)=[1]; v(i,2)=[-1]; v(i,3)=[-1]; else pp(i,1)==8 v(i,1)=[-1]; v(i,2)=[-1]; v(i,3)=[-1]; end if pp(i,2)==1 v(i,4)=[1]; v(i,5)=[1]; v(i,6)=[1]; elseif pp(i,2)==2 v(i,4)=[-1]; v(i,5)=[1]; v(i,6)=[1]; elseif pp(i,2)==3 v(i,4)=[1]; v(i,5)=[-1]; v(i,6)=[1]; elseif pp(i,2)==4 v(i,4)=[1]; v(i,5)=[1]; v(i,6)=[-1]; elseif pp(i,2)==5 v(i,4)=[-1]; v(i,5)=[-1]; v(i,6)=[1]; elseif pp(i,2)==6 v(i,4)=[-1]; v(i,5)=[1]; v(i,6)=[-1]; elseif pp(i,2)==7 v(i,4)=[1]; v(i,5)=[-1]; v(i,6)=[-1]; else pp(i,2)==8 v(i,4)=[-1]; v(i,5)=[-1]; v(i,6)=[-1]; end if pp(i,3)==1 v(i,7)=[1]; v(i,8)=[1]; v(i,9)=[1]; elseif pp(i,3)==2 v(i,7)=[-1]; v(i,8)=[1]; v(i,9)=[1]; elseif pp(i,3)==3 v(i,7)=[1]; v(i,8)=[-1]; v(i,9)=[1]; elseif pp(i,3)==4 v(i,7)=[1]; v(i,8)=[1]; v(i,9)=[-1]; elseif pp(i,3)==5 v(i,7)=[-1]; v(i,8)=[-1]; v(i,9)=[1]; elseif pp(i,3)==6 v(i,7)=[-1]; v(i,8)=[1]; v(i,9)=[-1]; elseif pp(i,3)==7 v(i,7)=[1]; v(i,8)=[-1]; v(i,9)=[-1]; else pp(i,3)==8 v(i,7)=[-1]; v(i,8)=[-1]; v(i,9)=[-1]; end if pp(i,4)==1 v(i,10)=[1]; v(i,11)=[1]; v(i,12)=[1]; elseif pp(i,4)==2 v(i,10)=[-1]; v(i,11)=[1]; v(i,12)=[1]; elseif pp(i,4)==3 v(i,10)=[1]; v(i,11)=[-1]; v(i,12)=[1]; elseif pp(i,4)==4 v(i,10)=[1]; v(i,11)=[1]; v(i,12)=[-1]; elseif pp(i,4)==5 v(i,10)=[-1]; v(i,11)=[-1]; v(i,12)=[1]; elseif pp(i,4)==6 v(i,10)=[-1]; v(i,11)=[1]; v(i,12)=[-1]; elseif pp(i,4)==7 v(i,10)=[1]; v(i,11)=[-1]; v(i,12)=[-1]; else pp(i,4)==8 v(i,10)=[-1]; v(i,11)=[-1]; v(i,12)=[-1]; end end for i=1:10 out(i)=sum(v(i,:)'); end aout = postmnmx(out,mint,maxt); y=(aout-b)/m;