path: root/cyclone/sickle/curve.gp

h(p) = (((p + 1e-20) * 1.2) ** .41) * .91
f(p) = h(p) / (1 - h(p))
g(x, p) = (exp(f(p) * x) - 1) / (exp(f(p)) - 1)
splot [0:1] [0:1] g(x, y)
pause -1
plot [0:1] g(.25, x), g(.5, x), g(.75, x)
pause -1
plot [0:1] g(x, 0), g(x, .25), g(x, .5), g(x, .75), g(x, 1)
print g(.5, 0)