|
|
1.1 root 1: /*
2: * return 1 if args are within epsilon else 0
3: */
4:
5: double epsilon = 1.0e-20;
6:
7: approx(a, b)
8: double a, b;
9: {
10: if(0 > (a -= b))
11: a = -a;
12: return(a <= epsilon);
13: }
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.