#!/usr/bin/env python

import math
import random

def mean():
	t, alist, mu = 0, [], 0
	for i in range(1,501):
		alist.append(random.randint(1,501))
		t=t+i
	mu  = (t / len(alist))
	#STD(mu,alist)

def STD():
	# this is the mean section
	t, alist, mu = 0, [], 0
	for i in range(1,501):
		alist.append(random.randint(1,501))
		t=t+i
	mu  = (t / len(alist))
	print "This computed Mean is: " + str(mu)
	
	# this is the STD portion
	ans, t, std = 0, 0, 0
	for i in alist:
		t = t + math.pow(i - mu, 2)
	ans = (t / (len(alist)-1))
	std = math.sqrt(ans)
	print "This computed Standard Deviation is: " + str(std)
	
	# this is test for skewness
	t, skew = 0, 0
	for i in alist:
		t = t + math.pow(i - mu, 3)
	skew = ( t / ( ( len(alist) - 1) * math.pow(std,3)) )
	print "This computed Coefficient of Skewness measure of symmetry is: " + str(skew)
	
	# this is test for Kurtosis
	t, kurt = 0, 0
	for i in alist:
		t = t + math.pow(i - mu, 4)
	kurt = ( t / ( ( len(alist) - 1) * math.pow(std,4)) )
	print "This computed Coefficient of Kurtosis measure of bight to STD: " + str(kurt)
	
	# this is Coefficient of Excessive Kurtosis compared with NOrmal
	print "Coefficient of Excessive Kurtosis compared with Normal: " + str(kurt - 3)
	

def main():
	mean()
	STD()

# added proper main function
if __name__ == "__main__":
	main()