import numpy as np
import matplotlib.pyplot as plt
 
def scatterFile(file, what, color):
	d = np.loadtxt(file, skiprows=1, delimiter=',')
 
	# arrange data
	sec = []
	vals = []
	for i in range(0, len(d)):
		# take a single second
		sec.append(d[i][0])
		# take the avarage of values for that second
		vals.append( np.average([d[i,j] for j in range(1,201)]) )
	print(sec)
	print(vals)
 
	# show data on plot
	plt.scatter(sec,vals,'r')
 
 
# Main:
def main():
	title = 'Concentration of g_i (red) and r_i (blue)'
	plt.title(title)
	plt.xlabel('time (s)')
 
	file  = 'g_internal.csv'#'prova.csv'
	what  = 'g_i'
	color = 'r'
	scatterFile(file,what,color)
 
	file  = 'r_internal.csv'#'prova.csv'
	what  = 'r_i'
	color = 'b'
	scatterFile(file,what,color)
	plt.show()
	print("done")

aW1wb3J0IG51bXB5IGFzIG5wCmltcG9ydCBtYXRwbG90bGliLnB5cGxvdCBhcyBwbHQKCmRlZiBzY2F0dGVyRmlsZShmaWxlLCB3aGF0LCBjb2xvcik6CglkID0gbnAubG9hZHR4dChmaWxlLCBza2lwcm93cz0xLCBkZWxpbWl0ZXI9JywnKQoKCSMgYXJyYW5nZSBkYXRhCglzZWMgPSBbXQoJdmFscyA9IFtdCglmb3IgaSBpbiByYW5nZSgwLCBsZW4oZCkpOgoJCSMgdGFrZSBhIHNpbmdsZSBzZWNvbmQKCQlzZWMuYXBwZW5kKGRbaV1bMF0pCgkJIyB0YWtlIHRoZSBhdmFyYWdlIG9mIHZhbHVlcyBmb3IgdGhhdCBzZWNvbmQKCQl2YWxzLmFwcGVuZCggbnAuYXZlcmFnZShbZFtpLGpdIGZvciBqIGluIHJhbmdlKDEsMjAxKV0pICkKCXByaW50KHNlYykKCXByaW50KHZhbHMpCgoJIyBzaG93IGRhdGEgb24gcGxvdAoJcGx0LnNjYXR0ZXIoc2VjLHZhbHMsJ3InKQoKCiMgTWFpbjoKZGVmIG1haW4oKToKCXRpdGxlID0gJ0NvbmNlbnRyYXRpb24gb2YgZ19pIChyZWQpIGFuZCByX2kgKGJsdWUpJwoJcGx0LnRpdGxlKHRpdGxlKQoJcGx0LnhsYWJlbCgndGltZSAocyknKQoKCWZpbGUgID0gJ2dfaW50ZXJuYWwuY3N2JyMncHJvdmEuY3N2JwoJd2hhdCAgPSAnZ19pJwoJY29sb3IgPSAncicKCXNjYXR0ZXJGaWxlKGZpbGUsd2hhdCxjb2xvcikKCglmaWxlICA9ICdyX2ludGVybmFsLmNzdicjJ3Byb3ZhLmNzdicKCXdoYXQgID0gJ3JfaScKCWNvbG9yID0gJ2InCglzY2F0dGVyRmlsZShmaWxlLHdoYXQsY29sb3IpCglwbHQuc2hvdygpCglwcmludCgiZG9uZSIp