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