#include "real.h"
#include <stdlib.h>

void plainmc(int n,real a[],real b[],real (*func)(real*),
	int N,real *sum,real *sum2,void (*rnd)(int,real*,real*,real*))
{
int i; real f,*x; x=(real*)malloc(n*sizeof(real));
*sum=0; *sum2=0;
for(i=1;i<=N;i++){(*rnd)(n,a,b,x); f=(*func)(x); *sum+=f; *sum2+=f*f;}
}