#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 3687,3324,25
//        387, 2525,60
//        3901,1708,25
// matrix 4040,3700
//             3900    
//         32  2880
//             2970    
//         3300 3870
//         3645 4032
//
//         3613 3766
//          +66  +30
  fp=fopen("ubad_2.old","rb");
  fread(h,36,80,fp); fread(a,4032,4096,fp); fclose(fp);
  ix=3687; iy=3324; r=25;
  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=387; iy=2525; r=60;
  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=3901; iy=1708; r=25;
  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;
  }
  for(j=3700;j<=3900;j++)for(i=4040;i<4096;i++)a[j][i]=1; 
  for(j=2880;j<=2970;j++)for(i=0;i<32;i++)a[j][i]=1; 
  for(j=3870;j<=4030;j++)for(i=3300;i<3645;i++)a[j][i]=1; 
  for(j=3766;j<=3796;j++)for(i=3613;i<3679;i++)a[j][i]=1; 
  fp=fopen("ubad_2.fit","wb");
  fwrite(h,36,80,fp); fwrite(a,4032,4096,fp); fclose(fp);
}