// Multi-evaluate func meval( void ) { f = open("meval.scm","w"); rstab,"get_sources.scm",2,ip,jp,typ="ii"; nrows = numberof(ip); im = readfits("m1.fits+0"); for( i = 1; i <= nrows; i++ ) { evaluate, im, [ip(i),jp(i)],5,9,p=1; w0; im(ip(i)-3:ip(i)+3,jp(i)-3:jp(i)+3) = 0; ans = rdline(prompt="Next ... "); } }