Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #include <stdio.h>
- double sqrt(double x) {
- long long i, r;
- double x2=x*0.5, y=x;
- i = *(long long*)&y;
- i = 0x5fe6eb50c7b537a9 - (i>>1);
- y = *(double*)&i;
- for(r=0 ; r<10 ; r++) y = y * (1.5 - (x2*y*y));
- return x * y;
- }
- int main() {
- double n;
- while(1) {
- scanf("%lf", &n);
- printf("sqrt = %.10lf\n", sqrt(n));
- }
- return 0;
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement