function answer=rklesseq2(m,M) % This program tells whether interval matrix px3 [m,M] with nonnegative first column % contains a matrix A such that A^w linear combination of A^v A^1 with v w = 2 3 or 3 2 % that is whether mrk(m,M) <=2 % answer=1 corresponds positive answer % answer=0 negative answer answer=0; a=sist(m(:,1),m(:,2),M(:,3),M(:,1),M(:,2),m(:,3)); b=sist(m(:,1),-M(:,2),M(:,3),M(:,1),-m(:,2),m(:,3)); c=sist(-M(:,1),m(:,2),M(:,3),-m(:,1),M(:,2),m(:,3)); d=sist(-M(:,1),-M(:,2),M(:,3),-m(:,1),-m(:,2),m(:,3)); if a+b+c+d>=1 answer=1; return else a1=sist(m(:,1),m(:,3),M(:,2),M(:,1),M(:,3),m(:,2)); b1=sist(m(:,1),-M(:,3),M(:,2),M(:,1),-m(:,3),m(:,2)); c1=sist(-M(:,1),m(:,3),M(:,2),-m(:,1),M(:,3),m(:,2)); d1=sist(-M(:,1),-M(:,3),M(:,2),-m(:,1),-m(:,3),m(:,2)); if a1+b1+c1+d1>=1 answer=1; end end