def gaussian(x, mu, sig):
    return (np.exp(-((mu - x)*(x - mu))/(2*sig*sig)))/(np.sqrt(2*(np.pi)*sig*sig))

#=======================================================================================================================
# Start program
#=======================================================================================================================
dataset/ETHPy150Open neuropoly/spinalcordtoolbox/dev/spinal_level/sct_extract_spinal_levels.py/gaussian/swapped operands of `-`