This is the function in question. It calculates the Pearson correlation coefficient for p1 and p2, which is supposed to be a number between -1 and 1.
When I use this
Integer division is confusing it. It works if you make n a float:
n
n=float(len(si))