#include <stdio.h>
#include <stdlib.h>
#include <math.h>

main() {
double x, func = 0.0, func2 = 0.0, resultF, errorF, resultErr;
int try, ntry;
printf("Input the number of MC trials\n");
scanf("%d",&ntry);
srand((unsigned)time((long *)0));
for (try=0; try<ntry; try++) {
x = rand()/(float)RAND_MAX;
func = func + x*x; /* function to be evaluated is parabola of x^2 */
func2 = func2 + x*x*x*x;
errorF = sqrt(func - func2);
}
resultF = func/ntry;
resultErr = errorF/ntry;
printf("MC estimate for function = %f\n", resultF);
printf("MC error for function = %f\n", resultErr);
}