#include <math.h>
#include <stdio.h>
#include <fcntl.h>

#define CYGWIN 1

int fd;

iwrite(fd,n,l)
int fd,l;
unsigned int n;
{
   write(fd,&n,l);
}

writebin()
{
#ifdef CYGWIN
   fd = open("out.bin",O_WRONLY|O_CREAT|O_TRUNC|O_BINARY,S_IRUSR|S_IWUSR|S_IRGRP);
#else
   fd = open("out.bin",O_WRONLY|O_CREAT|O_TRUNC|O_BINARY);
#endif
   if (fd < 0) return(-1);

   int i,j,a, wtlen , sam;

   a = (int) (pow(2,23) -1)/2;
   
   /* wtlen = 1024*128; */
   wtlen = 1024*128;
   for (i=0; i<wtlen;i++) {
      sam = 0;
      for (j=1; j<48; j++) {
         sam +=  (int) (1.0/j * a*sin( j *  2.0 *3.1415926535*((double) i)/((double) wtlen) ));

      }
/*     printf("%d\n",sam); */
     iwrite(fd,sam,4);
   }

   close(fd);
}


main()
{
   writebin();
}