#include <stdio.h>
#define n1 4096
#define n2 4032
char a[n2][n1],h[36][80],b[40];
FILE *fp;
main()
{
  int i,j,k,ix,iy,r;
// circle 2699,89,15
//        3977,761,90
  fp=fopen("ubad_3.old","rb");
  fread(h,36,80,fp); fread(a,4032,4096,fp); fclose(fp);
  ix=2699; iy=89; r=15;
  for(j=iy-r;j<=iy+r;j++)for(i=ix-r;i<=ix+r;i++){
    k=(j-iy)*(j-iy)+(i-ix)*(i-ix); if(k<=r*r)a[j][i]=1;
  }
  ix=3977; iy=761; r=90;
  for(j=iy-r;j<=iy+r;j++)for(i=ix-r;i<=ix+r;i++){
    k=(j-iy)*(j-iy)+(i-ix)*(i-ix); if(k<=r*r)a[j][i]=1;
  }
  fp=fopen("ubad_3.fit","wb");
  fwrite(h,36,80,fp); fwrite(a,4032,4096,fp); fclose(fp);
}