from perceptivo.psychophys.oracle import reference_audiogram
from perceptivo.psychophys.model import Gaussian_Process
from perceptivo.types.psychophys import Sample

oracle = reference_audiogram(scale=3)
model = Gaussian_Process(amplitude_range=(5,35))

for i in range(100):
    sound = model.next()
    sample = Sample(response=oracle(sound), sound=sound)
    model.update(sample)

model.plot()