diff --git a/dice.py b/dice.py
new file mode 100644
index 0000000..40ffc6b
--- /dev/null
+++ b/dice.py
@@ -0,0 +1,76 @@
+import random
+import math 
+
+number_of_side1  =  int(input("How many sides on dice 1? "))
+number_of_side2  =  int(input("How many sides on dice 2? "))
+number_of_side3  =  int(input("How many sides on dice 3? "))
+number_of_side4  =  int(input("How many sides on dice 4? "))
+add = int(input("minimum: 0/1 "))
+
+add1=add*min([1,number_of_side1])
+add2=add*min([1,number_of_side2])
+add3=add*min([1,number_of_side3])
+add4=add*min([1,number_of_side4])
+
+print()
+print("Input rolling method")
+print("0: single dice")
+print("1: highest of four")
+print("2: lowest of four")
+print("3: sum of four")
+print("4: difference of four")
+print("5: product of four")
+method=int(input())
+print()
+exp=int(input("Input number of dice rolls 10^"))
+
+results = []
+sum_of_results = 0.0
+maxi=0
+result = 0
+
+trials = 0
+while trials < 10**exp:
+	if method==1:
+		result = max([int(random.random()*number_of_side1)+add1,int(random.random()*number_of_side2)+add2,int(random.random()*number_of_side3)+add3,int(random.random()*number_of_side4)+add4])
+	elif method==2:
+		result = min([int(random.random()*number_of_side1)+add1,int(random.random()*number_of_side2)+add2,int(random.random()*number_of_side3)+add3,int(random.random()*number_of_side4)+add4])
+	elif method==3:
+		result = (int(random.random()*number_of_side1)+add1)+(int(random.random()*number_of_side2)+add2)+(int(random.random()*number_of_side3)+add3)+(int(random.random()*number_of_side4)+add4)
+	elif method==4:
+		result = abs((int(random.random()*number_of_side1)+add1)-(int(random.random()*number_of_side2)+add2)-(int(random.random()*number_of_side3)+add3)-(int(random.random()*number_of_side4)+add4))
+	elif method==5:
+		result = (int(random.random()*number_of_side1)+add1)*(int(random.random()*number_of_side2)+add2)*(int(random.random()*number_of_side3)+add3)*(int(random.random()*number_of_side4)+add4)
+	else:
+		result = int(random.random()*number_of_side1)+add1
+	
+	maxi = max(maxi,result)
+	sum_of_results += result
+	results.append(result)
+	trials += 1
+
+print(f"The average result is about {sum_of_results/trials}")
+
+ordered = [0]*(maxi+1)
+
+for result in results:
+	ordered[result] += 1
+
+maxi = 0
+result = 0	
+for number in ordered:
+	print(f"{result:03d} occured {number} times with a probability of {number*100/trials}%")
+	maxi = max(maxi, number)
+	result += 1
+
+result = 0	
+for number in ordered:
+	boxes = int(number*50/maxi)
+	box = 0
+	str = f"{result:03d}: "
+	while box < boxes:
+		str += "#"
+		box += 1
+	
+	print(str)
+	result += 1
\ No newline at end of file