Find PI using Monte Carlo method We implement three version of this method: Classic and simple iteration of all samples Divide samples in N piece and use N routine to compute Divide samples in NumCPU piece and use NumCPU thread to compute