/* ** constant -- program to calculate scaling constant for ** the CORDIC algorithm ** ** written by Joseph B. Evans, 11/15/91 ** evans@shannon.tisl.ukans.edu */ #include main(ac, av) int ac; char **av; { double scale; long i; int bits; if (ac == 2) { bits = atoi(av[1]); } else { printf("usage: %s bits\n",av[0]); exit(1); } for (i = 1, scale = 1.0; i <= bits; i++) { scale *= sqrt(1.0 + pow(2.0,(double)(-2*i))); printf("n = %d\t 1/K = %.9f\n",i,scale); } }