This is a simple synthesizer based on NumPy, which allows you to create surprisingly complex sounds with surprisingly simple code. Uses the Python library `sounddevice` for playback, so you'll need to `pip install sounddevice`

sample.py

import time
import spatial as sp
import sounddevice as sd

def play(sound, sample_rate=sp.SAMPLE_RATE):
  """Helper to play serially"""
  sd.play(sound)
  time.sleep(len(sound) / sample_rate)

space = sp.space(3)
play(sp.sin(space, 440)) #concert A sin wave for three seconds
play(sp.sin(space, 440) * (sp.sin(space, 1.) / 2. + 0.5)) #beating w/ a sin wave
play(sp.sigmoid(sp.sin(space,60), sp.sin(space,60)+220)*0.5) #buzzing by modulating frequency
play(sp.sigmoid(sp.sin(space,60), sp.sin(space,60)+220)*(sp.sin(space, 2.) / 2. + 0.5)*0.5) #buzzing by modulating frequency, with a beat

spatial.py

import numpy as np

SAMPLE_RATE=44100

def smoothstep(edge0, edge1, x):
  """broadcastable hermite interpolation, e.g. a sigmoid function"""
  t = np.clip((x - edge0) / (edge1 - edge0), 0.0, 1.0)
  return t * t * (3.0 - 2.0 * t)

def space(duration, sample_rate=SAMPLE_RATE):
  """Create a samplespace of `duration` seconds"""
  return np.mgrid[0:int(sample_rate*duration)] / sample_rate

def sin(space, freq,shift=0):
  """given an arbitrary space, modulate it using a sin wave of `freq` frequency, shifted by `shift` distance"""
  return np.sin(2*np.pi*space*freq+shift)

def sigmoid(space, freq, shift=0):
  """given an arbitrary space, modulate it using a sigmoid wave of `freq` frequency, shifted by `shift` distance"""
  return smoothstep(0, 1., sin(space, freq, shift)) + smoothstep(-1., 0., sin(space, freq, shift))