Help with if statement

Jikosan jikosan at myrealbox.com
Sun Jan 11 10:49:41 CET 2004


I got it working. I had originally misread the prompt.

Here's the code I finally got working. It calculates PI using the Monte
Carlo method.
I have to graph log N vs log PI, which explains the last few lines of the
program.

Gene

#HW1 - 1.5.2
#This code was written in the Python programming language.

import random   #Loads the random number generator library.
import math     #Loads the math library.

file=open('output.txt', 'w')    #Identifies output file.
a, b = -1.0, 1.0                #Sets the boundary conditions for (x,y).
for N in range(0, 100000, 25):  #4000 trials, 1<= N <= 100000 at steps of
25.
    if N == 0:
        N = 1
    onecheck = 0
    for n in range(0,N):        #This for loop generates (x,y) sets N times.
        xrand = random.random()
        yrand = random.random()
        x = a * (1.0-(xrand)) + b*(xrand)
        y = a * (1.0-(yrand)) + b*(yrand)
        c = math.sqrt(math.pow(x,2.0) + math.pow(y,2.0))    #Calculates the
distance to origin.
        if c < 1:
            onecheck = onecheck+1
    pi = 4.0*(onecheck*math.pow(N,-1))  #This is to calculate P(N), the
ratio of points whose
                                        #distance is less than 1 to the
total number of trials.
    if pi == 0.0:   #Prevent 'log 0' calculation errors
        pi == 1.0

#Convert pi and N to log10 and write to a file to graph on Excel.
    file.write(str(math.log(N, 10)))
    file.write('    ')
    file.write(str(math.log(pi, 10)))
    file.write('\n')
file.close()





More information about the Python-list mailing list