/*** integrate.c ***/
#include <stdio.h>
#include <stdlib.h>
#include <math.h>


int main(int argc, char *argv[])
{
  double x, y;                                                    /* x, y */
  double sum;                                         /* hits below curve */
  int t;                                                  /* loop counter */
  int num_runs;                     /* number of generated random numbers */
  double max_x, max_y;               /* maxima of generated random values */

  num_runs = atoi(argv[1]);                            /* read parameters */
  sscanf(argv[2], "%lf", &max_y);

  sum = 0.0;
  max_x = 2;
  for(t=0; t<num_runs; t++)
  {
    x = max_x*drand48();           /* generate random number in [0,max_x] */ 
    y = max_y*drand48();           /* generate random number in [0,max_y] */
    if(y < -x*x+2*x+1)                       /* is y number below curve ? */
	sum += 1;          
  }

  printf("%e %f\n", (double) num_runs,                    /* print result */
	 max_x*max_y*sum/num_runs);

  return(0);
}