format long e n=100; for d1=-5:5 for d2=-5:5 for i=1:100 for j=1:100 k((i-1)*100+j)=(d1+i/n)+2*(d2+j/n)+0.03*sin((d1+i/n)+(d2+j/n))-4; t((i-1)*100+j)=5*(d1+i/n)+6*(d2+j/n)+0.07*cos((d1+i/n)-(d2+j/n))-8; if (abs(k((i-1)*100+j)) < 10^(-2)) if (abs(t((i-1)*100+j)) < 10^(-2)) x1=d1+(i/n) x2=d2+(j/n) end end end end end end