#include <stdio.h>
#include <math.h>
int main(void) {
double x =1e15;
double x1 =1e16;
double result1
= sqrt(x1
+ 1) - sqrt(x1
);
printf("sqrt(x + 1) - sqrt(x)"); printf("result(10^15):%.25f\n",result
); printf("result(10^16):%.25f\n",result1
);
double x2 =1e15;
double x3 =1e16;
double result2a
= 1 / (sqrt(x2
+ 1) + sqrt(x2
)); double result2b
= 1 / (sqrt(x3
+ 1) + sqrt(x3
));
printf("1 / (sqrt(x + 1) + sqrt(x))"); printf("result2(10^15):%.25f\n",result2a
); printf("result2(10^16):%.25f\n",result2b
); return 0;
}
I2luY2x1ZGUgPHN0ZGlvLmg+CiNpbmNsdWRlIDxtYXRoLmg+CgppbnQgbWFpbih2b2lkKSB7CgkKCWRvdWJsZSB4ID0xZTE1OwoJZG91YmxlIHgxID0xZTE2OwogCglkb3VibGUgcmVzdWx0ID0gc3FydCh4ICsgMSkgLSBzcXJ0KHgpOwoJZG91YmxlIHJlc3VsdDEgPSBzcXJ0KHgxICsgMSkgLSBzcXJ0KHgxKTsKCQoJcHJpbnRmKCJzcXJ0KHggKyAxKSAtIHNxcnQoeCkiKTsKCXByaW50ZigicmVzdWx0KDEwXjE1KTolLjI1ZlxuIixyZXN1bHQpOwoJcHJpbnRmKCJyZXN1bHQoMTBeMTYpOiUuMjVmXG4iLHJlc3VsdDEpOwoJcHJpbnRmKCJcbiIpOwoJCglkb3VibGUgeDIgPTFlMTU7Cglkb3VibGUgeDMgPTFlMTY7CgkKCWRvdWJsZSByZXN1bHQyYT0gMSAvIChzcXJ0KHgyICsgMSkgKyBzcXJ0KHgyKSk7Cglkb3VibGUgcmVzdWx0MmI9IDEgLyAoc3FydCh4MyArIDEpICsgc3FydCh4MykpOwoJCglwcmludGYoIjEgLyAoc3FydCh4ICsgMSkgKyBzcXJ0KHgpKSIpOwoJcHJpbnRmKCJyZXN1bHQyKDEwXjE1KTolLjI1ZlxuIixyZXN1bHQyYSk7CglwcmludGYoInJlc3VsdDIoMTBeMTYpOiUuMjVmXG4iLHJlc3VsdDJiKTsKCXJldHVybiAwOwp9Cg==