util1_jan/coord/go.m

   1 format compact\r
   2 % get coordinates\r
   3 load c\r
   4 m=c(1,1);\r
   5 n=c(1,2);\r
   6 c(1,:)=[];\r
   7 c1=c(:,1);\r
   8 c2=c(:,2);\r
   9 c1=reshape(c1,m,n);\r
  10 c2=reshape(c2,m,n);\r
  11 for i=1:2\r
  12     y1=1+100*rand;\r
  13     y2=1+200*rand;\r
  14     [x1,x2]=bint2_inv(c1,c2,bint2(c1,y1,y2),bint2(c2,y1,y2));\r
  15     err=[y1-x1,y2-x2]\r
  16 end\r