1 2#include <stdio.h> 3#include <stdlib.h> 4#include <math.h> 5 6int main ( void ) 7{ 8 int i, j; 9 double s, r; 10 s = 0.0; 11 double* a1 = malloc(1000 * sizeof(double)); 12 double* a2 = malloc(1000 * sizeof(double)); 13 for (i = 0; i < 1000; i++) { 14 a1[i] = s; 15 s += 0.3374893482232; 16 a2[i] = s; 17 } 18 19 s = 0.0; 20 r = 0.0; 21 for (j = 0; j < 5000; j++) { 22 for (i = 0; i < 1000; i++) { 23 s += (a1[i] - a2[i]) * (a1[i] + a2[i]) - sqrt(r + 1.0); 24 r += 0.001; 25 } 26 } 27 printf("s = %f, r = %f\n", s, r ); 28 return 0; 29} 30