new_contents: import scipy.stats as stats
  import numpy as np
  
  def check_sequential(array):
  	'''
  	Input: list of numbers
  	Output: True if the numbers do not correspond to a normal distribution with p < 0.05
  	False otherwise
  	'''
  	zval, pval = stats.normaltest(array)
  	if pval < 0.05:
  		return True
  	else:
  		return False
  
  def pdf(input, array):
  	'''
  	Input: new value, list of numbers
  	Output: pdf * 100 that new value belongs to this list of numbers assuming the list is
  	modelled as a normal distribution
  	'''
  	return (stats.norm.pdf(input, loc=np.mean(array), scale=np.std(array)) * 100)
  	
  
subject: Create number_sanitize.py