/* parallel Monte Carlo area calculation */
/* slave program for third quadrant */

#include <stdio.h>
#include <math.h>
#include <pvm3.h>
#define steps 50000			/* number of points */
#define xmin -0.5
#define xmax 0.0
#define ymin -0.5
#define ymax 0.0

double f(double x)			/* the function for the */
{					/* border in this quadrant */
   return(-sqrt((1-x*x/0.25)*0.09));
}

main()
{
int i, ptid, bufid, dum, send[1], type;
double x, y;

srand48(pvm_mytid());			/* seed the number generator */
					/* not a good way of doing it */
ptid = pvm_parent();


do
{
   send[0]=0;
   for (i=1; i < =steps; i++)
   {
      x= drand48()*(xmax+xmin);		/* random points in the */
      y= drand48()*(ymax+ymin);		/* quadrant */
      if (f(x) <= y) send[0]++;     		/* point is inside the figure */
   }     
   pvm_initsend(PvmDataDefault);		/* send result back to master */
   pvm_pkint(send, 1, 1); 
   pvm_send(ptid, 3);
   bufid=pvm_recv(ptid, -1);			/* any message from master */
   pvm_bufinfo(bufid, &dum, &type, &dum);	/* more work ? */
}while(type);

pvm_exit;  					/* type=0 means we are done */ 
}