0001 function dif=compK1(arc,indQ1,indQ2,varargin) 0002 0003 KQ1=atgetfieldvalues(arc,indQ1,'PolynomB',{1,2}); 0004 KQ2=atgetfieldvalues(arc,indQ2,'PolynomB',{1,2}); 0005 0006 dif=KQ2(1)-KQ1(1); 0007 0008 return