### Print a help message
println("Plot data for the photonreceptor model using:")
println(" pyplot_photonsequence()")
println(" pyplot_bump(timepoints,amplitude = 4.0,shape = 3.0,scale =2.5)")

### Plot the photon sequence data using the PyPlot package
function pyplot_photonsequence()

    #calculate the data
    photons = photonsequence(10:100,1000)
    timepoints = dataindex(photons)

    #plot the data
    fig = PyPlot.figure("functionality/photoreceptor")
    PyPlot.plot(timepoints,datavalues(photons))
    PyPlot.xlabel("Time (s)")
    PyPlot.ylabel("Photons")
    PyPlot.title("Photon Sequence")
    PyPlot.grid("on")

    #return the figure object
    fig
end

function pyplot_bump(timepoints,amplitude = 4.0,shape = 3.0,scale = 2.5)

    #calculate the data
    bmp = bump(collect(timepoints),amplitude,shape,scale)

    #plot the data
    fig = PyPlot.figure("functionality/bump")
    PyPlot.plot(timepoints/1000.0,bmp)
    PyPlot.xlabel("Time (s)")
    PyPlot.ylabel("Current (nA)")
    PyPlot.title("Photon-Induced Current")
    PyPlot.grid("on")
end

pyplot_photonsequence()
pyplot_bump(1.0:1.0:30.0)

# File return statement
nothing
