- #1
alex steve
- 5
- 0
I am having an issue with finding the entropy in my program. I was asked to the find the entropy of the distribution as a function of time but i do not know where to start with entropy.
I understand entropy but putting it in my program is where I am stuck
Here is my code:
I understand entropy but putting it in my program is where I am stuck
Here is my code:
Python:
# -*- coding: utf-8 -*-
"""
Created on Thu Nov 12 11:15:44 2015"""import matplotlib.pyplot as plt
import random
def Function(D): #D = dimensions
sumOfSquare = 0.0
for i in range(0, len(D)):
sumOfSquare += D[i]**2
if sumOfSquare <=1:
return 1
else:
return 0
def MonteCarlo(f_n,dim):
intervalsForSphere = 1000000
integral = 0.0
for i in range(0, intervalsForSphere):
for j in range(0,len(dim)):
dim[j] = random.random()
integral += f_n(dim)
integral = (1-(-1))**len(dim)/intervalsForSphere * integral
return integral
print("10 dimensional unit circle ")
Ten_Dim= list(range(1,10+1))
ten_D_circle = MonteCarlo(Function,Ten_Dim)
print("area:",ten_D_circle)
AreaofCircle = []
x = []
for i in range(1,13):
D = list(range(1,i+1))
AreaofCircle.append(MonteCarlo(Function, D))
x.append(i)
plt.plot(x,AreaofCircle)
plt.xlim([0,13])
plt.xlabel("Dimensions")
plt.ylabel("area")
plt.title("Area of N-dimensional Unit Circle")
plt.show()