def test_plot_bounded_coal(self):
n = 1000
k = 6
T = 500
# plots should differ
p = plotfunc(lambda t: coal.prob_bounded_coal(t, k, n, T),
0, 1000, 10)
p.plotfunc(lambda t: coal.prob_coal(t, k, n),
0, 1000, 10)
def bspline2(pts, times, t):
p = len(times) - len(pts) - 1
b = []
for i in range(len(pts)):
b.append(bbasis(i, p, t, times))
return b
if __name__ == "__main__":
f = lambda x: hermite(0, 0, 1, -1, x)
g = lambda x: bezier(0, 1, 1, 0, x)
h = lambda x: catmullrom(-3, 0, 0, -3, x)
if False:
p = gnuplot.plotfunc(f, 0, 1, .1)
p.plotfunc(g, 0, 1, .1)
p.plotfunc(h, 0, 1, .1)
if False:
pts = [1, 2, 2, 3, 3, 1, -5, -4, -5, -3, -5, -2, -5, 0]
p = gnuplot.plotfunc(interpolate(pts), 0, len(pts)-1, .1)
def frange(n):
return map(lambda x: x/float(n), range(n+1))
pts = [0, 4, 0]
pts2 = [0, 4, 6]
p = gnuplot.Gnuplot()
